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TEMPERATURE  AND  VELOCITY  MEASUREMENTS 
AND  PREDICTIVE  MODEL  COMPARISONS 
IN  THE  NEAR- FIELD  REGION 
OF  SURFACE  THERMAL  DISCHARGES 


by 

R.  A.  Paddock,  A.  J.  Policastro, 
A.  A.  Frigo,  D.  E.  Frye, 
and  J.  V.  Tokar 


ABSTRACT 

Simultaneous  temperature  and  velocity  measurements 
were  made  in  the  near-field  regionof  the  surface  thermal  dis¬ 
charge  at  the  Point  Beach  Unit  1  and  Palisades  Nuclear  Power 
Plants  on  Lake  Michigan.  Data  collected  include  measurements 
of  temperature  and  velocity  at  the  0.5-,  1.0-,  1.5-,  2.0-,  2.5-, 
and  3.0-m  depths ,  along  with  measurements  of  ambient  lake  and 
meteorological  conditions.  Bottom  depth  was  also  measured  at 
various  locations.  Four  such  surveys  were  made  at  the  Point 
Beach  plant,  three  at  the  Palisades  plant. 

To  examine  the  jet  features  from  the  above  surveys  and 
facilitate  a  comparison  with  analytical  model  predictions,  a 
computer  program  was  written  to  smooth  the  data,  extracting 
such  jet  characteristics  as  trajectory,  centerline  temperature 
decay  and  temperature  half-widths,  centerline  velocity  decay 
and  velocity  half- widths ,  temperature  and  velocity  half-depths, 
and  isotherm  areas. 

Four  near- field  analytical  models  oftenused  in  environ¬ 
mental  impact  evaluations  of  power-plant  surface  discharges 
are  compared  to  the  jet  characteristics  determined  from  the 
smoothed  jet  data.  The  Pritchard  model  compares  rather  well 
with  these  limited  data  and  is  often  conservative  when  model-data 
discrepancies  exist.  The  Stolzenbach- Ha  rleman  and  Prych 
models  predict  too  rapid  a  temperature  and  velocity  decay  ac¬ 
companying  too  great  a  lateral  spread.  The  Motz- Benedict 
model  is  too  sensitive  to  an  entrainment  coefficient  E  with  lit¬ 
tle  consistent  data  available  for  its  determination  for  accurate 
prediction. 

Recommendations  for  future  research,  encompassing 
the  field-data  acquisition,  the  smoothing  procedure,  and  the 
presently  available  models  are  included. 


I.  INTRODUCTION 


Under  the  auspices  of  the  U.S.  Atomic  Energy  Commission,  the  Argonne 
Center  for  Environmental  Studies  has  been  studying  the  physical  effects  of 
heated  condenser  discharges  from  steam-electric  power  plants  on  the  Great  Lakes 
since  FY  1970.  Appendix  A  lists  the  reports  published  under  this  program  to 
the  present.  Two  of  the  primary  objectives  of  this  program  have  been  and 
continue  to  be  the  acquisition  of  reasonably  complete  prototype  the rmal- plume 
field  data  and  the  verification  of  analytical  predictive  plume  models.  To  this 
end,  field  data  in  the  jet  regime  of  two  nuclear  power  plants  on  Lake  Michigan 
have  been  collected  and  compared  to  models,  with  the  results  reported  herein. 

The  jet  regime  (the  near  field)  is  that  region  of  the  discharge  in  which 
the  heated  effluent  enters  the  receiving  body  of  water  possessing  a  velocity 
and  temperature  disparity  with  respect  to  the  receiving  body.  Thus,  as  a 
heated  effluent  enters  an  ambient  environment  from  a  particular  plant  outfall, 
viscous  shear  between  the  effluent  and  the  ambient  fluid  creates  turbulence 
in  the  contact  region.  This  turbulence  works  its  way  both  inward  toward  the 
jet  centerline  and  outward  toward  the  ambient  fluid,  with  a  resultant  net  out¬ 
ward  flux  of  momentum  and  heat  away  from  the  jet  axis.  Within  this  regime 
of  the  discharge,  it  is  the  mechanical  mixing  action  induced  by  the  kinetic 
energy  of  the  discharge  itself  that  dominates  the  eddy  transport  mechanisms. 

At  some  distance  from  the  outfall,  the  kinetic  energy  of  the  discharge  will  be 
sufficiently  dissipated  to  allow  the  natural  turbulence  existing  within  the  am¬ 
bient  receiving  water,  together  with  buoyant  forces,  to  dictate  plume  disper¬ 
sion.  It  is  nominally  assumed  that  the  effluent  is  no  longer  jetlike  in  character 
when  this  situation  is  reached. 

From  a  regulatory  point  of  view,  the  jet  regime  is  of  particular  in¬ 
terest.  It  is  often  within  this  region  that  outfall  architects  must  design  their 
discharges  to  meet  thermal  wate  r- quality  criteria  that  limit  the  temperature 
rise  in  the  thermal  plume  beyond  a  prescribed  distance  from  the  point  of 
discharge.  This  is  commonly  referred  to  as  a  mixing- zone  limitation.  Some 
states  have  adopted  very  restrictive  mixing- zone  criteria;  others  have  none 
at  all.  Therefore,  depending  on  the  nature  of  the  receiving  body,  the  size  of 
plant,  and  a  multitude  of  different  factors  including  the  thermal  criteria,  each 
plant  outfall  design  is  more  or  less  tailored  to  the  particular  siting  situation. 

On  the  Great  Lakes,  the  predominant  outfall  design  happens  to  be  a  shoreline, 
open,  rectangular  discharge  canal.  Several  more  recent  plant  designs  have 
used  more  sophisticated  offshore  multiorifice  submerged  discharges,  called 
diffuse  rs . 

The  literature  contains  numerous  models  that  attempt  to  predict  the 
behavior  of  shoreline  canal  discharges.  Some  of  these  models  are  qualitative 
in  nature;  others  profess  to  be  quantitative  as  well.  One  thing  all  these  models 
have  in  common  is  that  none  has  been  generally  verified  with  prototype  field 
data.  To  compound  the  problem,  a  survey  of  the  literature  reveals  surprisingly 
little  actual  jet- regime  field  data  with  which  models  can  be  tested.  Since  there 
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is  such  a  paucity  of  data  in  the  jet  regime,  a  field  program  was  developed  by 
the  Center  for  Environmental  Studies  specifically  to  acquire  prototype  data 
near  canal-type  discharges.  This  program  has  been  partially  described  in 
Refs.  21  and  25  of  Appendix  A.  A  description  of  experimental  methods  and 
detailed  information  concerning  the  results  of  seven  jet  studies  obtained  during 
1972  are  presented  herein. 

Since  it  would  have  been  difficult  to  compare  the  field  data  directly  to 
results  obtained  from  analytical  models,  a  data- smoothing  technique  was  de¬ 
veloped  to  help  in  this  endeavor.  The  smoothing  method  was  primarily  devel¬ 
oped  to  glean  as  much  information  from  the  experimental  data  as  possible, 
considering  the  limited  number  of  data  points  collected  using  the  present  field 
technique.  A  complete  description  of  the  smoothing  method  appears  in  Sec.  V. 
While  one  should  recognize  that  the  smoothing  method  has  some  obvious  limi¬ 
tations  and  biases,  it  nevertheless  has  worked  out  quite  well  for  the  purposes 
for  which  it  was  designed. 

Lastly,  the  results  of  the  smoothing  procedure  are  compared  to  four  ana¬ 
lytical  models  that  have  been  used,  in  some  cases  quite  extensively,  for  pre¬ 
dictive  purposes.  These  comparisons  and  a  discussion  of  them  appear  in 
Sec.  VII.  Note  that  the  success  or  apparent  lack  of  success  displayed  by  a 
particular  model  should  not,  at  this  point,  be  considered  as  a  total  test  of  the 
model.  Many  more  data  comparisons  under  different  outfall  situations  must 
be  made  before  any  model  can  be  realistically  evaluated. 

In  summary,  this  report  brings  together  details  concerning  the  acqui¬ 
sition,  smoothing,  and  model  analysis  of  a  relatively  unique  set  of  jet- regime 
plume  field  data.  We  hope  this  report  will  stimulate  more  interest  than  has 
been  shown  in  the  literature  in  attempting  to  validate  existing  predictive 
models.  In  our  opinion,  too  many  predictive  models  existing  in  the  literature 
have  not  been  adequately  tested.  Much  reliable  field  data  is  just  now  becom¬ 
ing  available,  and  it  should  be  the  immediate  goal  of  those  interested  in  ap¬ 
plying  predictive  models  to  test  these  models  with  actual  field  data.  Only  in 
this  way  can  we  hope  for  a  positive  improvement  in  the  existing  state  of  the 
art. 


II.  EXPERIMENTAL  TECHNIQUE 


A  two-  or  three-point  mooring 
-|-m  ( 1 8  -  ft)  cathedral-hull  fiberglass 


Fig.  1.  Boat  with  Current  Meter  Used 
in  Jet-regime  Studies 


system  was  used  to  hold  Argonne' s 
boat,  the  R.  V.  Aha,  steady  while  ob¬ 
taining  simultaneous  temperature  and 
velocity  measurements  in  the  near¬ 
field  region  of  the  thermal  plume  (see 
Fig.  l).  Anchors  were  located  on 
either  side  of  the  plume,  and  for  re¬ 
gions  very  near  the  outfall,  a  third  line 
was  sometimes  attached  to  the  outfall 
itself.  Transects  across  the  plume 
centerline  were  then  made  at  various 
distances  from  the  outfall.  The  posi¬ 
tion  of  the  boat  was  held  relatively 
constant  at  the  various  measuring 
stations,  and  the  position  of  each 
station  was  determined  by  using  a 
Motorola  Mini- Ranger  range  position¬ 
ing  system  (Fig.  2).  This  positioning 
system  consists  of  two  shore-based 
transponders  with  a  receiver-transmitter 
unit  and  range  console  on  board  the 
boat  that  displays  the  range  information 


Fig.  2.  Motorola  Mini-Ranger  Range  Positioning  System 
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from  each  transponder.  The  boat's  position  can  then  be  found  by  trilate ration. 
The  Mini- Ranger  is  powered  by  110  V  ac  (available  from  a  24- V  dc  high- 

efficiency  Flitetronics  PC  16  Air¬ 
craft  Static  Inverter).  The  usable 
range  of  the  system  with  omnidirec¬ 
tional  antennas  is  about  16  km. 

A  Bendix  Q- 1  5  geomagnetic 
ducted  current  meter  (Fig.  3),  with 
an  attached  YSI  thermistor,  was  used 
to  measure  the  velocity  and  tempera¬ 
ture  of  the  discharge  waters.  The 
Q- 1  5  has  a  five-bladed  impeller  that 
rotates  in  both  directions  and  is  en¬ 
closed  in  a  duct.  The  ducted  assembly 
is  aligned  with  the  current  by  a  vane 
of  adjustable  length.  The  effects  of 
wave  and  boat  motion  are  nulled  out 
by  electronic  averaging  (over  about 
2  5  sec)  of  the  number  of  turns  of  the 
impeller  and  by  the  presence  of  the 
duct.  Current  speed  and  direction 
are  displayed  on  deck  by  means  of  a 
readout  unit,  Bendix  Model  No.  S-232, 
which  is  connected  to  the  current 
meter  through  a  four- conducto r  cable. 
The  current  meter  is  powered  by 
six  9-  V  batteries.  The  meter  was 
lowered  over  the  side  of  the  boat  and 
suspended  at  0.5-m  intervals  to  a 
depth  of  3.0  m  or  to  the  bottom.  The 
first  time  the  experiment  was  being 
conducted,  it  was  discovered  that  time  variations  in  velocity  and  temperature 
occur.  Thus,  in  order  to  obtain  average  values  of  velocity  and  temperature 
along  with  any  variation,  strip  chart  recorders  were  connected  to  the  current- 
meter  and  thermistor  outputs. 


Fig.  3.  Bendix  Q-15  Geomagnetic  Ducted  Cur¬ 
rent  Meter.  ANL  Neg.  No.  190-568-11. 


Note  at  this  point  that  a  variety  of  factors  inherent  in  making  measure¬ 
ments  in  the  jet  regime  may  cause  uncertainties  in  the  data.  These  problems 
must  be  understood  if  proper  use  is  to  be  made  of  the  data. 

An  important  aspect  of  the  experimental  uncertainty  is  the  short-term 
variations  in  the  velocity  and  temperature  of  the  discharge  jet  (see  Ref.  30). 

The  cause  of  these  variations  is  not  clear  at  present,  but  they  may  be  due  to 
eddies  created  at  the  interface  between  the  jet  and  the  ambient  water,  to  surging 
which  is  apparent  in  the  discharge  canal,  or  to  other  factors.  These  fluctua¬ 
tions  appear  to  have  periods  ranging  from  a  few  seconds  to  a  few  minutes. 
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Since  point-by-point  measurements  in  the  jet  were  typically  made  over  a 
period  of  1-2  min,  it  is  apparent  that  an  unrepresentative  value  might  be  ob¬ 
tained  at  a  given  point.  The  scale  of  the  short-term  temperature  fluctuations 
is  on  the  order  of  several  Centigrade  degrees  or  less;  the  velocity  variations 
are  on  the  order  of  50%  of  the  mean  value  or  less.  These  fluctuations  were 
not  present  at  all  locations.  A  more  typical  value  for  the  uncertainty  in  the 
temperature  measurement  is  ±0.5C°;  a  typical  value  for  the  uncertainty  in  the 
current  speed  is  ±20%. 

Another  source  of  uncertainty,  in  terms  of  data  analysis,  is  the  ambi¬ 
guity  attached  to  the  values  of  ambient  current  and  temperature.  These  num¬ 
bers,  necessarily  assumed  to  be  constants  throughout  the  measurement  (which 
lasted  from  3  to  7  hr),  vary  not  only  in  time  but  in  space  as  well.  Ambient- 
current  measurements  were  typically  made  at  a  single  location  (at  several 
depths)  before  and  after  the  jet- regime  measurements.  Here  again,  current 
fluctuations  in  time  and  position  may  lead,  for  a  variety  of  reasons,  to  an  un¬ 
representative  value  for  ambient- cur  rent  speed.  (Direction  of  the  current  is 
thought  to  be  more  definite.)  Ambient- cur  rent  speeds  as  reported  may  have 
an  uncertainty  of  20-  50%;  lower  current  speeds  are  the  most  uncertain.  On 
some  occasions,  ambient-temperature  measurements  are  as  difficult  to  pin 
down  as  ambient- cur  rent  speed  and  are  somewhat  more  important  in  terms  of 
the  analysis  to  be  described.  We  chose  the  ambient  temperature  to  be  the 
water  temperature  (at  the  appropriate  depth),  which  appeared  not  to  be  in¬ 
fluenced  by  the  discharge  water,  yet  was  in  the  vicinity  of  the  discharge.  Un¬ 
fortunately,  on  days  when  upwelling,  downwelling,  shoreline  heating,  or  other 
disturbing  phenomena  occurred,  the  reported  values  of  ambient  water  tempera¬ 
ture  may  have  an  uncertainty  of  as  much  as  1-2C°. 

In  the  face  of  the  previously  discussed  uncertainties  in  the  data,  in¬ 
strumentation  accuracy  may  not  be  very  important,  but  for  completeness, 

Table  I  lists  the  instrument  specifications.  Of  special  interest  in  terms  of 
velocity  measurements  is  the  fact  that,  while  a  ducted  current  meter  is  used 
to  null  out  disturbing  vertical  motion,  shielding  of  the  impeller  occurs  if  the 
duct  is  not  aligned  with  the  flow.  The  importance  of  this  remains  unclear,  but 
because  the  meter  continuously  averages  the  speed  over  a  25-sec  period,  any 
shielding  would  result  in  lower  values  for  current  speed. 


TABLE  I.  Instrument  Accuracy 


Instrument 

Sensor 

Accuracy 

Threshold 

Range 

Time  Constant 

Resolution 

Remarks 

Bendix  0-15 
Current  Meter 

Speed:  impeller 

±4%  of 
full  scale 

3  cm/sec 

0-1.0  knot,  low  scale 
0-5.0  knots,  high  scale 

25  sec;  low  scale 

2.5  sec:  high  scale 

1  cm/sec;  low  scale 

5  cm/sec;  high  scale 

Direction;  compass 
with  vane 

±12° 

- 

0-360° 

5° 

Vane  has  adjustable  length.  0.3-3. 0  m 

Temperature 

Recorder9 

Thermistor 

+0.5C° 

0-30°C 

-2.5  sec 

0.2C° 

Consisting  of  a  Rustrak  Model  No.  2133 
temperature  recorder  and  YSI  No.  409 
thermistor  probe 

Temperature 

Recorder^ 

Thermistor 

a0.2C° 

- 

0-50°C 

-2.5  sec 

0.1C° 

Consisting  of  YSI  No.  709  probe  and 
digital  readout  built  at  Argonne 

Motorola 

Mini-Ranger 

+3  m 

0.1-35  km 

1  m 

Accuracy  applies  to  each  range  meas¬ 
urement.  System  accuracy  varies  with 
position  of  transponders  relative  to  boat 
For  most  of  these  measurements.  +3  m 
would  apply. 

aUsed  in  temperature  measurements  up  to  and  including  July  19,  1972. 
bUsed  in  temperature  measurements  after  July  19,  1972. 
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III.  DESCRIPTIONS  OF  POWER  PLANTS 

The  jet  regimes  of  the  thermal  plumes  were  surveyed  at  two  power 
plants  located  on  Lake  Michigan.  The  power  plants  studied  were  the  Point  Beach 
Nuclear  Power  Plant,  operated  by  the  Wisconsin  Electric  Power  Company  and 
the  Wisconsin  Michigan  Power  Company,  and  the  Palisades  Nuclear  Power 
Plant,  operated  by  Consumers  Power  Company.  Brief  descriptions  of  these 
plants  follow. 

A.  Point  Beach  Nuclear  Power  Plant 


The  Point  Beach  Nuclear  Power  Plant  is  in  the  town  of  Two  Creeks, 
Wisconsin,  on  the  western  shore  of  Lake  Michigan.  (Figure  4  is  an  aerial 
view  of  the  plant.)  The  plant  is  a  two-unit  steam- generating  station.  The 
nuclear  reactors  for  each  unit  are  pressurized  light- water-moderated  and 
-cooled  systems.  Each  unit  has  a  gross  capacity  of  523  MWe  and  a  net  ca¬ 
pacity  of  505  MWe.  The  water  intake  for  the  plant  consists  of  a  circular  crib 
533  m  from  the  shore.  Cooling  water  for  the  operation  of  the  power  plant  is 
drawn  from  Lake  Michigan  and  passes  through  the  cooling  condensers  at  a 
maximum  rate  of  25.1  m3/ sec  for  full-power  operation  of  each  individual  unit 
of  the  plant.  The  water  is  returned  to  the  lake  about  50  m  offshore  through 
two  10.7-m-wide  discharge  flumes  (one  flume  per  unit).  Water  depth  in  the 
flumes  is  about  4.2  m.  During  most  of  the  field  year,  the  second  unit  was  not 
operational.  Late  in  the  summer,  however,  the  second  unit  was  operating  at 
about  12%  power  and  50%  of  its  rated  discharge  flow. 

B.  Palisades  Nuclear  Power  Plant 


The  Palisades  Nuclear  Power  Plant  is  near  the  city  of  South  Haven, 
Michigan,  on  the  eastern  shore  of  Lake  Michigan.  (Figure  5  is  an  aerial  view 
of  the  plant.)  This  plant  uses  a  pre s surized- wate r  reactor  to  produce  a  maxi¬ 
mum  gross  output  of  714  MWe.  During  the  1972  field  year,  the  plant  was  oper¬ 
ating  at  a  net  generating  capacity  of  about  420  MWe.  The  cooling  water  is  taken 
from  Lake  Michigan  through  an  intake  crib  located  6.1  m  below  the  lake's  sur¬ 
face,  1.8  m  from  the  lake  bottom,  and  1000  m  from  the  shoreline.  For  the 
Palisades  plant,  as  presently  constructed,  the  cooling  water  passes  through 
the  cooling  condenser  at  a  maximum  flow  rate  of  25.6  m3/ sec  and  is  returned 
to  Lake  Michigan,  via  a  32.9-m-long  discharge  canal  at  the  shoreline.  The 
canal  is  11.3  m  wide  at  the  shoreline  outlet  and  diverges  to  a  width  of  28.3  m 
at  the  point  of  discharge.  At  this  point,  the  water  has  an  average  depth  of 
about  2. 1  m. 
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Fig.  4.  Aerial  View  of  Point  Beach  Nuclear  Power  Plant.  ANL  Neg.  No.  190-499. 
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Fig.  5.  Aerial  View  of  Palisades  Nuclear  Power  Plant 
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IV.  RESULTS  OF  FIELD  MEASUREMENTS 

During  the  1972  field  year,  seven  jet-regime  studies  were  conducted 
at  the  two  power  plants  described  in  Sec.  III.  Specifically,  these  studies  were 
at  the  Point  Beach  Nuclear  Power  Plant  Unit  1  on  May  18,  May  23,  July  13, 
and  September  9,  1972,  and  at  the  Palisades  Nuclear  Power  Plant  on  June  14, 
July  19,  and  October  10,  1972.  In  addition,  one  survey  was  conducted  late  in 
the  1971  field  year  (November  3,  1971)  as  a  preliminary  feasibility  study  of 
the  technique  (see  Appendix  B).  Data  collected  include  measurements  of 
velocity  and  temperature  in  the  near -field  region  of  the  thermal  plume  at  the 
0.5-,  1.0-,  1.5-,  2.0-,  2.5-,  and  3.0-m  depths,  along  with  measurements  of 
ambient  lake  and  meteorological  conditions.  Bottom  depth  was  also  measured 
at  various  locations.  From  the  bottom-depth  data,  approximate  depth  contours 
were  drawn  near  the  outfalls  and  are  shown  in  Figs.  6  and  7.  The  points  indi¬ 
cate  positions  at  which  data  were  taken.*  Results  of  the  jet-regime  measure¬ 
ments  are  shown  in  Figs.  8-48  for  the  dates  indicated.  The  figures  show 
station  locations  at  which  jet  velocities  and  temperatures  were  measured. 

The  velocity  is  represented  vectorially  at  each  station  location.  In  addition, 
the  current  speed,  current  direction,  and  temperature  at  each  station  are 
listed  in  a  table  on  each  figure.  Current  direction  is  given  in  degrees  as 
measured  from  magnetic  north.  Also  listed  on  the  figures  are  the  ambient 
lake  and  meteorological  data,  as  well  as  the  plant  operating  data.  Tempera¬ 
ture  and  velocity  centerlines  and  widths  are  shown  in  most  cases.  The 
widths  represent  the  lateral  distance  from  the  centerline  at  which  the  ap¬ 
propriate  parameter  has  reached  a  value  halfway  between  the  centerline 
value  and  the  ambient  value.  The  mathematical  fitting  technique  used  to 
obtain  the  centerlines  and  widths  is  described  in  Sec.  V.  Centerlines  and 
widths  are  not  shown  for  any  depths  for  the  jet-regime  studies  conducted  at 
the  Palisades  Plant  on  June  14  and  July  19,  1972.  Typically,  the  Palisades 
outfall  produces  a  very  wide  jet,  and  the  data  on  these  dates  did  not  lend 
themselves  to  the  type  of  analysis  necessary  for  determining  centerlines 
and  widths . 


♦Dotted  lines  indicate  estimated  contours  for  which  no  data  were  available. 
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60s 


Fig.  6.  Approximate  Depth  Contours  near  Point  Beach  Outfall.  ANL  Neg.  No.  190-878. 
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Fig.  7.  Approximate  Depth  Contours  near  Palisades  Outfall.  ANL  Neg.  No.  190-880 
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PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAD:  495  MWe 

DISCHARGE  FLOWRATE  24.7  m3/sec 

OUTFALL  TEMPERATURE  17.7°C 

INTAKE  TEMPERATURE  6.7«C 

AMBIENT  WATER  TEMPERATURE  9.2  •  9.6°C 

DRY  BULB  TEMPERATURE:  10.7  -  12.2°C 

RELATIVE  HUMIDITY  75% 

WIND  SPEED  &  DIRECTION  0  -  3.0  m/sec.  110° 
AMBIENT  CURRENT  SPEED  &  DIRECTION  0 
LAKE  SURFACE  CONDITIONS  CALM,  0  -  0.2  m  WAVES 
SKY  CONDITIONS:  CLEAR 


100 


VELOCITY  SCALE  -  cm/sec 


100 
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Fig.  8.  Jet-regime  Study  for  0.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  18,  1972,  1115-1440  Hours.  ANL  Neg.  No.  190-761  Rev.  1. 
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PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAD.  495  MWe 
DISCHARGE  FLOW  RATE  24  7  ni3/sec 
OUTFALL  TEMPERATURE  !7.7°C 
NTAKE  TEMPERATURE  6  7°C 
AMBIENT  WATER  TEMPERATURE  8  5°C 
DRY  BULB  TEMPERATURE  1C  7  -  12  2°C 
RELATIVE  HUMIDITY  75% 

WIND  SPEED  &  DIRECTION  0-3.0  m/sec  110° 
AMBIENT  CURRENT  SPEED  &  DIRECTION  0 
LAKE  SURFACE  CONDITIONS  CALM.  0  -  0.2  m  WAVES 
SKY  CONDITIONS  CLEAR 
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VELOCITY  SCALE  -cm/sec 
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Fig.  9.  Jet-regime  Study  for  l.O-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  18,  1972,  1115-1440  Hours.  ANL  Neg.  No.  190-895. 
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PLANT  LOAD.  495  MWe 
DISCHARGE  FLOWRATE  24  7  m3/sec 
OUTFALL  TEMPERATURE  17  7°C 
INTAKE  TEMPERATURE  6  7°C 
AMBIENT  WATER-  TEMPERATURE  8  3°C 
DRY  BULB  TEMPERATURE  10.7  -  12.2°C 
RELATIVE  HUMIDITY  75% 

UNO  SPEEO  &  DIRECTION  0  -3  0  m/sec.  110° 
AMBIENT  CURRENT  SPEED  &  DIRECTION  0 
LAKE  SURFACE  CONDITIONS  CALM,  0  -  0.2  m  WAVES 
SKY  CONDITIONS:  CLEAR 
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Fig.  10.  Jet-regime  Study  for  1.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  18,  1972,  1115-1440  Hours.  ANL  Neg.  No.  190-893. 
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PLANT  8  METEOROLOGICAL  OAT4 


PLANT  LOAD.  495  MWe 

DISCHARGE  FLOW  RATE  24.7  m3/sec 

OUTFALL  TEMPERATURE  !7.7»C 

INTAKE  TEMPERATURE  6  7'C 

AMBIENT  WATER  TEMPERATURE  7.7  8.3'C 

ORY  3ULB  TEMPERATURE'  1C7-12.2'C 

RELATIVE  HUM'ICITY  75V 

WIND  SPEED  8  DIRECTION  0  3  0  m  sec.  110' 
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SKY  CONDITIONS  CLEAR 


0  100  200 


VELOCITY  SCALE  cm/sec 
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Fig.  11.  Jet-regime  Study  for  2.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  18.  1972,  1115-1440  Hours.  ANL  Neg.  No.  190-899. 
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PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  495  MWe 

DISCHARGE  FLOW  RATE  24  7  m3/sec 

OUTFALL  TEMPERATURE  17  7°C 

INTAKE  TEMPERATURE  6.7°C 

AMBIENT  WATER  TEMPERATURE  7.6  -  7.8°C 

DRY  3ULB  TEMPERATURE  10  7  12  ?°C 

RELATIVE  HUMIDITY  75%  ' 

WIND  SPEED  &  DIRECTION  0  •  3.0  m/sec.  110°C 
AMBIENT  CURRENT  SPEED  &  DIRECTION  0 
LAKE  SURFACE  CONDITIONS  CALM,  0  -  0.2  m  WA\£ 
SKY  CONDITIONS:  CLEAR 
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Fig.  12.  Jet-regime  Study  for  2.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  18,  1972,  1115-1440  Hours.  ANL  Neg.  No,  190-898. 
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Fig.  13.  Jet-regime  Study  for  3.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  18,  1972,  1115-1440  Hours.  ANL  Neg.  No.  190-889. 
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PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAO.  497  MWe 
DISCHARGE  FLOW  RATE  25.1  m3  sec 
OUTFALL  TEMPERATURE  21  6°C 
INTAKE  TEMPERATURE:  11  1°C 
AMBIENT  WATER  TEMPERATURE  14  3°C 
DRY  BULB  TEMPERATURE  14  5°C 
RELATIVE  HUMIDITY.  83#* 
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Fig.  14.  Jet-regime  Study  for  0.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  23,  1972,  0945-1700  Hours.  ANL  Neg.  No.  190-760  Rev.  1. 
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PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD.  497  MWe 
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Fig.  15.  Jet-regime  Study  for  l.O-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  23.  1972,  0945-1700  Hours.  ANL  Neg.  No.  190-891. 
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Fig.  16.  Jet-regime  Study  for  1.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  23.  1972,  0945-1700  Hours.  ANL  Neg.  No.  190-906. 
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DISCHARGE  FLOW  RATE:  25.1  m3  sec 
OUTFALL  TEMPERATURE  216°C 
INTAKE  TEMPERATURE  11  1°C 
AMBIENT  WATER  TEMPERATURE  11.6«»C 
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RELATIVE  HUMIDITY  830« 
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Fig.  17.  Jet-regime  Study  for  2.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  23,  1972,  0945-1700  Hours.  ANL  Neg.  No.  190-894. 
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110 

21.0 

5 

72.3 

100 

21.0 

6 

350 

110 

17.6 

7 

7.8 

140 

12.4 

8 

LAKE  BOTTOM 

9 

LAKE  BOTTOM 

10 

44.5 

90 

19.3 

11 

24.5 

60 

17.3 

12 

14.5 

20 

15.3 

13 

7.2 

350 

14.1 

14 

LAKE  BOTTOM 

15 

NOT  MEASURED 

16 

24  5 

85 

176 

17 

0 

- 

12.3 

18 

LAKE  BOTTOM 

19 

LAKE  BOTTOM 

20 

15.0 

75 

17.3 

21 

16.7 

75 

17.1 

22 

5.6 

200 

11.6 

23 

5.6 

225 

11.3 

24 

0 

- 

11.3 

25 

0 

- 

13.3 

26 

0 

- 

14  8 

27 

0 

- 

14.7 

28 

0 

- 

12.8 

29 

0 

- 

12.8 

(TEMPERATURE) 


-  TEMPERATURE  WIOTH 

■  ■  ■  ■  VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD.  497  MWe 


DISCHARGE  FLOW  RATE  25.1  m3  sec 
OUTFALL  TEMPERATURE  21  6°C 
INTAKE  TEMPERATURE  11  1°C 
AMBIENT  WATER  TEMPERATURE  11  3°C 
DRY  BULB  TEMPERATURE  14  5°C 
RELATIVE  HUMIDITY:  83c« 

WIND  SPEED  &  DIRECTION  0- 1.5  m/ sec. 
AMBIENT  CURRENT  SPEED  &  DIRECTION 
LAKE  SURFACE  CONDITIONS.  CALM 
SKY  CONDITIONS:  CLEAR 


VELOCITY  SCALE  cm  sec 


120° 

0 


0  50  100  150 


DIMENSION  SCALE  melets 


Fig.  18.  Jet-regime  Study  for  2.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
May  23,  1972,  0945-1700  Hours.  ANL  Neg.  No.  190-892. 
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STATION 

NUMBER 

CURRENT 
SPEED 
(cm  sec' 

CURRENT 

DIRECTION 

(°. 

TEMPERATURE 

CD 

1 

LAKE  BOTTOM 

2 

LAKE  BOTTOM 

3 

43.9 

80 

19.3 

4 

528 

95 

20.9 

5 

50.0 

100 

208 

6 

18.3 

120 

16  4 

7 

LAKE  BOTTOM 

8 

LAKE  BOTTOM 

9 

LAKE  BOTTOM 

10 

LAKE  BOTTOM 

11 

NOT  MEASURED 

12 

LAKE  BOTTOM 

13 

LAKE  BOTTOM 

14 

LAKE  BOTTOM 

15 

NOT  MEASURED 

16 

NOT  MEASURED 

17 

NOT  MEASURED 

18 

LAKE  BOTTOM 

19 

LAKE  BOTTOM 

20 

6J 

85 

16.3 

21 

16.7 

75 

16  1 

22 

2.8 

190 

11.2 

23 

5.6 

225 

11  6 

24 

0 

- 

10.8 

25 

0 

- 

12.1 

26 

0 

- 

14.3 

27 

0 

- 

14  3 

28 

0 

- 

12.8 

29 

0 

- 

12.6 

(TEMPERATURE) 


- TEMPERATURE  WIDTH 

-  VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAO.  497  MWe 
DISCHARGE  FLOW  RATE  25.1  m3  sec 
OUTFALL  TEMPERATURE  21i°C 
INTAKE  TEMPERATURE  ll.CC 
AMBIENT  WATER  TEMPERATURE  10  8°C 
DRY  BULB  TEMPERATURE  14.5°C 
RELATIVE  HUMIDITY:  83°« 

WIND  SPEED  &  DIRECTION  0-1.5  m/sec.  120° 
AMBIENT  CURRENT  SPEED  &  DIRECTION  0 
LAKE  SURFACE  CONDITIONS  CALM 
SKY  CONDITIONS:  CLEAR 


100 


VELOCITY  SCALE  cm/sec 


0  50  100  150 

I _ I _ I _ I 

DIMENSION  SCALE  meteis 

Fig.  19.  Jet-regime  Study  for  3.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 

May  23.  1972,  0945-1700  Hours.  ANL  Neg.  No.  190-890. 
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STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

(°) 

TEMPERATURE 

(°C) 

1 

19.5 

180 

17.0 

2 

62.3 

no 

19.0 

3 

67.8 

95 

20.2 

4 

51.1 

80 

18.5 

5 

17.8 

65 

16.7 

6 

5.6 

60 

16.5 

7 

30.0 

80 

16.7 

8 

44.5 

95 

18.8 

9 

51.1 

110 

19.2 

10 

23.4 

125 

17.0 

11 

8.9 

135 

17.0 

12 

6.7 

75 

15.3 

13 

26.1 

115 

16.3 

14 

30.0 

120 

16.5 

15 

24.4 

140 

16.0 

16 

13.9 

170 

15.5 

17 

11.1 

145 

15.7 

(TEMPERATURE) 


—  —  —  TEMPERATURE  WIDTH 
-  VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAD  500  MWe 

DISCHARGE  FLOW  RATE  24.3  -  25  1  m3  'sec 

OUTFALL  TEMPERATURE  20.3°C 

INTAKE  TEMPERATURE:  10.0°C 

AMBIENT  WATER  TEMPERATURE:  13.0  -  13.5°C 

DRY  BULB  TEMPERATURE  17.3  -  18  0°C 

RELATI  V E  HUMIDITY:  80°c 

WIND  SPEED  &  DIRECTION  0  -  2.0  m/sec;  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  6.6  cm  sec,  180°#1215hrs 

2.2  cm  sec,  245°#  1720  hrs 

LAKE  SURFACE  CONDITIONS  CALM  0  -  0.1  m  WAVES 
SKY  CONDITIONS:  CLEAR 


100 


200 


VELOCITY  SCALE  •  cmsec 


100 


150 


DIMENSION  SCALE  •  meters 


Fig.  20.  Jet-regime  Study  for  0.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
July  13,  1972,  1308-1706  Hours.  ANL  Neg.  No.  190-762  Rev.  2. 
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STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

(°C) 

1 

12.2 

115 

16  7 

2 

62.3 

110 

19  0 

3 

600 

95 

20.0 

4 

37.3 

80 

18  0 

5 

14.5 

130 

16.7 

6 

2.2 

325 

16.0 

7 

20.6 

90 

15.5 

8 

42.3 

95 

18  3 

9 

40.0 

110 

18.3 

10 

19.5 

115 

16.3 

11 

5.0 

160 

162 

12 

2.2 

115 

14.0 

13 

16.7 

120 

15.0 

14 

22.2 

120 

15.0 

15 

16.7 

145 

14.7 

16 

3.9 

195 

14.0 

17 

106 

145 

13.7 

(TEMPERATURE) 


TEMPERATURE  WIDTH 
VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAD  500  MWe 

DISCHARGE  FLOW  RATE  24  3  25  1  m3/sec 

OUTFALL  TEMPERATURE  20  3° C 

INTAKE  TEMPERATURE  10.0°C 

AMBIENT  WATER  TEMPERATURE  12  2  •  12.7°C 

DRY  BULB  TEMPERATURE  17  3  18  0°C 

RELATIVE  HUMIDITY  80% 

WIND  SPEED  &  DIRECTION  0  -  2  0  m/sec.  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  6  6  cm  sec.  180° •  1215  hrs 
2  2  cm  sec,  245°  •  1720  hrs 

LAKE  SURFACE  CONDITIONS  CALM  0  -  0  1  m  WAVES 
SKY  CONDITIONS  CLEAR 


100  200 

J_ I 


VELOCITY  SCALE  •  cm  sec 


100 

I 


150 


DIMENSION  SCALE  meters 


Fig.  21.  Jet-regime  Study  for  l.O-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
July  13,  1972,  1308-1706  Hours.  ANL  Neg.  No.  190-900  Rev.  1. 
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STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

(°) 

TEMPERATURE 

(°C) 

1 

12.2 

45 

16.0 

2 

32.2 

110 

17.0 

3 

67.8 

95 

18.5 

4 

40.0 

90 

18.5 

5 

14.5 

180 

16.0 

6 

0 

- 

15.5 

7 

16.7 

125 

15.2 

8 

36.7 

95 

17.7 

9 

37.3 

110 

18.3 

10 

11.1 

115 

16.0 

11 

3.9 

250 

15.7 

12 

6.7 

145 

12.5 

13 

111 

125 

14.5 

14 

13.3 

130 

14.0 

15 

9.5 

150 

13.5 

16 

0 

- 

12.7 

17 

6.7 

150 

12.5 

t 

(VELOCITY)  £ 

!3  ^  (TEMPERATURE) 


TEMPERATURE  WIOTH 
VELOCITY  WIOTH 


/ 


/ 


4 


PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  500  MWe 

DISCHARGE  FLOW  RATE  24  3-25  1  m3/sec 

OUTFALL  TEMPERATURE  20  3° C 

INTAKE  TEMPERATURE  10.0°C 

AMBIENT  WATER  TEMPERATURE  11  4  -  11  7°C 

DRY  BULB  TEMPERATURE  17  3  18  0°C 

RELATIVE  HUMIDITY  80% 

WIND  SPEED  &  DIRECTION  0  -  2.0  m/sec;  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  6  6  cm  sec.  180°  *  1215  hrs 
2  2  cm /sec.  245°  ©  1 720  hf s 

LAKE  SURFACE  CONDITIONS  CALM  O-Olm  WAVES 
SKY  CONDITIONS  CLEAR 


0  100  200 


VELOCITY  SCALE  -  cm  sec 


50  100  150 


DIMENSION  SCALE  meters 


Fig.  22.  Jet-regime  Study  for  1.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
July  13,  1972,  1308-1706  Hours.  ANL  Neg.  No.  190-903  Rev.  1. 
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STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

(°C) 

1 

16.7 

25 

14.5 

2 

40.0 

110 

18.5 

3 

77.8 

95 

19.5 

4 

30.0 

100 

17.5 

5 

LAKE  BOTTOM 

6 

11.1 

185 

14.3 

7 

13.9 

170 

13.5 

8 

33.4 

95 

17.7 

9 

28.9 

125 

16.5 

10 

6.7 

40 

15.0 

11 

LAKE  BOTTOM 

12 

5.6 

180 

11.5 

13 

8.3 

135 

140 

14 

11.1 

150 

13.3 

15 

11 

170 

12.3 

IS 

5.6 

320 

11.7 

17 

1.1 

180 

11.5 

- TEMPERATURE  WIDTH 

-  VELOCITY  WIDTH 


PLANT  &  METCOROLOGICAL  DATA 

PLANT  LOAD  500  MWe 
DISCHARGE  FLOW  RATE  24  3  •  25  1  m3'sec 
OUTFALL  TEMPERATURE  20  3°C 
INTAKE  TEMPERATURE  10  0°C 
AMBIENT  WATER  TEMPERATURE  11.0  -  11  6°C 
DRY  BULB  TEMPERATURE  17  3  18  0°C 
RELATIVE  HUMIDITY  80% 

WIND  SPEED  &  DIRECTION  0  2  0  m/sec.  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  6  6cm/sec.  180° •  1215 

2.2  cm/sec.  245°  *1720 

LAKE  SURFACE  CONDITIONS  CALM  0  0  1m  WAVES 
SKY  CONDITIONS:  CLEAR 


100 


200 

J 


VELOCITY  SCALE  -  cm  sec 


100 

JL 


150 

J 


DIMENSION  SCALE  •  melers 


Fig.  23.  Jet-regime  Study  for  2.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1) 
July  13.  1972,  1308-1706  Hours.  ANL  Neg.  No.  190-905  Rev.  1. 
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STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

(#C) 

1 

20.0 

15 

12.7 

2 

23.4 

70 

15.0 

3 

51.1 

95 

19.0 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

12.2 

180 

12.0 

7 

13.9 

195 

12.0 

8 

31.7 

95 

17.7 

9 

13.9 

150 

15  5 

10 

10.0 

30 

11.7 

11 

LAKE  BOTTOM 

12 

3.9 

180 

11.3 

13 

3.9 

175 

12.5 

14 

5.6 

190 

12.0 

15 

1.1 

230 

11.3 

16 

5.6 

330 

11.0 

17 

0 

- 

10.7 

^VELOCITY)  l 


(TEMPERATURE) 


- TEMPERATURE  WIDTH 

-  VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  OATA 

PLANT  LOAD  500  MWe 

DISCHARGE  FLOW  RATE  24  3  -  25  1  m3  sec 

OUTFALL  TEMPERATURE  20  3° C 

INTAKE  TEMPERATURE  10  0°C 

AMBIENT  WATER  TEMPERATURE  10  6  -  11  3°C 

DRY  BULB  TEMPERATURE  17  3  -  18  0°C 

RELATIVE  HUMIDITY  80^ 

WIND  SPEED  &  DIRECTION  0-2.0m'sec,  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  9  4  cm  sec.  180°  •  1215  hrs 

5.5  cm/sec.  180°  #  1720  hrs 

LAKE  SURFACE  CONDITIONS  CALM  0  0.1  m  WAVES 
SKY  CONDITIONS:  CLEAR 


200 


VELOCITY  SCALE  •  cm  sec 


100 


150 

j 


DIMENSION  SCALE  -  meters 


Fig.  24.  Jet-regime  Study  for  2.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
July  13,  1972,  1308-1706  Hours.  ANL  Neg.  No.  190-896  Rev.  1. 


42 


STATION 

NUMBER 

CURRENT 
SPEEO 
(cm  sect 

CURRENT 

DIRECTION 

(°> 

TEMPERATURE 

CC) 

1 

LAKE  BOTTOM 

2 

150 

35 

15.0 

3 

LAKE  BOTTOM 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

LAKE  BOTTOM 

7 

LAKE  BOTTOM 

8 

LAKE  BOTTOM 

9 

LAKE  BOTTOM 

10 

LAKE  BOTTOM 

11 

LAKE  BOTTOM 

12 

2.8 

185 

11.3 

13 

7.2 

205 

11.5 

14 

5.6 

290 

11.0 

IS 

1.1 

250 

11.0 

16 

3.9 

330 

11.0 

w 

0 

- 

10.3 

PLANT  l  METEOROLOGICAL  DATA 


PLANT  LOAD'  500  MWe 
OISCHARGE  FLO*  RATE  24  3  -  25.1  m3/sec 
OUTFALL  TEMPERATURE  20  3”C 
INTAKE  TEMPERATURE  10.0-C 
AMBIENT  *ATER  TEMPERATURE  10.3  -  11  1* 
ORY  BULB  TEMPERATURE:  17.3  •  18.0-C 


100 


150 

J 


DIMENSION  SCALE  -  meters 


Fig.  25.  Jet-regime  Study  for  3.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1) 
July  13.  1972,  1308-1706  Hours.  ANL  Neg.  No.  190-875  Rev.  1. 


43 


(VELOCITY) 


STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

(°) 

TEMPERATURE 

(°C) 

1 

8.3 

220 

17.2 

2 

25.0 

70 

19.6 

3 

55.6 

90 

23.0 

4 

55.6 

110 

23.8 

5 

41.7 

105 

21.0 

6 

11.1 

110 

19.1 

7 

9.5 

20 

18.8 

8 

26.7 

80 

20.5 

9 

27.8 

115 

20.7 

10 

29.5 

115 

20.2 

11 

23.9 

145 

20.1 

12 

8.9 

165 

20.0 

13 

13.3 

110 

19.3 

n 

12.8 

75 

19.5 

15 

12.2 

110 

19.9 

16 

6.7 

no 

19.6 

17 

0 

- 

18.9 

(TEMPERATURE) 

/ 


- —  TEMPERATURE  WIDTH 

-  VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAD:  493  MWe 
DISCHARGE  FLOWRATE:  24.7  m3/sec 
OUTFALL  TEMPERATURE:  24.5°C 
INTAKE  TEMPERATURE:  13.9°C 
AMBIENT  WATER  TEMPERATURE:  16.3  •  16.6°C 
DRY  BULB  TEMPERATURE:  14.4  -  16.7°C 
RELATIVE  HUMIDITY:  74  89% 

WIND  SPEED  &  DIRECTION  0-6.0  m/sec,  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION:  15.0  cm/sec;  180°  •  0950  hrs 

0  •  1440  hrs 

LAKE  SURFACE  CONDITIONS:  SLIGHT  CHOP;  0-0.3m  WAVES 
SKY  CONDITIONS:  CLEAR 


0  100 

I  L 


VELOCITY  SCALE  ■  cm/sec 


ISO 


DIMENSION  SCALE  •  meters 


Fig.  26.  Jet-regime  Study  for  0.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
September  9,  1972,  1045-1420  Hours.  ANL  Neg.  No.  190-763  Rev.  1. 
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/ 

/ 

/ 


/ 


/ 

/ 

/ 


/ 


STATION 

NUMBER 

CURRENT 
SPEED 
' cm  sect 

CURRENT 

DIRECTION 

TEMPERATURE 

(°C> 

1 

7.8 

200 

16.3 

2 

13.9 

90 

19.5 

3 

47.3 

85 

22.8 

4 

38.9 

90 

21.5 

5 

27.8 

no 

20.5 

6 

2.8 

90 

18.6 

7 

2.8 

30 

18.3 

8 

228 

80 

19.9 

9 

22.2 

no 

19.2 

10 

21.1 

115 

19.8 

11 

19.5 

125 

19.5 

12 

3.9 

180 

18.5 

13 

2.2 

125 

17.8 

14 

0 

- 

17.7 

15 

7.2 

125 

18.2 

16 

6.7 

125 

18.1 

17 

0 

- 

18.9 

13 


/ 


- TEMPERATURE  WIDTH 

-  VELOCITY  WIDTH 


PLANT  S  METEOROLOGICAL  DATA 


PLANT  LOAD  493  MWe  —  _ _  __ 

DISCHARGE  FLO*  RATE:  24,7  m3/sec 
OUTFALL  TEMPERATURE  24VC 
INTAKE  TEMPERATURE  I3.9“C 
AMBIENT  WATER  TEMPERATURE:  15.2  •  15.5-C 
DRY  BULB  TEMPERATURE:  14  4  •  16.7“C 
RELATIVE  HUMIDITY:  74  -89* 

WINDSPEED&  DIRECTION:  0  -  6.0  m/sec,  135- 

AMBIENT  CURRE  NT  SPEED  &  DIRECTION  15.0  cm/sec.  180“  »  0950  hrs 

0  -  1440  hrs 

LAKE  SURFACE  CONDITIONS  SLIGHT  CHOP  0  0  3  m  WAVES 
SKY  CONDITIONS  CLEAR 


0  100  200 


I  1 

VELOCITY  SCALE  cm  sec 


5C  100  150 


DIMENSION  SCALE  meters 


Fig.  27.  Jet-regime  Study  for  l.O-m  Depth  at  Point  Beach  Power  Plant  (Unit  1) 
September  9,  1972,  1045-1420  Hours.  ANL  Neg.  No.  190-904. 
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/ 


/ 


/ 


STATION 

NUMBER 

CURRENT 
SPEED 
(cm  sec) 

CURRENT 

DIRECTION 

(°) 

TEMPERATURE 

(°C) 

1 

LAKE  BOTTOM 

2 

8.3 

110 

19.6 

3 

41.7 

100 

21.3 

4 

41.7 

100 

22.2 

5 

26.1 

100 

18.1 

6 

9.5 

340 

18  4 

7  / 

0.6 

35 

17.7 

8' 

13.9 

90 

18  8 

/  9 

14.5 

105 

18.8 

'  10 

10.6 

145 

17.9 

11 

11.1 

135 

18.3 

12 

2.8 

220 

16.8 

13 

0 

— 

16.0 

14 

0 

— 

15.9 

15 

2.8 

170 

15.8 

16 

1.7 

170 

16.8 

17 

0 

- 

18.3 

/ 


(TEMPERATURE) 


TEMPERATURE  WIDTH 
VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAD  493  MWe 

DISCHARGE  FLOW  RATE  247  m3/sec 

OUTFALL  TEMPERATURE  24.5°C 

INTAKE  TEMPERATURE  13.9"C 

AMBIENT  WATER  TEMPERATURE  14  4  14  6"C 

DRY  BULB  TEMPERATURE.  14  4  16  7“C 

RELATIVE  HUMIOITY  74  89% 

WIND  SPEED  &  DIRECTION  0  6.0  m/sec,  135" 

AMBIENT  CURRENT  SPEEDS  DIRECTION  15  0  cm/sec,  180- •  0950  his 

0  #  1440  hrs 

LAKE  SURFACE  CONDITIONS:  SLIGHT  CHOP,  0-0  3  m  WAVES 
SKY  CONDITIONS:  CLEAR 


200 


VELOCITY  SCALE  cm  sec 


0 

L 


5C 


10C 


DIMENSION  SCALE  meters 


150 


Fig.  28. 


Jet-regime  Study  for  1.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
September  9,  1972,  1045-1420  Hours.  ANL  Neg.  No.  190-897. 
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STATION 

NUMBER 

CURRENT 
SPEED 
icm  seci 

CURRENT 

DIRECTION 

n 

TEMPERATURE 

(°C) 

1 

LAKE  BOTTOM 

? 

LAKE  BOTTOM 

3 

44.5 

100 

21.8 

4 

44.5 

80 

21.6 

5 

15.0 

110 

17.7 

6 

15.0 

340 

16.7 

7 

6.7 

165 

16.6 

8 

1.7 

100 

179 

9 

0.6 

115 

17.6 

10 

12.2 

135 

17.7 

11 

0.6 

135 

16.7 

12 

5.6 

325 

16.1 

13 

0 

- 

15.5 

14 

0 

- 

14.8 

15 

1.1' 

180 

14.8 

16 

/  0 

- 

15.7 

17  ✓ 

2.2 

180 

16.1 

»  / 


/ 

/ 

/ 


/ 


(TEMPERATURE) 


TEMPERATURE  WIDTH 
VELOCITY  WIDTH 


PLANT  &  METEOROLOGICAL  DATA 

PLANT  LOAD:  <93  MWe 

DISCHARGE  FLOW  RATE  24  7  m3/sec 

OUTFALL  TEMPERATURE  24.5-C 

INTAKE  TEMPERATURE  13  9“C 

AMBIENT  WATER  TEMPERATURE  M  l  •  14  4°C 

DRY  BULB  TEMPERATURE:  14  4  •  16.7-C 

RELATIVE  HUMIDITY  74  89% 

WIND  SPEED  S  DIRECTION:  0  •  6.0  m/sec;  135" 

AMBIENT  CURRENT  SPEEDS  DIRECTION:  15  0  cm/sec.  180”» 0950 hrs 

0  #  1440  Ins 

LAKE  SURFACE  CONDITIONS:  SLIGHT  CHOP,  0  0.3  m  WAVES 
SKY  CONDITIONS:  CLEAR 


100  200 

J_ I 


VELOCITY  SCALE  cm  sec 


0 

L 


50  100  150 


DIMENSION  SCALE  meters 


Fig.  29.  Jet-regime  Study  for  2.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
September  9,  1972,  1045-1420  Hours.  ANL  Neg.  No.  190-902. 
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STATION 

NUMBER 

CURRENT 
SPEED 
(cm  seci 

CURRENT 

DIRECTION 

(°> 

TEMPERATURE 

(°C) 

1 

LAKE  BOTTOM 

2 

LAKE  BOTTOM 

3 

LAKE  BOTTOM 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

L/KE  BOTTOM 

7 

7.2 

165 

15.6 

8 

6.1 

165 

15.6 

9 

0 

- 

16.2 

10 

4  4 

180 

16.2 

11 

5.6 

•290 

15.1 

12 

6.7 

345 

14.6 

13 

0 

- 

149 

14 

0 

- 

14.6 

15 

1.7 

215 

14  6 

16 

4  4 

235 

14.6 

17 

3.9 

200 

15.1 

/ 


/ 


I 

l 

\ 


/ 


9 


/ 


- TEMPERATURE  WIDTH 


PLANT  LOAD  493  MWe 

DISCHARGE  FLOW  RATE  24  7  m3  sec 

OUTFALL  TEMPERATURE  24  5°C 

INTAKE  TEMPERATURE  13  9*C 

AMBIENT  WATER  TEMPERATURE  13  9  14  3°C 

DRY  BULB  TEMPERATURE  14  4  •  16  7«C 

RELATIVE  HUMIDITY  74  89% 

WIND  SPEED  &  DIRECTION  0- 60m  sec,  135° 
AMBIENT  CURR ENT  SPEED  &  DIRECTION  11  1  cm/sec, 

2.8  cm/sec; 


180° #  0950  hrs 
1 70°  •  1440  hrs 


LAKE  SURFACE  CONDITIONS 
SKY  CONDITIONS  CLEAR 


SLIGHT  CHOP,  0  -0  3M  WAVES 


100 


200 


VELOCITY  SCALE  cm,  sec 


0.  50  100  150 


DIMENSION  SCALE  melets 


Fig,  30.  Jet-regime  Study  for  2.5-m  Depth  at  Point  Beach  Power  Plant  (Unit  1): 
September  9,  1972,  1045-1420  Hours.  ANL  Neg.  No.  190-885. 
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STATION 

NUMBER 

CURRENT 
SPEED 
(cm/ sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

(°C) 

1 

LAKE  BOTTOM 

2 

LAKE  BOTTOM 

3 

LAKE  BOTTOM 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

LAKE  BOTTOM 

7 

7.2 

180 

14.9 

8 

LAKE  BOTTOM 

9 

LAKE  BOTTOM 

10 

LAKE  BOTTOM 

11 

LAKE  BOTTOM 

12 

LAKE  BOTTOM 

13 

0 

— 

14.6 

14 

0 

- 

14.5 

15 

1.7 

245 

14.5 

16 

2.8 

235 

14.5 

17 

0 

- 

14.5 

1*5 


7 


9 


10 


PLANT  S  METEOROLOGICAL  DATA 


PLANT  LOAD  493  MWe 

DISCHARGE  FLO*  RATE:  24  7  m3/sec 

OUTFALL  TEMPERATURE:  24.5»C 

INTAKE  TEMPERATURE  13  9°C 

AMBIENT  WATER  TEMPERATURE  13  8  14  3»C 

DRY  BULB  TEMPERATURE  14  4-16.7°C 


0  SO  100  150 


DIMENSION  SCALE  -  meters 


Fig.  31.  Jet-regime  Study  for  3.0-m  Depth  at  Point  Beach  Power  Plant  (Unit  1) 
September  9.  1972,  1045-1420  Hours.  ANL  Neg.  No.  190-888. 
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STATION 

NUMBER 

CURRENT 
SPEED 
(cm  sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

(°C) 

1 

33.9 

285 

22.3 

2 

25.6 

300 

22.9 

3 

31.1 

320 

22.9 

4 

16.7 

335 

22.4 

5 

24.5 

345 

23.0 

6 

22.2 

350 

20.4 

7 

21.1 

320 

20.9 

8 

24.5 

320 

22.2 

9 

23.4 

340 

21.7 

10 

19.5 

355 

21.8 

11 

19.5 

10 

22.0 

12 

18.3 

355 

20.2 

13 

15.0 

10 

19.9 

14 

18.3 

10 

16.8 

15 

20.0 

360 

20.9 

16 

47.3 

320 

23.5 

17 

69.5 

320 

23  5 

18 

61.2 

320 

23.6 

19 

55.6 

330 

23.6 

PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  405  MWe 
DISCHARGE  FLOW  RATE  25  6  m3  sec 
OUTFALL  TEMPERATURE  23  6-C 
INTAKE  TEMPERATURE  13  9-C 
AMBIENT  WATER  TEMPERATURE  16.8°C 
DRY  BULB  TEMPERATURE  20.5-C 
RELATIVE  HUMIDITY  80% 

WIND  SPEED  &  DIRECTION  6  3  msec,  220- 
AMBIENT  CURRENT  SPEED  &  DIRECTION  26.7  cm/sec,  020- 
LAKE  SURFACE  CONDITIONS  CHOPPY  0  5  10  m  WAVES 
SKY  CONDITIONS;  CLOUDY 


12  \ 


13 


1 


0  100  200 


VELOCITY  SCALE  •  cm/sec 

0  SO  100 


DIMENSION  SCALE  meteis 


Fig.  32.  Jet-regime  Study  for  0.5-m  Depth  at  Palisades  Power  Plant: 

June  14, 1972, 1000-1348  Hours.  ANL  Neg.  No.  190-757  Rev.  1. 
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STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

CC) 

1 

21.1 

310 

22.4 

2 

17.8 

310 

22.7 

3 

16.7 

310 

22.7 

4 

13.3 

20 

21.6 

5 

17.8 

345 

226 

6 

22.2 

360 

23.3 

7 

21.7 

330 

19.7 

8 

21.1 

325 

21.3 

9 

16.7 

340 

20  6 

10 

13.3 

15 

19.8 

11 

13.3 

70 

16.1 

12 

18.3 

360 

176 

13 

16.1 

10 

16.8 

14 

189 

15 

165 

IS 

10.0 

360 

17.8 

16 

44.5 

340 

236 

17 

61.2 

325 

23  5 

18 

50.0 

300 

236 

19 

55.6 

330 

236 

PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  «5  Mi* 

DISCHARGE  FLO*  RATE  25  6  aVuc 
OUTFALL  TEMPERATURE  23VC 
INTAKE  TEMPERATURE  13V  C 
AMBIENT  (ATEN  TEMPERATURE  II  l*C 
DRY  BULB  TEMPERATURE  2t.S*C 
RELATIVE  HUMIDITY  IM 
UNO  SPEEO  I  DIRECTION  6.3  b/mc;  220* 

AMBIENT  CURRENT  SPEED  l  DIRECTION  26  7  ca/scc.  020* 
LAKE  SURFACE  CONDITIONS  CHOPPY.  0.5  •  1.0  m  NAVES 
SKY  CONDITIONS  CLOUOY 


12  \ 


0  100  200 

I _ I _ I 

VELOCITY  SCALE  -  cm/sec 

0  50  100  150 

I _ I _ I _ I 

DIMENSION  SCALE  ■  meters 


Fig.  33.  Jet-regime  Study  for  1.0-m  Depth  at  Palisades  Power  Plant: 
June  14,  1972,  1000-1348  Hours.  ANL  Neg.  No.  190-872. 
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STATION 

NUMBER 

CURRENT 
SPEED 
(cm  sec) 

CURRENT 

DIRECTION 

(#) 

TEMPERATURE 

(°C) 

1 

18.3 

315 

21.7 

2 

16.7 

325 

21J 

3 

14.5 

35 

22.2 

4 

15.6 

55 

21.3 

5 

11.1 

360 

22.6 

6 

20.6 

360 

23.3 

7 

22.2 

355 

18.3 

8 

22.2 

355 

19.1 

9 

11.7 

345 

17.3 

10 

9.5 

45 

16.0 

11 

LAKE  BOTTOM 

12 

16.7 

15 

16.2 

13 

14.5 

360 

16.5 

14 

18  9 

20 

16.3 

15 

1.1 

20 

16  4 

16 

25.0 

320 

23.5 

17 

55.6 

325 

23.5 

18 

55.6 

300 

23.6 

19 

44.5 

325 

23.6 

PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  405  MWe 
DISCHARGE  FLOW  RATE  25  6  m3sec 
OUTFALL  TEMPERATURE  23  6°C 
INTAKE  TEMPERATURE  13  9°C 
AMBIENT  WATER  TEMPERATURE  16.0°C 
DRY  BULB  TEMPERATURE  20  5° C 
RELATIVE  HUMIDITY  80#. 

WIND  SPEED  &  DIRECTION  6  3  m  sec.  220° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  26.7  cm/sec.  020° 
LAKE-SURFACE  CONDITIONS  CHOPPY.  05  10  m  WAVES 
SKY  CONDITIONS  CLOUDY 


0  100  200 

I _ I _ I 

VELOCITY  SCALE  •  cm/sec 

0  50  100  150 

I _ I _ I _ I 

DIMENSION  SCALE  meins 


Fig.  34.  Jet-regime  Study  for  1.5-m  Depth  at  Palisades  Power  Plant 
June  14,  1972.  1000-1348  Hours.  ANL  Neg.  No.  190-879. 
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STATION 

NUMBER 

CURRENT 
SPEEO 
(cm  sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

CC) 

1 

21  1 

10 

21.5 

2 

13.9 

10 

21.5 

3 

139 

90 

21.7 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

LAKE  BOTTOM 

7 

20.0 

10 

17.3 

8 

22.2 

20 

17.3 

9 

1.1 

350 

16.4 

10 

LAKE  BOTTOM 

11 

LAKE  BOTTOM 

12 

139 

20 

15.8 

13 

156 

360 

16.4 

14 

178 

15 

16  2 

15 

LAKE  BOTTOM 

16 

LAKE  BOTTOM 

17 

LAKE  BOTTOM 

18 

LAKE  BOTTOM 

19 

47  3 

330  |  236 

PUNT  4  METEOROLOGICAL  OATA 


PUNT  LOAD  405  MWe 
OISCHARGE  FLO*  RATE  25  6  m3  sec 
OUTFALL  TEMPERATURE  23  6"  C 
INTAKE  TEMPERATURE  13  9”C 
AMBIENT  *ATER  TEMPERATURE  15.8°C 
DRY  BULB  TEMPERATURE..  20,5°C 
REUTIVE  HUMIDITY.  80S 
WIND  SPEEO  &  DIRECTION  6  3  m  sec,  220“ 

AMBIENT  CURRENT  SPEED  &  DIRECTION  26  7  cm/sec.  020° 
LAKE  SURFACE  CONDITIONS  CHOPPY.  0  5  10  m  WAVES 
SKY  CONDITIONS  CLOUDY 


12  1 

0  100  200 

I _ I _ I 

VELOCITY  SCALE  -  cm/sec 

0  50  100  150 

I _ I _ I - i 

DIMENSION  SCALE  meters 


Fig.  35.  Jet-regime  Study  for  2.0-m  Depth  at  Palisades  Power  Plant: 

June  14, 1972, 1000-1348  Hours.  ANL  Neg.  No.  190-877  Rev.  1. 
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STATION 

NUMBER 

CURRENT 
SPEED 
(cm  sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

(°C) 

1 

LAKE  BOTTOM 

2 

LAKE  BOTTOM 

3 

LAKE  BOTTOM 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

LAKE  BOTTOM 

7 

LAKE  BOTTOM 

8 

LAKE  BOTTOM 

9 

LAKE  BOTTOM 

10 

LAKE  BOTTOM 

11 

LAKE  BOTTOM 

12 

11.7 

35 

15.6 

13 

111 

15 

162 

14 

13.9 

15 

16  2 

15 

LAKE  BOTTOM 

16 

LAKE  BOTTOM 

17 

LAKE  BOTTOM 

18 

LAKE  BOTTOM 

19 

LAKE  BOTTOM 

PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  405  MWe 
DISCHARGE  FLOW  RATE  25  6  m3  sec 
OUTFALL  TEMPERATURE  23  6°C 
INTAKE  TEMPERATURE  13  9°C 
AMBIENT  WATER  TEMPERATURE  15.6°C 
ORY  BULB  TEMPERATURE  20  5° C 
RELATIVE  HUMIDITY  80** 

WINO  SPEED  &  DIRECTION  6  3  m  sec,  220° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  24.5  cm/sec.  020° 
LAKE  SURFACE  CONDITIONS  CHOPPY.  0  5  10  m  WAVES 
SKY  CONDITIONS  CLOUDY 


13 


1 


0  100  200 

I _ I _ I 

VELOCITY  SCALE  cm/sec 

0  SO  100  ISO 


DIMENSION  SCALE  meleis 


Fig.  36.  Jet-regime  Study  for  2.5-m  Depth  at  Palisades  Power  Plant: 
June  14,  1972,  1000-1348  Hours.  ANL  Neg.  No.  190-887. 
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STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

C) 

TEMPERATURE 

(#C) 

1 

1.7 

230 

26  5 

2 

21.1 

285 

27.0 

3 

34.5 

280 

28  5 

4 

19.5 

295 

28.7 

5 

16.7 

325 

29.0 

6 

3.3 

18 

28.2 

7 

8.3 

280 

26.0 

8 

19.5 

275 

27,0 

9 

16.7 

300 

27.0 

10 

14.5 

330 

28.0 

11 

11.7 

330 

27.5 

12 

6.7 

255 

25.0 

13 

156 

280 

26.8 

14 

15.6 

290 

25.5 

15 

11.7 

300 

25.5 

16 

12.2 

300 

25.8 

17 

10.6 

305 

26  6 

18 

545 

300 

286 

19 

72.3 

315 

28.7 

20 

61.2 

300 

28.7 

7  •- 


10  v 


11  v 


PLANT  t  METEOROLOGICAL  OATA 


PLANT  LOAD  420  MWe 

DISCHARGE  FLOW  RATE.  25.6  m3/sec 

OUTFALL  TEMPERATURE  28  7* C 

INTAKE  TEMPERATURE  20  0”C 

AMBIENT  WATER  TEMPERATURE  23.5  -  2«.5*C 

DRY  BULB  TEMPERATURE  22.7  -  25.6”  C 

RELATIVE  HUMIDITY:  86% 

WIND  SPEED  &  DIRECTION:  0  •  2.0  m/sec;  135“ 

AMBIENT  CURRENT  SPEED  &  DIRECTION  8.6  cm/sec,  180” 
LAKE  SURFACE  CONDITIONS:  CALM 
SKY  CONDITIONS:  CLOUOY 


12  «■ 


13  — 


14  — 


15  - 


16  — 


17- 


0  100  200 


VELOCITY  SCALE  -  cm/sec 
0  50  100 


DIMENSION  SCALE  •  melers 


Fig.  37.  Jet-regime  Study  for  0.5-m  Depth  at  Palisades  Power  Plant: 

July  19.1972,0922-1414  Hours.  ANL  Neg.  No.  190-758  Rev.  1. 
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PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  420  MWe 

DISCHARGE  FLOW  RATE  25  6  m3/sec 

OUTFALL  TEMPERATURE  28.7°C 

INTAKE  TEMPERATURE  20  0*C 

AMBIENT  WATER  TEMPERATURE  22.1  •  23.7-C 

DRY  BULB  TEMPERATURE  22  7  25  6°C 

RELATIVE  HUMIDITY  86% 

WIND  SPEED  &  DIRECTION  0  2  0  m/sec.  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  8.6  cm/sec.  180- 
LAKE  SURFACE  CONDITIONS  CALM 
SKY  CONDITIONS:  CLODUY 


12  • 


13  • 


16* 


0  100  200 


VELOCITY  SCALE  •  cm/sec 
0  50  100 


DIMENSION  SCALE  •  meters 


Fig.  38.  Jet-regime  Study  for  l.O-m  Depth  at  Palisades  Power  Plant: 
July  19,  1972,  0922-1414  Hours.  ANL  Neg.  No.  190-873. 
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STATION 

NUMBER 

CURRENT 
SPEED 
(cm  sec) 

CURRENT 

DIRECTION 

(°) 

TEMPERATURE 

CC) 

1 

LAKE  BOTTOM 

2 

.AKE  BOTTOM 

3 

6  7 

300 

27.5 

4 

15.6 

295 

28.4 

5 

4  4 

320 

27.6 

6 

LAKE  BOTTOM 

7 

3.3 

200 

22.5 

8 

56 

190 

222 

9 

1.1 

235 

22.3 

10 

0.6 

280 

22.0 

11 

0 

- 

22.0 

12 

0 

- 

22.0 

13 

3.3 

250 

21.5 

14 

2.2 

270 

21.6 

15 

1.7 

270 

224 

16 

11 

240 

22.0 

17 

2.8 

270 

23.1 

18 

345 

300 

28.7 

19 

57.0 

345 

28  7 

20 

48  8 

300 

28  7 

8  t 


9» 


PUNT  &  METEOROLOGICAL  DATA 


PUNT  LOAD  420  Ml e 

DISCHARGE  FLOW  RATE  25  6  m3,sec 

OUTFALL  TEMPERATURE  28.7*C 

INTAKE  TEMPERATURE  20  0°C 

AMBIENT  WATER  TEMPERATURE  21.3  •  22.5*C 

DRY  BULB  TEMPERATURE  22  7  -  25  6*C 

REUTIVE  HUMIDITY  86% 

WIND  SPEED  &  DIRECTION  0  2  0  m  sec,  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  8  6  cm /sec.  180° 
UKE  SURFACE  CONDITIONS  CALM 
SKY  CONDITIONS:  CLOUDY 


0 

L 


13. 


0  100  200 

1 _ I _ I 

VELOCITY  SCALE  ■  cm/sec 

50  100 


DIMENSION  SCALE  meteis 


ISO 


Fig.  39.  Jet-regime  Study  for  1.5-m  Depth  at  Palisades  Power  Plant 
July  19.  1972,  0922-1414  Hours.  ANL  Neg.  No.  190-884. 


57 


N 


il 


STATION 

NUMBER 

CURRENT 

SPEED 

(cm/sec) 

CURRENT 

DIRECTION 

(•) 

TEMPERATURE 

(°C) 

1 

LAKE  BOTTOM 

2 

LAKE  BOTTOM 

3 

8.3 

30  J_ 26.0 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

LAKE  BOTTOM 

7 

4  4 

180 

220 

8 

5.6 

180 

22.3 

9 

1.7 

235 

22.3 

10 

4  4 

200 

22.0 

11 

0.6 

250 

21.6 

12 

0 

21.5 

13 

2.2 

255 

21.4 

14 

1.1 

220 

21.5 

15 

1.7 

200 

21.4 

16 

3.9 

215 

21.8 

17 

2.8 

205 

22.0 

18 

LAKE  BOTTOM 

19 

48  9 

330 

28 .3 

20 

51.5 

300 

28  7 

10  . 


PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAO  420  MWe 
DISCHARGE  FLOW  RATE  25  6  m3/sec 
OUTFALL  TEMPERATURE  28.7°C 
INTAKE  TEMPERATURE  20.0°C 
AMBIENT  WATER  TEMPERATURE  21,4°C 
DRY  BULB  TEMPERATURE:  22  7  25  6°C 
RELATIVE  HUMIDITY  86% 

WIND  SPEED  &  DIRECTION  0  2.0  m/sec ,  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  8.6  cm/sec,  180° 
LAKE  SURFACE  CONDITIONS  CALM 
SKY  CONDITIONS  CLOUDY 


0  100  200 

I  I J 

VELOCITY  SCALE  -  cm/sec 

0  50  100  150 

I _ I _ I _ I 

DIMENSION  SCALE  •  meters 


Fig.  40.  Jet-regime  Study  for  2.0-m  Depth  at  Palisades  Power  Plant: 

July  19.1972.0922-1414  Hours.  ANLNeg.No.  190-871  Rev.  1. 
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STATION 

NUMBER 

CURRENT 
SPEED 
i  cm  sed 

CURRENT 

DIRECTION 

TEMPERATURE 

(°C> 

1 

LAKE  BOTTOM 

2 

LAKE  BOTTOM 

3 

LAKE  BOTTOM 

4 

LAKE  BOTTOM 

5 

LAKE  BOTTOM 

6 

LAKE  BOTTOM 

7 

LAKE  BOTTOM 

8 

22 

140 

21  8 

9 

3  3 

170 

22  1 

10 

5  6 

185 

22  0 

11 

0  6 

220 

21.5 

12 

0  6 

205 

21  0 

13 

11 

255 

21  0 

14 

1  1 

185 

21  4 

15 

1  1 

190 

21  0 

16 

2.8 

210 

21  8 

17 

3.3 

205 

21.5 

18 

LAKE  BOTTOM 

19 

LAKE  BOTTOM 

2C 

LAKE  BOTTOM 

II  • 


PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  420  MWe 
DISCHARGE  FLOW  RATE  25  6  m3  sec 
OUTFALL  TEMPERATURE  28  7° C 
INTAKE  TEMPERATURE  20  0°C 
AMBIENT  WATER  TEMPERATURE  21  0°C 
DRY  BULB  TEMPERATURE  22  7  25  6° C 
RELATIVE  HUMIDITY  863« 

WINO  SPEEO  &  DIRECTION  0  20  in  sec.  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  5  0  cm' sec.  200° 
LAKE  SURFACE  CONDITIONS  CALM 
SKY  CONDITIONS  CLOUDY 
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Fig.  41.  Jet-regime  Study  for  2.5-m  Depth  at  Palisades  Power  Plant: 
July  19,  1972,  0922-1414  Hours.  ANL  Neg.  No.  190-881. 
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PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD  420  MWe 
DISCHARGE  FLOW  RATE  25  6  m3  sec 
OUTFALL  TEMPERATURE  28.7“C 
INTAKE  TEMPERATURE  20  0-C 
AMBIENT  WATER  TEMPERATURE  20.6°C 
DRY  BULB  TEMPERATURE  22  7  25  6°C 
RELATIVE  HUMIDITY  86% 

WIND  SPEED  &  DIRECTION  0  2.0  m/sec,  135° 

AMBIENT  CURRENT  SPEED  &  DIRECTION  5.0  cm/sec,  200“ 
LAKE  SURFACE  CONDITIONS  CALM 
SKY  CONDITIONS:  CLOUDY 
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Fig.  42.  Jet-regime  Study  for  3.0-m  Depth  at  Palisades  Power  Plant 
July  19,  1972,  0922-1414  Hours.  ANL  Neg.  No.  190-882  Rev.  1, 
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PLANT  &  METEOROLOGICAL  DATA 
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PUNT  LOAD:  434  We 

DISCHARGE  FLOWRATE:  25.6  m3/»ec 

OUTFALL  TEMPERATURE:  22J*C 

INTAKE  TEMPERATURE:  12.2*C 

AMBIENT  WATER  TEMPERATURE:  13.0  ■  13.3*C 

DRY  BULB  TEMPERATURE:  6  5  •  11.1'C 


REUTIVE  HUMIDITY  67% 

WIND  SPEEO  i  DIRECTION:  0  -  4.5  a/sec ;  150  •  170* 
AMBIENT  CURRENT  SPEED  A  DIRECTION:  0.9  ca/sec;  345* 
UKE  SURFACE  CONDITIONS:  CALM,  0  -  0.3  a  WAVES 
SKY  CONDITIONS:  PARTLY  CLOUDY 
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43.  Jet-regime  Study  for  0.5-m  Depth  at  Palisades  Power  Plant: 
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Fig.  44.  Jet-regime  Study  for  l.O-m  Depth  at  Palisades  Power  Plant: 
October  10,  1972,  1025-1550  Hours.  ANL  Neg.  No.  190-883. 
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PUNT  &  METEOROLOGICAL  DATA 


PUNT  LOAD  434  MWe 

DISCHARGE  FLOW  RATE:  25.6  m3  sec 

OUTFALL  TEMPERATURE  22  Z»C 

INTAKE  TEMPERATURE  12.2”C 

AMBIENT  WATER  TEMPERATURE.  12  7  •  13.0“C 

DRY  BULB  TEMPERATURE:  8.9  -  1 1.1-C 

REUTIVE  HUMIDITY:  675. 

WIND  SPEED  &  DIRECTION.  0  •  4.5  m  sec,  150  •  170” 
AMBIENT  CURRENT  SPEED  &  DIRECTION  8.9  cm/sec;  345“ 
UKE  SURFACE  CONDITIONS:  CALM,  0-0.3  m  WAVES 
SKY  CONDITIONS  PARTLY  CLOUDY 
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Fig.  45.  Jet-regime  Study  for  1.5-m  Depth  at  Palisades  Power  Plant 
October  10.  1972,  1025-1550  Hours.  ANL  Neg.  No.  190-901. 
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PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD:  434  MWe 

DISCHARGE  FLOW  RATE  25.6  m3'sec 

OUTFALL  TEMPERATURE  22  2°C 

INTAKE  TEMPERATURE:  12.2°C 

AMBIENT  WATER  TEMPERATURE  12.8  -  13.0°C 

DRY  BULB  TEMPERATURE:  8.9  -  1U°C 

RELATIVE  HUMIDITY  67% 

WIND  SPEED  &  DIRECTION:  0-4.5  m/sec;  150  -  170° 
AMBIENT  CURRENT  SPEEO&  DIRECTION  8.9  cm/sec.  345° 
LAKE  SURFACE  CONDITIONS:  CALM.  0-0.3  m  WAVES 
SKY  CONDITIONS:  PARTLY  CLOUDY 
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Fig.  46.  Jet-regime  Study  for  2.0-m  Depth  at  Palisades  Power  Plant 
October  10.  1972,  1025-1550  Hours.  ANL  Neg.  No.  190-876. 
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PLANT  4  METEOROLOGICAL  DATA 


PUNT  LOAD  434  Mfe 

DISCHARGE  FLON  RATE  25.6  m3  set 

OUTFALL  TEMPERATURE  22  2”C 

INTAKE  TEMPERATURE  12.2'C 

AMBIENT  NATER  TEMPERATURE  12  7  13  0*C 

DRY  BULB  TEMPERATURE  6.9  -  ll.l-C 

REUTIVE  HUMIDITY  67\ 

»IND  SPEED  4  DIRECTION  0-4.5  m  sec,  150  -  170” 
AMBIENT  CURRENT  SPEED  4  DIRECTION  2  6  cm/sec.  250" 
UKE  SURFACE  CONDITIONS:  CALM,  0-0.3  m  NAVES 
SKY  CONDITIONS  PARTLY  CLOUOY 
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Fig.  47. 


Jet-regime  Study  for  2.5-m  Depth  at  Palisades  Power  Plant: 
October  10,  1972,  1025-1550  Hours.  ANL  Neg.  No.  190-886. 
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PLANT  &  METEOROLOGICAL  DATA 


PLANT  LOAD:  434  MWe 

DISCHARGE  FLOW  RATE:  25.6  m3/sec 

OUTFALL  TEMPERATURE  22  2°C 

INTAKE  TEMPERATURE:  12.2°C 

AMBIENT  WATER  TEMPERATURE:  12.7  -  13  0°C 

DRY  BULB  TEMPERATURE:  8.9  -  11.1°C 

RELATIVE  HUMIDITY  67% 

WIND  SPEED  &  DIRECTION  0-4.5  m/sec,  150  -  170° 
AMBIENT  CURRENT  SPEED  &  DIRECTION  2.8  cm/sec,  250° 
LAKE  SURFACE  CONDITIONS:  CALM,  0-0.3  m  WAVES 
SKY  CONDITIONS  PARfLY  CLOUDY 
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Fig.  48. 


Jet-regime  Study  for  3.0-m  Depth  at  Palisades  Power  Plant: 
October  10,  1972,  1025-1550  Hours.  ANL  Neg.  No.  190-874. 


V.  EXTRACTION  OF  JET  CHARACTERISTICS 
FROM  FIELD  MEASUREMENTS 


In  order  to  examine  the  features  of  the  temperature  and  velocity  dis¬ 
tributions  in  the  near -field  region  of  the  discharges,  and  additionally  to  facil¬ 
itate  comparisons  with  analytical -model  predictions,  it  was  decided  to  assume 
a  functional  form  for  the  distributions  and  fit  these  distributions  to  the  field 
measurements.  The  choice  of  the  functional  forms  was  based  on  preliminary 
examination  of  the  Point  Beach  data.  The  functions  chosen  use  Gaussian  dis¬ 
tributions  for  the  lateral  excess -temperature  and  excess -velocity  profiles. 

The  excess  temperature  0  is  defined  as  the  difference  between  the  measured 
temperature  T  and  the  ambient  temperature  T^  of  the  lake  water  at  that 
depth.  The  excess  velocity  u^  is  defined  as  the  component  of  the  measured 
velocity  parallel  to  the  centerline  of  the  jet,  minus  the  component  of  the 
ambient  current  parallel  to  the  centerline.  The  ambient  current  u^  was 
taken  to  be  parallel  to  shore  for  this  purpose.  The  widths  of  these  profiles 
are  assumed  to  grow  linearly  with  distance  from  the  outfall,  s,  as  measured 
along  the  jet  centerline.  The  centerline  temperature  excess  9C  and  the  cen¬ 
terline  velocity  excess  uc£  are  assumed  to  decay  as  the  inverse  of  the  square 
root  of  the  distance  from  the  outfall.  A  quadratic  form  was  chosen  for  the 
trajectories  of  the  centerline. 

The  temperature  excess  at  a  lateral  distance  T]  from  the  temperature 
centerline  at  a  distance  s  from  the  outfall  is  thus  given  by 

/I  \[VWT(s)l2 

e(Ti,s)  =  ec(s)  exp  {-(in  z)[Ti/wT(s)]2 }  = 

whe  re 


ab0 

90A  for  s  £  — ~ 

A 


and 


WT(s)  =  Cb0  +  VS. 

In  the  above  expression,  90  is  the  initial  temperature  excess  at  the 
outfall,  b0  is  the  width  of  the  outfall,  and  W ^  is  the  half -width  of  the  excess 
temperature  distribution.  The  half -width  is  defined  as  the  lateral  distance 
at  which  the  excess  temperature  drops  off  to  one -half  the  centerline  value. 
The  dimensionless  parameters  A,  a,  C,  and  y  are  free  parameters  to  be 
determined  by  fitting  to  the  field  data. 
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To  specify  the  centerline,  a  coordinate  system  is  setup  with  the  origin 
at  the  center  of  the  outfall.  The  positive  y  axis  is  perpendicular  to  the  average 
shoreline  in  the  off-shore  direction,  the  positive  z  axis  is  vertically  upward, 
and  the  positive  x  axis  is  chosen  to  be  orthogonal  to  the  other  two  so  as  to 
form  a  right-handed  system  in  the  conventional  sense  (x,y,z).  The  tempera¬ 
ture  centerline  trajectory  is  most  easily  specified  in  parametric  form: 

kt 

xcT  =  ?  cos  (RT^o)  -  sin  (rt^o) 

and 

kt 

Y c t  =  ?  sin  (RtPo)  +  ^  cos  (RTPo) 


for 


§  2:  0. 


The  trajectory  represents  a  parabola  passing  through  the  origin  at  an 
angle  of  R^p^o  with  respect  to  the  positive  x  axis.  The  parametric  variable  § 
is  introduced  only  to  simplify  the  form  of  these  equations.  The  angle  of  the 
outfall  with  respect  to  the  positive  x  axis  is  f30>  and  R-p  and  K-j  are  two  dimen¬ 
sionless  free  parameters  to  be  determined  by  fitting  to  the  field  data. 


The  case  for  which  R-p  =  1.0  and  Kp,  =  0.0  corresponds  to  a  jet 
centerline  that  is  directed  straight  out  from  the  outfall.  Positive  values 
of  Krp  correspond  to  the  jet  bending  to  the  left,  and  negative  values  corre¬ 
spond  to  bending  to  the  right.  The  quantity  Rrpj50  determines  the  initial  angle 
of  the  centerline  at  the  origin  (discharge). 


The  velocity  excess  (assumed  parallel  to  the  centerline)  at  a  lateral 
distance  T|  from  the  velocity  centerline  is  given  by  the  following  expression: 

ue(T|,s)  =  uCE(s)(y 

where 


Pb0 

u0B  for  s  £  — — 

B2 


and 


Wu(s)  =  Db0  +  6s. 
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These  velocity -exces s  equations  have  the  same  form  as  the  expressions  for 
the  temperature  excess,  except  that  u0  is  the  outfall  velocity  averaged  over 
the  entire  outfall.  The  velocity  centerline  has  a  form  analogous  to  that  for 
the  temperature  centerline: 

K 

xcu  =  5  cos  (RuM  -  §2  sin  (Ru3o) 


and 


Ku 

ycu  =  5  sin  (RuE>o)  +  52  cos  (RUP„). 

The  parameters  B,  j3,  D,  6,  Ru,  and  Ku  are  determined  by  fitting  to  the  data. 

The  temperature  and  velocity  functions  described  above,  each  with 
only  six  free  parameters,  are  fairly  restrictive.  However,  due  to  the  limited 
amount  of  data  available,  it  was  felt  that  more  general  functions  with  more 
free  parameters  would  be  unwarranted.  As  a  consequence,  only  the  general 
trends  of  the  resulting  fitted  function  can  be  considered  significant.  The 
details  are  artifacts  of  the  functions  chosen. 

The  fitting  procedure  was  carried  out  on  the  computer  by  a  minimiza¬ 
tion  technique.  The  root -mean -square  deviations  of  the  functions  from  the 
data  are  defined  as  follows: 


(aT;A,o,C,  y,RT,KT)  = 


N 

I 

i=i 


(0  Di "  0Fi) 

N  -  6 


-,1  /2 


and 


(ou;B,  j3,D,  6,RU,KU) 


where  N  is  the  number  of  data  points,  is  the  measured  temperature 

excess  at  the  ith  data  point,  is  the  temperature  excess  calculated  from 

the  function,  u^.  is  the  component  of  the  measured  velocity  excess  parallel 
to  the  velocity  centerline,  and  u^'  is  the  velocity  excess  calculated  from 
the  function.  The  6  in  the  denominator  is  introduced  to  account  for  the 
six  degrees  of  freedom  associated  with  the  six  free  parameters  of  the  fitting 
function.  The  final  values  of  the  parameters  were  then  chosen  to  be  those 
values  that  minimize  CTj  and  cu.  The  FORTRAN  program  JETFIT  used 
for  this  fitting  procedure  is  listed  in  Appendix  C. 

The  sets  of  measurements  taken  during  the  four  jet  studies  at  the 
Point  Beach  outfall  and  one  of  the  studies  at  the  Palisades  outfall  (Octo¬ 
ber  10,  1972)  were  each  fitted  independently.  The  measurements  of  the 


other  two  Palisades  surveys  did  not  lend  themselves  to  this  fitting  procedure 
because  the  limited  number  of  data  points  did  not  appear  to  characterize  a 
major  portion  of  the  jet  region;  i.e.,  the  jet  seems  to  be  significantly  wider 
than  the  region  surveyed. 

The  parameters  resulting  from  the  fitting  procedure  and  a  tabulation 
of  the  deviations  of  the  fitted  functions  from  the  measurements  are  included 
in  Appendix  D.  The  average  deviation  of  the  function  from  the  data  varied 
from  about  0.2  to  1.1  C°,  with  an  overall  average  deviation  for  all  five  studies 
of  0.7C°.  The  average  velocity  deviation  varied  from  1  to  8  cm/sec,  with  an 
overall  average  deviation  of  4  cm/sec.  These  overall  average  deviations 
correspond  to  about  7%  of  the  average  initial  temperature  excess  and  outfall 
velocity,  respectively.  Figures  8-48  show  the  temperature  and  velocity 
centerlines  and  the  temperature  and  velocity  half -widths  resulting  from  this 
fitting  procedure,  along  with  the  actual  field  measurements.  Note  that  the 
field  measurements  were  made  in  groups  corresponding  approximately  to 
perpendicular  transects  of  the  jet.  Some  of  the  fitting  results  for  the 
Point  Beach  studies  are  shown  together  in  Figs.  49-52  so  that  they  may  be 
compared  with  each  other.  The  circles  on  the  curves  of  0C/  90,  plotted  as 
a  function  of  distance  from  the  outfalls  (Fig.  49),  indicate  the  approximate 
location  of  the  measurement  groupings.  The  results  for  the  one  Palisades 
case  are  presented  in  Sec.  VII  along  with  the  model  comparisons. 

Figure  49  shows  the  centerline  temperature  and  velocity  decays 
resulting  from  the  fits  to  the  measurements  taken  at  a  depth  of  0.5  m  for 
all  four  Point  Beach  dates.  The  drop-off  rates  are  approximately  the  same 
for  each  survey  and  for  both  temperature  and  velocity,  except  for  the  velocity 
results  of  July  13,  1972.  This  difference  is  presently  unexplained.  Indeed, 
the  differences  shown  in  the  figures,  with  that  one  exception,  are  probably 
within  the  errors  associated  with  the  experimental  measurements  and  the 
fitting  procedure.  On  all  four  dates,  the  power  plant  was  operating  at  essen¬ 
tially  the  same  power  level  and  discharge  flow  rate.  The  temperature  dis¬ 
tribution  across  the  outfall  has  been  measured  and  found  to  be  quite  uniform. 
Therefore,  the  excess  temperature  ratio  0C /  90  plotted  in  the  upper  half  of 
Fig.  49  must  be  1.0  at  s  =  0.  This  restriction  was  not  placed  on  the  fitting 
function,  and  so  not  much  significance  should  be  attached  to  the  details  con¬ 
tained  in  the  first  100  m  of  the  fitting  results.  The  initial  velocity  distri¬ 
bution  at  the  Point  Beach  outfall  has  also  been  measured.*  It  was  not  found 
to  be  uniform.  In  fact,  velocities  as  much  as  60%  above  the  average  outfall 
velocity  were  observed.  This  fact  shows  up  in  the  lower  half  of  Fig.  49, 
where  the  velocity  ratio  exceeds  1.0  for  small  values  of  s.  Again,  not  much 
significance  should  be  attached  to  the  details  of  the  results  in  this  initial 
100  -m  region. 

Figure  50  presents  the  corresponding  temperature  and  velocity  half¬ 
widths  for  the  same  four  dates.  It  is  evident  that  the  temperature  distribu¬ 
tion  is  wider  than  the  velocity  distribution.  Indeed,  on  the  average,  it  is 
approximately  twice  as  wide. 

*A  typical  outfall  velocity  distribution  for  Point  Beach  appears  in  ANL/ES-16  (Ref.  19  of  Appendix  A). 

Measurements  on  several  other  dates  exist  and  will  be  published  in  the  future. 
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Fig.  49.  Centerline  Excess-temperature  and  -velocity  Fig.  50.  Half-widths  of  Temperature  and  Velocity  Distri 

Decays  Resulting  from  Fits  to  Point  Beach  Jet  butions  Resulting  from  Fits  to  Point  Beach  Jet 

Data  at  0.5-m  Depth.  ANL  Neg.  No.  190-950.  Data  at  0.5-m  Depth.  ANL  Neg.  No.  190-961. 


71 


(SJ343LU)  J-M  (SJStauj)  nM 


Fig.  51.  Centerline  Temperature  Excess  and  Velocity  Excess  as  aFunc-  Fig.  52.  Half-widths  of  Temperature  and  Velocity  Distributions  as 

tion  of  Depth  Resulting  from  Fits  to  Point  Beach  Jet  Data  a  Function  of  Depth  at  s  =  150  m,  Resulting  from  Fits  to 

(normalized  to  the  0.5-m  results).  ANL  Neg.  No.  190-953.  Point  Beach  Jet  Data.  ANL  Neg.  No.  190-945. 


Figure  51  shows  the  centerline  temperature  and  velocity  excesses  as 
a  function  of  depth  relative  to  the  0.5-m  results.  Because  of  the  nature  of  the 
functional  form  chosen,  these  quantities  are  independent  of  distance  from  the 
outfall  after  the  initial  region,  which  was  always  less  than  about  100  m.  Within 
the  accuracy  of  these  results,  the  excesses  are  fairly  constant  with  depth. 
However,  the  temperature  data  of  September  9,  1972,  do  show  an  exception. 
Because  there  is  no  clear  drop-off  with  depth,  a  characteristic  depth  of  the 
jet  could  not  be  extracted  from  the  JETFIT  results. 

Figure  52  shows  the  widths  as  a  function  of  depth  at  a  distance  of 
150  m  from  the  outfall.  This  distance  was  chosen  because  it  lies  approxi¬ 
mately  in  the  middle  of  the  range  for  which  measurements  were  made.  These 
widths  show  a  decreasing  trend,  tending  toward  zero  at  a  depth  of  2.5 -3.0  m. 
The  temperature  widths  are  approximately  constant  for  the  first  1.0 -1.5  m 
and  then  decrease  while  the  velocity  widths  decrease  more  uniformly. 


In  summary,  the  Point  Beach  jet  studies  indicate  that  the  centerline 
excess  temperature  and  excess  velocity  ratios  decay  at  about  the  same  rate, 
with  jet  longitudinal  distance,  for  all  four  dates  and  for  depths  down  to  2  m 
or  more.  This  decay,  beyond  about  the  first  50  m,  is  characterized  by 


In  the  above  expression,  4.2  is  the  average  value  of  O'  and  g  for  all  cases 
except  for  the  0.5-m  velocity  results  of  July  13,  1972,  and  the  1 . 5  -,  2 .0  -  , 
and  2.5 -m  temperature  results  of  September  9,  1972.  The  average  deviation 
of  the  actual  values  of  O'  and  $  from  this  value  of  4.2  is  0.7.  The  tempera¬ 
ture  half -widths  at  the  0.5-m  depth  grow  at  a  rate  of  3-6  m  for  each  10  m 
from  the  outfall;  the  velocity  half-widths  grow  at  only  about  2  m  per  10  m. 
Both  these  rates  decrease  with  depth,  tending  toward  zero  at  about  2.5 -3.0  m. 
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VI.  MATHEMATICAL  MODELING  OF  NEAR -FIELD  REGION 
OF  SURFACE  THERMAL  DISCHARGES 

The  utility  of  analytical  near -field,  surface -dis charge  models  is  pri¬ 
marily  centered  around  three  considerations.  The  first  of  these  is  a  desire 
to  assess,  and  therefore  avoid,  possible  recirculation  of  heated  water  into 
the  plant  intake.  The  second  is  the  necessity  to  develop  an  appropriate  dis¬ 
charge  design  to  satisfy  temperature  standards.  The  third  is  a  desire  to 
help  predict  possible  biological  effects  relating  to  changes  in  the  physical 
and  chemical  properties  of  the  water. 

Most  presently  operating  power  plants  use  a  surface  channel  or  canal 
to  discharge  their  heated  condenser -cooling  water,  and  many  new  plants  still 
plan  for  similar  discharge  designs.  Since  thermal  discharges  of  cooling  water 
from  power  plants  will  undoubtedly  increase  greatly  in  both  magnitude  and 
number  during  the  next  decade,  the  adequacy  of  predictive  models  for  such 
discharges  becomes  an  important  factor  in  design. 

Table  II  summarizes  the  basic  characteristics  of  15  near-field 
analytical  models1-15  for  surface  thermal  discharges  presently  available 
in  the  literature.  Table  III  summarizes  those  complete -field  models16-24 
that  have  separate  jet -regime  analyses.  Due  to  limitations  of  time  and 
space,  the  field  of  24  models  had  to  be  narrowed  down  considerably.  Con¬ 
sidered  in  this  report  and  compared  with  the  Point  Beach  Unit  1  and  Palisades 
data  are  the  more  promising  and  widely  used  models  of  Pritchard  (Model  No.  1 ), 
Motz  and  Benedict,  Stolzenbach  and  Harleman,  and  Prych.  Some  of  the  reasons 
for  not  including  comparisons  with  the  other  models  are  as  follows: 

1.  Hoopes  et  al.:  model  too  sensitive  to  the  main  parameter,  wind 
speed;  generally  unsatisfactory  comparison  with  earlier  data  taken  at 
Point  Beach  (Ref.  21  of  Appendix  A). 

2.  Hayashi  and  Shuto:  mainly  historical  nature  of  model;  no  jet 
entrainment  simulated  in  this  model. 

3.  Wada  (Model  No.  l)  :  computer  program  unavailable. 

4.  Carter:  mainly  historical  nature  of  model;  jet  model  based  on 
hydraulic  data  taken  in  too  small  a  basin;  improved  data  appeared  recently 
in  Ref.  25. 

5.  Koh  and  Fan  (2D  Model):  no  lateral  entrainment  simulated  in 
this  two-dimensional  (longitudinal,  vertical)  model;  as  with  the  Wada 
Model  No.  2,  it  is  basically  a  tool  to  develop  greater  physical  insight. 

6.  Koh  and  Fan  (axisymmetric  model):  discharge  outfall  assumed 
in  this  model  is  circular,  not  rectangular  as  exists  at  Point  Beach  and 
Palisade  s . 

7.  Barry  and  Hoffman:  model  authors  prefer  to  make  changes  in 
model  development  before  further  application. 


8.  McLay  et  al. :  modification  on  Hoopes  et  al.  model  which  includes 
a  region  of  flow  establishment  and  a  constant  spreading  rate;  not  used  because 
of  generally  unsatisfactory  experience  with  computer  code. 

9.  Stefan:  no  region  of  flow  establishment,  which  severely  restricts 
application  of  model;  also,  model  comparison  with  Stefan's  tank  data  not 
satisfactory  (possibly  due  also  to  stratification  of  ambient  water  in  tank). 

10.  Paul  and  Lick:  new  model  with  computer  code  incomplete  at 
time  of  present  study. 

11.  Engelund  and  Pedersen:  model  applicable  only  for  surface  dis¬ 
charges  having  large  initial  densimetric  Froude  numbers. 

12.  Waldrop  and  Farmer:  model  is  of  too  recent  an  origin. 

13.  Wada  (Model  No.  2):  no  computer  code  available;  no  lateral  en¬ 
trainment  simulated  in  this  two-dimensional  (longitudinal,  vertical)  model. 

14.  Sundaram  et  al. :  mainly  historical  nature  of  model;  jet- 
trajectory  formulas  not  applicable  until  reasonably  distant  from  outfall. 

15.  Elliott  and  Harkness:  model  authors  are  presently  improving 
their  near-field  analysis. 

16.  Giles  et  al.  :  model  authors  will  be  making  further  (though 
minor)  modifications  in  the  model;  the  hybrid  computer  required  for  model 
application  is  too  large  for  Argonne  facilities. 

17.  Loziuk  et  al. :  model  authors  are  presently  modifying  the 
hydraulics  of  their  jet  analysis;  computer  code  presently  unavailable. 

18.  Brady  and  Geyer:  model  is  of  too  recent  an  origin  for  application 
in  this  report. 

19-  Till:  model  is  of  too  recent  an  origin. 

20.  Pritchard  (Model  No.  2):  model  is  of  too  recent  an  origin. 

The  four  analytical  models  to  be  considered  in  this  report  and  com¬ 
pared  with  the  present  field  data  are  summarized  next. 

A.  Motz  -  Benedict  Model 

Motz  and  Benedict  have  developed  a  two-dimensional  model  for  the 
velocity  and  temperature  distribution  of  a  heated  surface  jet.5  Ambient 
receiving-water  turbulence  and  buoyancy  are  assumed  to  be  of  minimal 
influence  with  regard  to  plume  dynamics  and  heat  transfer  within  the  jet 
regime.  The  authors  built  their  integral  analysis  upon  the  framework  es¬ 
tablished  primarily  by  Morton,26  Fan,27  Hoopes  et  al.,1  and  Carter.4  Con¬ 
servation  equations  of  mass,  x  and  y  momentum,  and  energy,  along  with 
two  equations  of  jet  bending,  lead  to  a  system  of  six  ordinary  differential 
equations  which  must  be  solved  numerically  to  yield  the  jet  trajectory, 
centerline  temperature  and  velocity,  and  profile  width.  The  Motz  -  Benedict 
model  assumes  that  only  two  factors  affect  the  flux  of  momentum  at  any 


TABLE  II.  Summary  of  Characteristics  of  Jet  Models 
(ANL  Neg.  No.  190-990  Rev.  1) 
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TABLE  III.  Summary  of  Characteristics  of  Complete -field  Models 

(ANL  Neg.  No.  190-989  Rev.  1) 
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lateral  slice  of  the  jet:  The  first  is  the  entrainment  of  lateral  crossflow 
momentum,  and  the  second  is  a  net  pressure  force  caused  by  eddying  of 
the  ambient  fluid  in  the  lee  of  the  jet  and  by  distortion  of  the  jet  boundaries. 
Explicit  expressions  must  be  assumed  for  both  the  entrainment  function  (to 
determine  mass  and  momentum  conservation)  and  the  drag  force  per  unit 
length  (to  determine  momentum  conservation).  The  assumed  form  of  the 
entrainment  velocity  is 

Vj_  =  E(U  -  Ua  cos  0), 

whe  re 

v^  =  inflow  velocity  (of  entrainment), 

U  =  local  jet  centerline  velocity, 

Ua  =  ambient  current  velocity, 

0  =  local  angle  between  jet  centerline  and  ambient  current, 

and 

E  =  entrainment  coefficient. 


The  form  of  the  assumed  drag  force  is 


^D^az  s^n  ^ 


D 


where 


Fq  =  drag  force  operating  normal  to  jet  axis, 

Cq  =  experimentally  determined  drag  coefficient, 

and 

z  =  jet  depth. 


This  form  for  is  based  on  the  assumption  that  the  interaction  of 

the  jet  with  the  ambient  current  can  be  treated  as  if  the  jet  were  a  solid  body 
and  the  resulting  pressure  gradients  can  be  represented  by  a  drag  force. 
Utilization  of  the  model  requires  the  specification  of  the  drag  and  entrainment 
coefficients,  and  E,  respectively. 

In  the  region  of  established  flow,  the  authors  make  similarity  assump¬ 
tions  for  velocity  and  temperature: 

u(s ,  71)  =  U(s)  exp(-7)2/b2) 


and 


T ' (s ,  T])  =  T (s )  exp(  - T)2/b2) 
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where 

u  =  local  jet  velocity, 

T'  =  local  jet  temperature, 

U  =  jet  centerline  velocity, 

T  =  jet  centerline  temperature, 
s  =  distance  along  jet  centerline, 

T]  =  lateral  distance  from  jet  centerline, 

and 

b  =  characteristic  width  of  lateral  profiles. 

The  assumed  forms  for  u,  T',  v^,  and  Fq  are  used  in  the  integral 
equations  of  conservation  to  yield  the  jet  characteristics  in  the  established 
flow  regime. 

To  account  for  the  zone  of  flow  establishment,  Motz  and  Benedict 
have  developed  phenomenological  relationships  for  the  length  of  the  region 
of  flow  establishment  (s^)  and  the  initial  angle  of  the  jet  at  the  end  of  the 
region  of  flow  establishment  ({30).  These  relationships  are  based  on  a  series 
of  laboratory  hydraulic  studies.  No  details  of  the  temperature  and  velocity 
distribution  in  the  region  of  flow  establishment  were  determined,  only  its 
length  and  final  angle  with  the  shore.  These  studies  yielded  values  for  s^, 
and  E  in  terms  of  the  ratio  A  of  ambient  to  initial  jet  velocity  and 
the  actual  initial  angle  of  the  jet  with  respect  to  shore,  j3g.  (The  authors, 
however,  recommend  a  value  for  of  0.5  be  used  for  most  cases.)  Also 

needed  is  the  width  of  the  jet  at  the  end  of  the  region  of  flow  establishment, 
2b0.  The  authors  suggest  a  value  calculated  by  equating  the  heat  flux  at  the 
point  of  discharge  with  the  heat  flux  at  the  end  of  flow  establishment.  The 
actual  width  at  the  point  of  discharge  is  denoted  by  2bg. 

B .  Stolzenbach-Harleman  Model 

The  mathematical  model  of  Stolzenbach  and  Harleman9  predicts  the 
distributions  of  temperature  and  velocity  within  a  completely  determined  jet 
structure  for  a  near -field  region,  defined  by  the  predominance  of  initial  jet 
momentum  over  the  effects  of  ambient  lake  turbulence. 

Wherever  possible,  the  model  synthesizes  previous  knowledge  of 
buoyant  and  nonbuoyant  jets.  The  heated  discharge  is  assumed  to  be 
structured  in  its  physical  characteristics  as  well  as  its  assumed  velocity 
and  temperature  distribution,  basically  like  a  classical,  turbulent,  non¬ 
buoyant  jet.  Just  as  for  nonbuoyant  jets,  the  authors  assume  an  initial 
core  region  void  of  shear  followed  by  the  main  turbulent  region.  The  basic 
jet  structure  assumed  by  Stolzenbach  and  Harleman  appears  in  Fig.  53. 
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Fig.  53.  Geometrical  Characteristics  of  Jet  Assumed  in  Stolzenbach-Harleman  Model 


Using  an  integral  method  as  in  most  jet  analyses,  velocities  and  tem¬ 
peratures  at  each  longitudinal  cross  section  are  presumed  to  be  related  to 
centerline  values  by  similarity  profiles.  Figure  54  illustrates  the  similarity 
profiles  assumed.  These  are  represented  mathematically  by 

u  =  uc(x)Fy(y)Fz(z)  -  V  cos  0 


and 


AT  -  ATc(x)T  (y)T  (z), 


whe  re 


y 


u,  AT  =  local  value  of  velocity  and  temperature  excess, 
x,y,z  =  longitudinal,  lateral,  and  vertical  coordinates, 
V  cos  0  =  component  of  ambient  velocity  normal  to  jet, 
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y 
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The  form  of  the  similarity  functions  (from  Abramovich  )  are 

f  =  (1  -  G3'2)2, 


and 


t  =  l  -  c3/2. 

The  variables  b,  s,  h,  and  r  are  defined  in  Figs.  53  and  54. 

Horizontal  and  vertical  entrainment  of  ambient  fluid  is  related  to  the 
jet  centerline  velocity  by  appropriate  entrainment  coefficients.  The  lateral- 
entrainment  coefficient  is  determined  by  nonbuoyant  jet  theory  alone  and  is 
constant;  the  vertical -entrainment  coefficient  is  related  to  the  local  tem¬ 
perature  gradient  by  the  local  Richardson  number,  using  the  experimental 
results  of  Ellison  and  Turner,29  and  is  chosen  to  reduce  to  the  nonbuoyant 
value  when  no  density  gradients  exist. 
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Fig.  54.  Velocity  and  Temperature  Characteristics  of 
Jet  Assumed  in  Stolzenbach-Harleman  Model 

Buoyant  convection  is  incorporated  through  the  pres  sure -gradient  terms 
in  the  equations  of  motion  and  through  the  ve rtical -entrainment  coefficient. 
Buoyancy  effects  generally  reduce  vertical  entrainment  and  enhance  lateral 
spreading. 

To  treat  the  interaction  of  these  complex  phenomena  encompassing  jet 
momentum,  entrainment,  buoyancy,  and  ambient  crosscurrent,  the  authors 
developed  their  model  from  the  steady,  time -ave raged  differential  equations  of 
mass,  momentum,  and  conservation  of  heat  energy  by  dropping  negligible  terms, 
assuming  a  form  for  some  of  the  unknown  variables,  and  finally  integrating  the 
simplified  equations  over  the  four  assumed  regions  of  the  jet  longitudinal  cross 
section. 

Boundary  conditions  for  the  differential  equations  of  mass,  momentum, 
and  energy  (along  with  jet  geometry)  required  assumptions  for  the  boundary 
values  of  heat  and  mass  fluxes  as  well  as  for  internal  lateral  and  vertical 
velocity  distributions. 


The  set  of  conservation  differential  equations  for  a  buoyant  deflected 
jet  are  integrated  over  the  jet  cross  section.  The  continuity  and  x-momentum 
equations  are  each  integrated  over  all  of  the  four  regions  defined  in  Fig.  53, 
yielding  eight  equations.  The  y-momentum  equation  is  integrated  over  half 
the  (symmetric)  cross  section  to  yield  the  ninth  equation.  The  tenth  is  a 
jet -bending  equation  obtained  from  equating  the  rate  of  entrainment  of 
lateral  momentum  to  the  rate  of  jet  deflection.  (No  drag  force  on  the  jet 
is  assumed.) 


In  this  manner,  the  coupled  flow  and  energy  equations  reduce  to  a 
system  of  simultaneous,  first-order,  nonlinear  ordinary  differential  equa¬ 
tions  in  the  single  variable  x  (longitudinal  distance  along  jet  centerline 
from  the  outfall),  which  are  then  solved  numerically.  The  solution  to  the 
model  equations  yields  three-dimensional  temperature  and  velocity  pre¬ 
dictions  as  well  as  the  jet  physical  characteristics. 


The  complete  solution  may  be  written  in  nondimensional  form  as 


T  -  T  a  " 
To  -  Ta 


u 

u0 


/  K  V  \ 

>  functions  [  IF0.  A,  - ,  — ,  and 

^  pcpu0  u0  / 


iet  characteristics  = 


whe  re 


IFn  = — ' —  —  •  =  discharge  densimetric  Froude  Number, 

/Ap„ 

y  —  sh° 

A  =  h0/b0  =  dis charge -channel  aspect  ratio, 

K  r 

-  =  surface -heat -los  s  parameter, 

pcpu0 

V  =  ambient  crossflow  velocity  (possibly  a  function  of  offshore 
distance), 

u0  =  channel  outlet  velocity, 

a0  =  angle  of  discharge  with  respect  to  shoreline, 

Ta  =  uniform  lake  water  temperature, 

T0  =  initial  discharge  temperature, 

Apo  =  density  difference  between  discharged  and  ambient  water, 
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and 


pa  =  density  of  ambient  lake  water, 
g  =  acceleration  due  to  gravity, 

2b0  =  width  of  discharge  canal, 
h0  =  depth  of  discharge  canal, 

K  =  surface -heat-transfer  coefficient, 
p  =  density  of  water, 

c  =  specific  heat  of  water. 

P 


Thus  the  only  s ite -dependent  input  parameters  required  are: 


IF q,  A, 


K  V 
pcpu0’  u0 


and  ct0. 


C.  Prych  Model 


1  ? 

The  Prych  model  is  based  on  a  three-dimensional  integral  analysis 
of  a  turbulent,  buoyant,  horizontal,  surface  jet  into  a  large,  deep,  uniform, 
turbulent  flowing  receiving  water.  Integral  equations  of  conservation  are 
written  for  mass,  horizontal  x  and  y  momentum,  and  energy,  as  well  as  equa¬ 
tions  for  the  jet  trajectory.  Figure  55  illustrates  the  coordinate  systems  and 
the  jet  region.  The  model  uses  Gaussian  similarity  assumptions  for  tem¬ 
perature  and  velocity: 

t(s,n,z)  =  T'(s)  exp[ -n1 2/B2(s )]  exp[ -z2 /H2  (s  )] 

u(s,n,z)  =  U(s)  exp[  -n2/B2(s)]  exp[ -z2/H2(s )]  +  V  cos  9(s), 


whe  re 


t,  u  =  local  values  of  temperature  excess  and  velocity, 


s ,  n 


z 

T 1 (s ),  U (s ) 


B (s ),  H(s) 

and 


curvilinear  coordinates  along  jet  centerline  parallel  and 
perpendicular  to  the  jet  trajectory,  respectively, 

vertical  distance  from  receiving -water  surface, 

centerline  values  of  excess  temperature  and  excess 
velocity, 

local  characteristic  width  and  depth  of  jet, 


=  component  of  ambient  velocity  parallel  to  jet  trajectory. 


V  cos  9 
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Fig.  55.  Definition  Sketch  for  Coordinate  System  and  Jet  Region 
Assumed  in  Prych  Model.  ANL  Neg.  No.  190-1045. 

The  rate  of  increase  of  the  local  jet  flow  due  to  incorporation  of  am¬ 
bient  water  is  included  in  the  model  by  simulating  entrainment  due  to  jet  mix¬ 
ing  and  also  entrainment  caused  by  turbulence  in  the  ambient  fluid.  Horizontal 
and  vertical  rates  of  entrainment  are  calculated  separately  for  both  processes 
and  are  then  summed.  Prych  computes  jet  entrainment  in  a  manner  similar  to 
Stolzenbach  and  Harleman.  Most  notable  is  the  requirement  that  the  entrain¬ 
ment  coefficient  decreases  with  increasing  fluid  stability  in  the  density- 
stratified  jet.  Prych  calculates  the  contribution  of  ambient  mixing  to  horizontal 
and  vertical  entrainment  as  3.5 e/a,  where  e  is  the  ambient  turbulent  diffusion 
coefficient  (horizontal  or  vertical)  and  a  is  the  standard  deviation  of  the  lateral 
distribution  of  a  tracer  (excess  temperature)  in  a  two-dimensional  jet.  Prych 
assumes  a  =  2  s/H  for  the  vertical  direction. 

Four  force  terms  appear  in  the  x-  and  y -momentum  equations: 

1.  Ambient  momentum  flux  in  the  direction  of  the  ambient  current. 

2.  Pressures  due  to  the  density  differential  between  the  fluid  in  the 
jet  region  and  the  ambient  fluid.  This  force  is  assumed  hydrostatic  and  is 
calculated  by  integrating  the  excess  pressure  force: 


APo 

g  —  (s ,  n,  z 
1  o 


)  dz  1 


over  each  vertical  cross-sectional  face. 


3.  A  form  drag  on  the  jet  due  to  the  pressure  difference  between 
offshore  and  lee  sides  of  the  jet,  represented  by 

Fd  =  Hv|v|  sin2  0, 

where  =  form-drag  coefficient. 

This  force  is  computed  in  the  same  way  one  computes  form  drag 
on  a  solid  body. 
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4.  Interfacial  shear  stress  between  the  jet  fluid  and  the  underlying 
ambient  fluid.  Prych  approximates  this  force  by  modifying  an  expression  for 
the  shear  stress  at  the  base  of  a  turbulent  boundary  layer  on  a  smooth,  flat, 
solid  surface. 


Once  the  basic  conservation  equations  are  integrated  over  the  jet  cross 
section  using  the  above  assumptions,  a  system  of  ordinary  differential  equa¬ 
tions  results  that  are  solved  for  jet  trajectory,  jet  width,  jet  depth,  and  jet- 
centerline  temperature  and  velocity.  A  separate  region  of  flow  establishment 
is  included  with  the  same  forces  simulated  as  described  previously. 

The  basic  input  parameters  to  the  model  are: 


and 


Q0  =  discharge  flow  rate  from  the  outfall, 

T0  =  excess  temperature  of  jet  (above  ambient)  at  the  outfall, 

0O  =  angle  of  discharge  velocity  vector  with  positive  x  axis, 

2b0  =  width  of  outfall, 
h0  =  depth  of  outfall, 

V  =  ambient -current  velocity  (assumed  constant), 

AP0  =  difference  in  density  between  ambient  water  and  water  from 
outlet, 

K  =  surface -heat -transfer  coefficient, 

E0  =  jet -entrainment  coefficient  (horizontal), 

Cq  =  form-drag  coefficient, 

Cp  =  (interfacial)  shear-stress  coefficient. 


D.  Pritchard  Model 


t  1  7 

Pritchard's  model  is  basically  a  synthesis  of  previous  theoretical  - 
and  phys  ical -modeling  results  for  buoyant  and  nonbuoyant  jets,  complemented 
by  results  the  author  has  gleaned  from  field  data  obtained  from  existing  power 
plants  sited  on  bays,  estuaries,  and  large  lakes. 

The  model  is  simple  and  considers  plume  dispersion  to  be  governed 
solely  by  momentum-jet  entrainment,  turbulent  diffusion,  and  surface  heat 
loss.  An  integral  technique  is  used  in  which  the  plume  velocity  and  excess 
temperature  (above  ambient)  are  assumed  to  have,  at  each  longitudinal  posi¬ 
tion,  a  "tophat"  distribution  laterally  and  vertically.  Buoyancy -induced 
convective  motions  are  not  explicitly  considered.  No  ambient  current  is 
assumed  to  exist  in  the  theoretical  development;  yet,  the  author  expects 
that  the  predicted  cente rline -tempe rature  decay  and  areas  within  isotherms 


will  still  be  accurate  in  the  presence  of  a  current  less  than  10%  of  the  initial 
discharge  velocity.  In  the  model,  environmental  changes  are  reflected  solely 
in  the  surface -heat -transfer  coefficient,  K.  Entrainment  is  accounted  for  by 
the  specification  of  a  fixed  inverse  spreading -rate  parameter,  n,  normally 
taken  to  be  about  6.  The  lake  bottom  is  assumed  to  have  no  effect  other  than 
upon  the  author's  choice  of  a  plume  depth  and,  when  necessary,  upon  the 
length  and  depth  of  an  initial  region  of  vertical  entrainment. 

The  model  handles  both  the  jet  and  far -field  regions;  it  predicts  a 
two-dimensional  temperature  field  and  areas  within  isotherms  down  to  a 
0.56C°  (l.OF0)  temperature  excess.  The  parameters  required  for  the  appli¬ 
cation  of  this  model  are: 

b0  =  width  of  rectangular  outfall, 
h0  =  depth  of  rectangular  outfall, 

Qjl*  =  excess  heat-rejection  rate  of  power  plant  based  on  90, 

90  =  initial  excess  temperature  of  the  jet, 

and 

K  =  surface -heat -loss  coefficient. 

In  spite  of  the  simplifications  made  in  the  model  development,  the 
author  claims  the  model  to  be  conservative  in  many  respects  and  simple  to 
apply. 


The  author's  theoretical  development  is  carried  out  in  four  consecutive 

stages : 


1.  Horizontal  spreading  is  considered,  neglecting  vertical  diffusion 
and  surface  heat  loss  to  the  atmosphere.  A  two-dimensional  temperature 
field  is  determined  for  the  jet  and  far-field  regions  by  the  integral  equations 
of  conservation.  Centerline  temperature  and  velocity  is  found  to  drop  off  as 
1  /  VW  in  the  jet  regime,  where  s  is  the  distance  along  the  centerline  after  a 
constant  region  of  flow  establishment  of  length  6b0.  For  the  far-field  region, 
excess  temperature  is  assumed  to  drop  off  as  l/s. 

2.  Vertical  entrainment  is  then  considered  independent  of  horizontal 
spreading  and  surface  heat  loss.  The  depth  is  assumed  constant,  except  pos¬ 
sibly  for  a  small  region  in  the  vicinity  of  the  outfall,  where  the  depth  grows 
slowly  in  a  linear  fashion.  When  vertical  spreading  is  assumed  to  occur,  a 
correction  of  the  two-dimensional  temperature  field  is  made  to  account  for 
the  additional  dilution.  For  Lake  Michigan,  Pritchard  suggests  that  vertical 


*Qj_j  would  be  identically  equal  to  the  total  condenser  heat-rejection  rate  if  the  condenser  intake  temperature 
were  identically  equal  to  the  ambient  temperature. 


spreading  be  allowed  for  until  the  plume  depth  reaches  3.05  m  (10  ft).  For  a 
greater  initial  depth,  the  jet  is  assumed  to  remain  at  that  constant  depth 
throughout  the  complete  field. 

3.  A  second  correction  of  the  temperature  field  may  be  needed, 
depending  on  the  temperature  of  the  water  entrained  into  the  plume.  Due 
to  possible  recirculation  of  condenser  cooling  water,  the  diluting  water 
mixed  into  the  plume  may  have  an  excess  temperature  above  ambient. 

Once  this  additional  correction  on  the  temperature  field  is  made,  the  areas 
within  isotherms  can  then  be  calculated  for  the  condition  of  mixing  alone. 

4.  Surface  heat  losses  are  then  included  in  the  analysis  as  a  cor¬ 
rection  to  the  areas  derived  in  stage  3.  This  surface -heat-loss  correction 
yields  the  final  two-dimensional  temperature  field  and  the  isotherm  areas. 


VII.  MODEL  COMPARISONS  TO  DATA 


Table  IV  summarizes  the  basic  discharge  and  environmental  parame¬ 
ters  required  for  application  of  the  analytical  models  to  the  four  Point  Beach 
cases  and  the  single  Palisades  case.  Parameters  such  as  entrainment  and 
drag  coefficients  were  chosen,  based  upon  the  recommendations  of  the  model 
authors.  In  some  cases  (especially  with  the  Motz-Benedict  model),  no  clear- 
cut  choice  exists  for  some  of  the  required  parameters.  This  problem  is 
discussed  in  more  detail  in  Sec.  VIII. C. 2  below. 

The  data,  JETFIT  -  smoothed  results  and  model  predictions  are  dis¬ 
cussed  and  compared  with  respect  to  some  of  the  major  jet  characteristics 
in  the  following  paragraphs. 

A.  Jet  Trajectories 


Figures  56-60  show  the  jet  trajectories  resulting  from  JETFIT  and 
the  four  model  predictions.  Before  the  analytical  models  are  compared  to 
the  jet  trajectory  data,  a  few  comments  should  be  made  regarding  the  JETFIT 
trajectories  themselves.  Figure  56  shows  the  jet  centerline  to  bend  gradually 
toward  the  north-shore  direction  (negative  x  direction).  In  a  similar  fashion, 
Fig.  57  shows  the  jet  to  bend  gradually  toward  the  south-shore  direction 
(positive  x  direction).  In  both  cases,  the  ambient  current  has  been  assigned 
a  nominal  value  of  zero,  based  on  field  measurements.  Possible  causes  of 
the  seemingly  anomalous  behavior  may  be  combinations  of  the  following: 

1.  A  small  but  undetectable  ambient  current,  directed  north  on 
May  18  and  south  on  May  23,  may  have  been  present.  The  threshold  of  the 
instrument  used  for  measuring  ambient  current  is  about  3.0  cm/sec;  con¬ 
sequently,  any  current  below  or  near  this  value  is  virtually  undetectable. 
Further,  the  ambient-current  conditions  were  measured  at  a  limited  number 
of  specific  locations.  A  small  current  could  possibly  have  been  detected  at 
some  other  location,  since  near-shore  current  measurements  are  known  to 
be  spatially  and  temporally  unsteady. 

2.  Local  gyres  might  exist  in  the  region  of  the  outfall.  The  intru¬ 
sion  of  the  33-m  discharge  structure  into  a  small  ambient- current  field  will 
cause  local  eddies  and  gyres  to  form  in  the  vicinity  of  the  discharge.  Data 
taken  near  the  outfall  may  be  influenced  by  them.  Also,  a  return  current  on 
the  lee  side  of  the  jet  (when  an  ambient  current  exists)  or  on  both  sides  of  the 
jet  (when  no  ambient  current  is  present)  is  required  to  provide  water  to  the 
jet  as  it  disperses  due  to  jet  entrainment  and  ambient  mixing.  When  an 
ambient  current  exists,  the  region  between  the  outer  boundary  of  the  jet  on 
the  lee  side  and  the  shoreline  will  be  likely  to  contain  eddies  of  continually 
recirculating  water. 

3.  Lake-bottom  irregularities  in  the  vicinity  of  the  discharge  might 
influence  the  jet  trajectory.  The  contours  in  Fig.  6  represent  average  values 
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Procedure  and  Model  Calculations  for  Point  Beach:  Procedure  and  Model  Calculations  for  Point  Beach: 

May  18,  1972.  ANL  Neg.  No.  190-935.  May  23,  1972.  ANL  Neg.  No.  190-936. 


Fig.  58.  Centerline  Trajectories  Resulting  from  the  Fitting  Fig.  59.  Centerline  Trajectories  Resulting  from  the  Fitting 

Procedure  and  Model  Calculations  for  Point  Beach:  Procedure  and  Model  Calculations  for  Point  Beach: 

July  13,  1972.  ANL  Neg.  No.  190-940.  September  9,  1972.  ANL  Neg.  No.  190-948. 
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Fig.  60.  Centerline  Trajectories  Resulting  from  the  Fitting 
Procedure  and  Model  Calculations  for  Palisades: 
October  10,  1972.  ANL  Neg.  No.  190-957. 


of  bottom  soundings  taken  at 
several  different  dates  in  1972. 
Continual  operation  of  the 
Point  Beach  Unit  1  discharge  has 
apparently  dug  out  a  small  channel 
directly  off  the  discharge.  A  small 
ridge  appears  north  of  the  discharge 
canal,  as  well  as  a  small  bar  just 
south  of  the  outfall.  Changes  in  the 
local  bottom  contours  are  possible, 
due  to  occasional  storms  in  the  area 
and  daily  variation  in  both  ambient 
current  and  discharge  conditions. 

4.  The  fitting  of  symmetric 
profiles  of  temperature  and  velocity, 
as  used  in  JETFIT,  to  distributions 
that  are  probably  unsymmetric  when 
an  ambient  current  exists,  may  in¬ 
fluence  the  apparent  location  of  the 
centerline.  Figure  61  illustrates 
a  somewhat  more  realistic  profile 
of  centerline  excess  temperature 
©/ ©o  and  jet- centerline -velocity 
ratio  u/u0  (in  the  lateral  direction) 
for  a  bent  jet  near  the  outfall.  The 
profiles  are  based  on  hydraulic- 


model  measurements  by  Carter25  for  an  outfall  angle  P0  of  60°  and  a  ratio  of 
outfall  to  ambient  current,  R,  of  2.0.  The  JETFIT  data  discussed  here  are 
for  R  =  10  to  00  (i.e.,  weak  currents),  which  should  reduce  but  not  eliminate 
the  asymmetry  of  the  velocity  and  temperature  profiles  shown  in  the  figure. 
Any  asymmetry  due  to  an  ambient  current  should  bias  the  centerline  fit  to  the 


Fig.  61 

Idealized  Surface  Profile  of  Excess  Temperature  0/6^, 
and  Velocity  u/uq  across  (y  direction)  Bent  Jet  near 
Orifice  (based  on  measurements).  R  =  2  (nominal) 
and  Sq  =  60°  (Ref.  25). 


lee  direction  for  the  temperature  centerline  and  to  the  offshore  direction  for 
the  velocity  centerline.  This  might  explain  the  relative  orientation  of  velocity 
and  temperature  centerlines  for  May  23,  July  13,  and  September  9 
Point  Beach  data  (Figs.  57-59).  The  velocity  centerline  is  also  upstream  of 
the  temperature  centerline  for  the  fitting  results  at  the  other  depths  on  these 
dates;  discrepancies  when  occurring  were  very  small  and  could,  however, 
just  as  easily  have  been  caused  by  insufficient  data.  Far-field  temperature 
measurements  taken  several  hours  later  at  1645-1  751  hours  on  May  18  (not 
presented  here)  indicate  a  southerly  directed  plume  under  the  influence  of 
a  small  ambient  current.  The  results  of  May  18  seem  to  be  largely  biased 
by  data  measurements  at  Station  No.  5,  where  the  plume  appears  to  be 
directed  northward.  There  is  the  possibility  of  a  small  wind-induced  cur¬ 
rent  directed  north  existing  near  the  third  transect  offshore  in  opposition  to 
a  southerly  moving  ambient  current.  Also ,  note  that  Station  No.  7  with  zero  ve¬ 
locity  was  eliminated  from  the  JETFIT  calculation,  since  it  was  not  thought 
to  be  a  point  in  the  thermal  plume.  This  would  bias  the  velocity  centerline 
in  a  southerly  direction  on  May  18.  More  data  are  necessary  to  more  ac¬ 
curately  define  the  centerlines  in  this  case. 

5.  Uncertainties  in  the  data  (up  to  ±0.5C'J  for  temperature,  ±20%  for 
velocity  magnitude,  ±12°  of  velocity  direction)  and  the  limited  amount  of  data 
available  clearly  influence  the  fitting  results.  The  temperatures  recorded 
were  averaged  over  1  -  to  2-min  time  intervals;  thermal  fronts  have  been 
observed  (Ref.  30)  that  can  change  the  plume  temperature  at  a  fixed  point  at 
Point  Beach  up  to  3C  in  1-2  sec.  The  appearance  and  frequency  of  such 
fronts  (estimated  to  be  from  30  sec  to  several  minutes)  will  be  important  in 
any  given  averaged  measurement.  Determination  of  ambient  temperatures 
and  currents  is  also  important  to  the  analysis,  yet  are  difficult  to  determine 
in  the  field. 

Several  of  the  general  comments  given  in  paragraphs  2-5  above  also 
apply  to  the  jet  data  of  July  13  and  September  9.  On  September  9,  Unit  2  at 
Point  Beach  was  operating  at  12%  power  with  only  one  of  the  two  available 
pumps.  This  corresponds  to  a  flow  rate  of  about  14  m3/sec  with  a  tempera¬ 
ture  excess  above  ambient  of  about  2C°.  Clearly,  the  velocity  distribution 
should  be  affected  by  this  second  discharge,  with  the  temperature  field  per¬ 
turbed  only  slightly.  The  nonlinear  interaction  of  flow  and  temperature  fields 
complicates  any  further  evaluation  of  resultant  effects. 

Analytical  model  predictions  for  May  18  and  May  23  at  Point  Beach 
were  based  on  a  zero  -  ambient- current  value.  Consequently,  a  straight- 
line  trajectory  was  predicted  by  all  models.  The  JETFIT  results  for  these 
dates  yielded  only  slight  differences  from  a  straight-line  trajectory,  prob¬ 
ably  due  to  the  reasons  given  previously.  Models  will  be  compared  to  the 
temperature  centerline  of  the  jet,  since  it  is  generally  of  greater  interest 
than  the  velocity  centerline. 


The  Stolzenbach-Harleman  model  and  the  Prych  model  appear  most 
accurate  for  the  dates  with  a  measured  ambient  current.  Prych  predicts 
greater  bending  than  Stolzenbach  and  Harleman  for  both  dates.  The  Pritchard 
model  suggests  only  a  small  change  in  trajectory  from  a  straight  line  for 
situations  in  which  a  small  ambient  current  exists.  Consequently,  a  straight- 
line  trajectory  is  shown  for  Pritchard  in  Figs.  58  and  59-  The  Motz -Benedict 
model  overpredicts  jet  bending  for  both  dates.  As  stated  above,  the  model  is 
very  sensitive  to  the  value  of  the  entrainment  coefficient  chosen. 

Interpreting  the  Palisades  results  is  difficult  for  several  reasons. 
First,  the  lake  bottom  in  the  vicinity  of  the  discharge  is  very  shallow,  as  may 
be  seen  in  Fig.  7.  This  may  inhibit  any  significant  vertical  entrainment.  For 
the  first  100  m,  the  lake  bottom  is  at  approximately  the  depth  of  the  outfall. 

At  the  north  end  (positive  y-direction)  of  the  diverging  Palisades  outfall,  lake 
water  occasionally  enters  the  discharge  channel.  The  discharge  cross  section 
is  highly  irregular,  as  sketched  in  Fig.  62.  The  discharge  width,  b0,  was 
chosen  as  28.3  m,  since  it  is  the  channel  width  at  the  lake  opening.  The  dis¬ 
charge  depth  h0  was  chosen  to  be 


hn  = 


discharge  area 


the  calculated  h0  was  then  averaged  for  the  three  dates  shown  in  Fig.  62. 

X  (meters) 


•  OCT  10,  1972 

Fig.  62.  Bottom  Depth  at  Palisades  Outfall  for  Three  Dates  of  Jet  Studies 


Considering  the  significant  bottom  interaction  that  undoubtedly  occurs 
in  this  case,  all  model  predictions  will  be  considered  qualitatively  only.  Zero 
ambient  current  on  October  10  clearly  implies  no  bending  for  model  predic¬ 
tions;  JETFIT  results  indicate  little  deviation  from  this  position  (see  Fig.  60). 


B.  Centerline  Temperature  Decay  and  Temperature  Half-widths 


The  centerline  excess  temperatures  resulting  from  JETFIT  and  the 
four  model  predictions  are  shown  in  Figs.  63-67.  The  assumed  parametric 
temperature  profile  to  which  the  jet  data  was  fitted  with  JETFIT  is  again 
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and  A,  Ot ,  C,  and  Y  are  fitted  parameters.  It  is  clear  from  Figs.  63-67  that 
the  temperature  excess  near  the  outfall  from  the  data- smoothing  procedure 
can  be  different  from  0O.  This  corresponds  to  A  /  1. 


The  values  of  A,  Ot ,  and  length  of  region  of  flow  establishment 
s0  =  Ot / A2  (in  units  of  the  outfall  width)  for  the  four  Point  Beach  dates  con¬ 
sidered  are: 


May  18 

May  23 

July  13 

September  9 

A 

0.84 

0.84 

0.92 

0.79 

ot 

3.8 

5. 6 

6.2 

4.5 

So 

5.38 

7.95 

7.32 

7.22 

The  values  of  s0  do  compare  with  values  in  the  range  from  5  to  7  usually 
quoted  in  the  literature.  However,  the  above  Point  Beach  numbers  certainly 
reflect  insufficiencies  in  the  data  and  inadequacies  in  the  fitting  procedure; 
consequently,  much  significance  should  not  be  attached  to  them. 

Since  0  ^  0O  (i.e.,  A  4  l)  at  the  outfall,  one  should  not  place  much 
significance  on  the  results  of  the  fitting  procedure  for  the  first  50-100  m. 

It  may  have  been  more  profitable  to  have  A  fixed  at  1.0  and  introduce  a  new 
parameter  e  as  the  power  of  ocb0/s  in  the  formula  for  0c(s).  (The  present 
formulation  restricts  e  to  a  value  of  1  /  Z . )  The  fit  would  then  involve  ot ,  e, 

C,  and  Y-  It  is  not  expected,  however,  that  the  present  JETFIT  results  would 
be  altered  significantly  if  such  changes  were  implemented  in  the  fitting  func¬ 
tion.  One  further  comment:  The  results  of  JETFIT  for  the  jet  character¬ 
istics  most  generally  agreed  with  what  one  would  expect  from  visual 
examination  of  the  data  along  each  transect  of  the  jet.  This  gives  additional 
confidence  to  the  JETFIT  results. 
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Fig.  63.  Centerline  Temperature  Excess  and  Velocity  Decays  Result-  Fig.  64.  Centerline  Temperature  Excess  and  Velocity  Decays  Result¬ 
ing  from  the  Fitting  Procedure  and  Model  Calculations  for  ing  from  the  Fitting  Procedure  and  Model  Calculations  for 

Point  Beach:  May  18,  1972.  ANL  Neg.  No.  190-952.  Point  Beach:  May  23,  1972.  ANL  Neg.  No.  190-949. 
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Fig.  65.  Centerline  Temperature  Excess  and  Velocity  Decays  Result-  Fig.  66.  Centerline  Temperature  Excess  and  Velocity  Decays  Result¬ 
ing  from  the  Fitting  Procedure  and  Model  Calculations  for  ing  from  the  Fitting  Procedure  and  Model  Calculations  for 

Point  Beach:  July  13,  1972.  ANL  Neg.  No.  190-958.  Point  Beach:  September  9,  1972.  ANL  Neg.  No.  190-951. 
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Fig.  67.  Centerline  Temperature  Excess  and  Velocity  Decays  Result¬ 
ing  from  the  Fitting  Procedure  and  Model  Calculations  for 
Palisades:  October  10,  1972.  ANL  Neg.  No.  190-960. 

The  Pritchard  model  appears  most  accurate  and  also  conservative  for 
the  four  Point  Beach  cases  of  Figs.  63-66.  Both  JETFIT  and  the  Pritchard 
model  have  centerline  temperature  decay  rates  proportional  to  the  -l/l  power 
of  the  centerline  distance;  hence,  both  curves  become  parallel  as  s  increases. 
The  Motz  -  Benedict  model  gives  reasonably  conservative  predictions  for  the 
no- current  cases  of  May  18  and  May  23  where  the  entrainment  coefficient  was 
chosen  to  be  0.05.  (Benedict  recommends  0.04  for  a  90°  outfall  and  suggests  a 
slight  increase  for  an  off-90°  discharge.)  The  model  does  poorly  on  July  13 
and  September  9,  apparently  due  to  the  choice  of  entrainment  coefficient.  The 
authors'  recommendation  for  the  value  of  the  entrainment  coefficient  for  cases 
with  ambient  current  less  than  0.2u0  is  not  clear.  The  model  overpredicts 
temperature  decay  for  both  these  cases.  The  Stolzenbach-Harleman  and  Prych 
models  both  predict  a  much  greater  temperature  decay  than  the  data  indicate. 
When  no  current  exists,  the  Stolzenbach-Harleman  decay  is  initially  very  abrupt, 
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yet  tends  to  level  off  after  about  125  m;  the  centerline  temperature  data 
eventually  drop  below  the  model  results.  For  the  current  cases  of  July  13 
and  September  9,  the  centerline  decay  of  Stolzenbach  and  Harleman  is  more 
regular,  yet  is  too  rapid.  The  Prych  model  predicts  too  great  a  temperature 
decay,  yet  does  so  at  a  rather  regular  rate  of  decrease,  seemingly  independent 
of  an  ambient  current.  Moreover,  Prych  predicts  lower  temperatures  than 
Stolzenbach  and  Harleman,  at  least  within  the  first  200-250  m  from  the  outfall. 
Beyond  this  point,  the  two  models  run  nearly  parallel. 

Figures  68-72  show  the  widths  of  the  excess  temperature  distribu¬ 
tions  resulting  from  JETFIT  and  the  model  calculation.  (Note  that  the  linearity 
of  the  JETFIT  results  is  a  consequence  of  the  linear  form  assumed  in  the 
fitting  function.)  A  very  rapid  rate  of  lateral  spreading  is  apparent  for  the 
Prych  and  Stolzenbach- Harleman  models.  Both  models  require  an  assumed 
form  for  a  lateral  spreading  velocity  that  is  instrumental  in  determining  the 


Fig.  68.  Half-widths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  the  Fitting  Procedure  and  Model  Calculations  for 
Point  Beach:  May  18,  1972.  ANL  Neg.  No.  190-947. 
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Fig.  69.  Half-widths  of  Temperature  and  Velocity  Distributions  Result-  Fig.  70.  Half-widths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  the  Fitting  Procedure  and  Model  Calculations  for  ing  from  the  Fitting  Procedure  and  Model  Calculations  for 

Point  Beach:  May  23,  1972.  ANL  Neg.  No.  190-946.  Point  Beach:  July  13,  1972.  ANL  Neg.  No.  190-942. 
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Fig.  71.  Half-widths  of  Temperature  and  Velocity  Distributions  Result-  Fig.  72.  Half-widths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  the  Fitting  Procedure  and  Model  Calculations  for  ing  from  the  Fitting  Procedure  and  Model  Calculations  for 

Point  Beach:  September  9.  1972.  ANL  Neg.  No.  190-941.  Palisades:  October  10,  1972.  ANL  Neg.  No.  190-959. 
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lateral  extent  of  the  jet.  Both  models,  on  the  basis  of  present  results,  used 
assumptions  that  were  apparently  too  theoretical  in  nature  and  apparently 
require  calibration  with  some  empirical  data  to  yield  more  accurate  results. 
The  overextended  lateral  spread  of  these  two  jet  solutions  is  a  probable 
cause  of  the  extremely  rapid  centerline  temperature  decay  required  by  energy 
conservation  at  each  jet  cross  section. 

The  Motz-Benedict  model  predicts  a  jet  that  is  by  far  too  narrow  for 
all  four  Point  Beach  cases.  The  model  is  very  sensitive  to  the  particular 
choice  for  the  value  of  some  of  the  input  parameters,  in  particular,  the 
entrainment  coefficient  and  the  jet  width  at  the  end  of  the  region  of  flow 
establishment.  A  more  judicious  choice  for  these  parameters  would  un¬ 
doubtedly  improve  the  comparisons.  However,  when  this  model  is  used  as 
a  predictive  tool,  it  is  difficult,  at  present,  to  know  what  values  to  choose 
other  than  those  recommended  by  the  authors  of  the  model. 

Pritchard  predicts  rather  accurate  temperature  decays  and  widths, 
except  for  the  temperature  width  on  September  9-  The  wider  width  of  the 
data  might  well  be  due  to  the  second  unit  operating  at  12%  power  on  that  date, 
adding  to  the  temperature  excesses  on  the  offshore  side  of  the  plume. 

The  Palisades  results  are  quite  interesting  (Figs.  67  and  72).  As 
expected,  the  temperature  does  not  drop  off  until  about  175  m  downstream, 
due  to  the  larger  region  of  flow  establishment  resulting  from  the  wider  out¬ 
fall  at  Palisades.  Again,  Prych  and  Stolzenbach  and  Harleman  are  still  too 
optimistic  in  their  temperature-decay  predictions  and  too  wide  in  their  lateral 
jet  spreading.  Pritchard  predicts  too  great  a  temperature  decrease  and  too 
narrow  a  jet.  This  is  probably  due  to  his  method  of  handling  vertical  entrain¬ 
ment,  which  is  very  sensitive  to  the  difference  between  the  outfall  depth  and 
the  critical  mixing  depth  of  the  lake.  Vertical  entrainment  is  not  a  factor  in 
Pritchard's  predictions  for  Point  Beach.  Motz  and  Benedict  are  conserva¬ 
tive  in  temperature  decay  for  Palisades  due  to  the  large  region  of  flow  estab¬ 
lishment  (5.2  times  the  full  width  of  the  outfall).  Lateral  widths  are  again 
very  much  underpredicted. 

C.  Centerline  Velocity  Decay  and  Velocity  Half-widths 


Figures  63-67,  for  uc/u0  derived  from  JETFIT,  indicate  a  significant 
deviation  of  the  initial  jet  velocity  from  the  calculated  average  channel  ve¬ 
locity.  For  July  13,  velocity  measurements  near  the  channel  outlet  indeed 
verified  a  velocity  greater  than  u0  for  the  first  meter  depth.  Other  measure¬ 
ments,  on  different  dates,  all  showed  values  greater  than  u0. 

The  Stolzenbach-Harleman  and  Prych  models  predict  too  rapid  a 
velocity- centerline  decay.  In  each  case,  the  Stolzenbach-Harleman  model 
predicts  an  increase  in  velocity  after  the  jet  leaves  the  outlet  up  to  60  m 
along  the  centerline.  This  is  possible  for  discharges  having  very  low  den- 
simetric  Froude  numbers.  A  light  fluid,  discharged  over  a  heavier  fluid, 
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accelerates  laterally  as  well  as  longitudinally  due  to  hydrostatic  pressure 
gradients.  For  low  densimetric  Froude  numbers  (near  1.0)  such  buoyant 
accelerations  become  dominant.  Although  increases  in  velocity  are  possible, 
increases  in  temperature  are  never  realistic.  If  such  hydrostatic  pressure 
gradients  were  removed  or  neglected  in  the  model,  a  momentum  jet  would 
occur.  For  discharges  having  high  densimetric  Froude  numbers,  lateral 
accelerations  greatly  dominate  longitudinal  ones  and  consequently  allows 
such  longitudinal  gradients  to  be  neglected.  The  densimetric  Froude  num¬ 
bers  for  the  four  Point  Beach  cases  are  about  Z.4,  which  may  be  a  little 
large  for  the  above  floating -plume  phenomenon  to  occur.  For  the  stagnant- 
lake  cases  of  Figs.  63  and  64,  the  Prych  model  has  a  sharp  decrease  in 
velocity  within  1Z5  m  of  the  outfall  and  then  tends  to  level  off.  The 
Stolzenbach-Harleman  model  is  more  regular  in  its  rapid  velocity  decline, 
crossing  the  Prych  curve  for  uc/u0  at  about  170  m  from  the  outfall.  The 
ambient- current  cases  realize  a  more  regular  velocity  decline,  with  the 
Prych  model  giving  consistently  lower  velocities. 

The  Motz-Benedict  and  Pritchard  models  appear  adequate  for  both 
stagnant- lake  cases  of  May  18  and  May  Z3,  and  the  current  case  of  September  9- 
The  data  of  July  13  reveal  higher  velocities  than  predicted  by  any  of  the  models. 
An  apparent  defect  in  the  Motz-Benedict  model  is  the  assumption  of  decay  of 
lateral  velocity  to  zero,  even  where  an  ambient  current  exists.  A  more 
realistic  assumption  would  require  a  decay  superimposed  on  Ua  cos  P  (the 
component  of  the  ambient  current  parallel  to  the  local  jet  centerline). 

Velocity  widths  as  predicted  by  the  Prych  and  Stolzenbach-Harleman 
models  are  much  too  large  for  the  four  Point  Beach  dates  (Figs.  68-71).  This 
again  reflects  on  their  assumptions  of  lateral  spreading  velocity.  The  Prych 
predictions  are  consistently  larger  than  those  of  Stolzenbach  and  Harleman. 

The  Stolzenbach-Harleman  predictions  appear  too  sensitive  to  ambient  cur¬ 
rent;  the  Prych  predictions  seem  very  insensitive.  Pritchard  has  no  predic¬ 
tions  for  velocity  width.  Examining  the  data  reveals  that  the  temperature 
width  is  approximately  twice  the  size  of  the  velocity  width  for  any  fixed  date 
at  Point  Beach. 

The  Palisades  data  indicate  a  very  slow  velocity  decline  after  a  rather 
rapid  drop  in  the  first  50  m  (see  Fig.  67).  This  rapid  drop  might  be  due  to 
the  long,  diverging  discharge  channel.  Irregularities  in  the  velocity  distribu¬ 
tion  across  the  channel  (lake  water  entering  the  channel  at  the  north  end 
combined  with  significant  channel-depth  irregularities)  may  be  a  contributing 
factor.  No  model  for  velocity  decay  appears  adequate  for  Palisades  on 
October  10.  Velocity  widths  of  Prych  and  Stolzenbach  and  Harleman  are 
again  too  large,  and  the  Motz-Benedict  predictions,  although  more  reasonable, 
are  too  small  (see  Fig.  7Z).  Due  to  site-dependent  irregularities  (diverging 
outfall,  shallow  bottom),  Palisades  is  not  the  ideal  site  to  evaluate  near-field 
models,  albeit  a  real  site. 
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D.  Temperature  and  Velocity  Half-depth 

Figures  73-76  compare  the  predictive  models  with  respect  to  half-depth. 
The  half-depth  is  the  vertical  distance  at  which  the  excess  temperature  and/or 
excess  velocity  drop  to  one-half  the  local  surface- centerline  value.  Also  shown 
is  an  indication  of  the  location  of  the  lake  bottom.  No  data  are  plotted,  since 
the  vertical  profile  information  was  insufficient  to  define  a  half-depth  in  most 
cases.  Either  the  3.0-m  limit  or  the  lake  bottom  was  reached  before  the  half¬ 
depth  or  the  temperature  or  velocity  distribution  was  reached. 
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Fig.  73.  Half-depths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  Model  Calculations  for  Point  Beach:  May  18,  1972 
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Fig.  74.  Half-depths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  Model  Calculations  for  Point  Beach:  May  23,  1972 
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Fig.  75.  Half-depths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  Model  Calculations  for  Point  Beach:  July  13,  1972 
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Fig.  76.  Half-depths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  Model  Calculations  for  Point  Beach:  September  9,  1972 


The  Motz-Benedict  and  Pritchard  models,  being  two-dimensional,  natu¬ 
rally  predict  constant-depth  jets.  However,  Prych  and  Stolzenbach  and  Harleman 
predict  a  buoyant  jet  that  spreads  into  a  thin  layer  (typically  1.0  to  1.5  m  thick), 
after  brief  contact  with  the  lake  bottom  for  a  distance  of  about  100  m.  The 
Stolzenbach-Harleman  model  predicts  several  small  ridges  in  the  jet-depth  con¬ 
tour  that  do  not  appear  reasonable.  The  sharp  decrease  in  jet-depth  within  the  first 
75  m  might  be  associated  with  the  increase  in  velocity  above  the  channel  velocity 
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near  the  outfall  predicted  by  Stolzenbach  and  Harleman.  The  Prych  model 
contains  no  mechanism  to  allow  for  a  difference  in  the  vertical  temperature 
and  velocity  spreading;  Stolzenbach  and  Harleman  employ  different  velocity 
and  temperature  similarity  functions.  Note  also  that  Prych  consistently 
predicts  thinner  plumes  for  the  four  Point  Beach  cases.  Bottom  interaction, 
as  predicted  by  the  above  three-dimensional  models,  appears  less  than  was 
actually  determined  by  field  measurements.  Tables  V-IX  summarize  the 
vertical  distributions  of  excess  temperature  and  excess  velocity  on  a  location- 
by-location  basis  for  the  four  dates  at  Point  Beach  and  one  at  Palisades.  The 
jet  appears  to  be  influenced  by  bottom  interaction  in  the  vicinity  of  the  outfall, 
followed  by  stratification  at  the  final  transect.  Changes  in  temperature  and 
velocity  excess  with  depth  tend  to  indicate  a  lens -shaped  plume,  with  a  greater 
vertical  spread  of  heat  than  momentum.  Estimates  of  the  plume  velocity  and 
temperature  half-depths  as  the  plume  stratifies  are  given  in  the  tables.  At 
Point  Beach,  the  plume  tends  to  stratify  to  a  temperature  depth  of  Z  m  from 
an  initial  temperature  depth  of  4.Z  m;  the  velocity  spreading  is  to  approxi¬ 
mately  1.5  m.  The  Palisades  plume  is  thinner:  It  reaches  an  asymptotic 
temperature  and  velocity  depth  of  about  0.5- 1.5  m  from  an  initial  discharge 
depth  of  Z.  1  m.  The  data  are  not  sufficiently  refined  to  permit  more  than 
just  a  gross  estimate.  Irregularities  in  the  decay  of  temperature  and  velocity 
excess  from  an  expected  monotonic  profile  are  due  perhaps  to: 

1.  Different  ambient  values  of  temperature  and  velocity  at  each 
depth. 

Z.  The  transient  thermal-front  phenomenon,  which  complicates 
the  concept  of  a  steady- state  temperature  at  a  single  location. 

This  irregular  behavior  is  real  and  cannot  be  ignored;  it  reflects  on  the 
transient  nature  of  the  jet  discharge  as  well  as  a  discharge  into  a  stratified 
ambient  environment. 

As  is  evident  from  the  data,  bottom  interaction  should  imply  greater 
lateral  spread,  since  vertical  entrainment  is  reduced;  even  the  observed 
spread  was  not  nearly  as  great  as  Prych  and  Stolzenbach-Harleman  predict 
for  no  interaction.  The  bottom  influence  appears  to  be  most  significant  for 
the  first  150  m  offshore  for  the  four  Point  Beach  cases  with  an  asymptotic 
approach  to  a  half-depth  of  temperature  and  velocity  of  1.0-1. 5  m. 

Figure  77  compares  model  predictions  for  plume  half-depth  along 
the  jet  centerline  with  actual  lake-bottom  depths  to  assist  in  determining 
bottom  interaction  for  October  10  at  Palisades.  Again,  each  model  predic¬ 
tion  in  Fig.  77  represents  the  half-depth  along  the  centerline  predicted  by 
that  model.  Predicted  model  centerlines  do  not  coincide,  but  are  sufficiently 
close  for  the  present  discussion.  The  lake-bottom  depth  shown  in  Fig.  77 
is  along  a  line  straight  out  from  the  outfall. 

As  can  be  seen  in  the  figure,  the  lake  depth  is  about  Z  m  to  a  distance 
of  1Z5  m  from  the  outfall;  then  the  bottom  drops  away  rather  sharply.  The 


TABLE  V.  DETERMINATION  OF  TEMPERATURE  AND  VELOCITY  HALF-DEPTHS 
AT  EACH  STATION  LOCATION  FROM  THE  POINT  BEACH  DATA  OF  MAY  18 ,  1972 
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TABLE  VI.  DETERMINATION  OF  TWERATURE  AND  VELOCITY  HALF-DEPTHS 
AT  EACH  STATION  LOCATION  FRCM  THE  POINT  BEACH  DATA  OF  MAY  23,  1972 
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TABLE  VII.  DETERMINATION  OF  TEMPERATURE  AND  VELOCITY  HALF-DEPTHS 
AT  EACH  STATION  LOCATION  FROM  THE  POINT  BEACH  DATA  OF  JULY  13,  1972 
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TABLE  VIII.  DETERMINATION  OF  TEMPERATURE  AND  VELOCITY  HALF-DEPTHS 
AT  EACH  STATION  LOCATION  FROM  THE  POINT  BEACH  DATA  OF  SEPTEMBER  9,  1972 
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TABLE  IX.  DETERMINATION  OF  TEMPERATURE  AND  VELOCITY  HALF-DEPTHS 
AT  EACH  STATION  LOCATION  FROM  THE  PALISADES  DATA  OF  OCTOBER  10,  1972 
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POINTS  OF  IRREGULAR  VERTICAL  EXCESS-VELOCITY  VARIATION:  18-25. 

FROM  POINTS  7-17,  THE  VELOCITY  HALF-DEPTH  OF  THE  STRATIFIED  JET  MAY  BE  ESTIMATED  AS  0.5-1. 5  METERS. 
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Fig.  77.  Half-depths  of  Temperature  and  Velocity  Distributions  Result¬ 
ing  from  Model  Calculations  for  Palisades:  October  10,  1972 

three-dimensional  Prych  and  Stolzenbach-Harleman  models  predict  a  buoyant 
surface  jet,  with  only  the  Stolzenbach-Harleman  temperature  half-depth  reach¬ 
ing  the  bottom.  On  a  point-by-point  basis  (see  Figs.  43-48),  the  temperatures 
at  Stations  9,  11.  14,  and  1  5  do  not  decline  to  half  their  near-surface  values 
before  the  lake  bottom  is  reached.  All  other  points  indicate  a  1.0-  to  1.5-m 
temperature  half-depth.  It  is  inferred  from  the  data  that  the  jet  hugs  the 
bottom  for  the  first  125  m  and  then  rises  to  a  1.0-  to  1.5-m  half-depth  as  the 
bottom  drops  away  and  jet  buoyancy  becomes  predominant.  A  point-by-point 
analysis  of  the  velocity  data  indicates  no  real  trend.  The  bottom  was  usually 
encountered  on  the  first  two  transects  before  the  half-depth  was  reached.  The 
velocity  half-depth  appears  to  be  1.0- 1.5  m  otherwise.  The  variation  of  ve¬ 
locity  data  with  depth  at  Palisades  is  just  too  erratic  for  any  consistent  con¬ 
clusions.  This  also  was  true  for  the  Palisades  data  of  June  14  and 
July  19. 

E.  Isotherm  Areas 


Isotherm  areas  were  calculated  based  upon  the  temperature -distribution 
function  resulting  from  the  fit  to  the  data  by  JETFIT.  The  data  were  not  suf¬ 
ficiently  detailed  to  contour  isotherms  and  determine  areas  directly.  Conse¬ 
quently,  the  JETFIT  area  calculations  should  be  used  with  caution.  It  is 
encouraging,  however,  to  note  that  isotherm  areas  for  the  four  dates  at 
Point  Beach  are  approximately  the  same;  this  is  to  be  expected,  since  plant 
conditions  were  nearly  the  same  and  ambient  currents  were  weak,  if  present 
at  all.  These  isotherm  areas,  along  with  the  model  predictions,  are  presented 
in  Figs.  78-82. 
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Fig.  82.  Isotherm  Areas  Resulting  from  the  Fit¬ 
ting  Procedure  and  Model  Calculations 
for  Palisades:  October  10,  1972.  ANL 
Neg.  No.  190-937. 


The  Pritchard  model  appears 
to  be  most  accurate  in  area  predic¬ 
tions  for  the  four  Point  Beach  dates. 
Any  discrepancies  appear  to  be  on 
the  conservative  side.  For  the  no¬ 
current  cases  (Figs.  78,  79,  and  82), 
the  areas  predicted  by  Stolzenbach 
and  Harleman  are  reasonable  but 
optimistic  for  the  region  in  which 
data  were  taken  (to  3-4C°  isotherm). 
Beyond  the  data,  however,  the  model 
predicts  sharply  increasing  areas, 
perhaps  due  to  the  large  lateral 
spreading  predicted.  The  current 
cases  of  July  13  and  September  9 
reveal  a  less  steep  increase  in  areas 
with  decreasing  temperature  excess; 
the  model  for  those  dates  tends  to 
underpredict  areas  for  the  region  in 
which  data  exist.  Too  sharp  a  dis¬ 
tinction  between  the  no -current  and 
weak-current  cases  appears  in  the 
area  predictions.  The  Stolzenbach- 
Harleman  model  as  derived  is  too 
complex  to  allow  any  speculation  as 
to  possible  reasons  for  that  behavior. 


The  Prych  model  gives  rea¬ 
sonable  areas,  but  they  are  consistently  lower  than  the  data.  Also,  the  pre¬ 
dictions  are  nearly  the  same  for  all  dates.  This  is  to  be  expected,  since  the 
plant  and  environmental  conditions  varied  only  slightly  from  date  to  date  at 
Point  Beach.  Motz- Benedict  predictions  are  consistently  lower  on  all  dates. 
A  different  choice  of  parameters  in  this  regard  could  affect  these  results 
greatly. 


The  Palisades  area  calculations  are  difficult  to  evaluate.  Most  of  the 
jet  data  taken  appears  to  be  in  the  region  of  flow  establishment.  Although  the 
actual  excess  temperature  was  9-2C°,  the  data  fit  indicated  a  temperature 
excess  of  7.8C°  for  a  good  200  m  offshore.  The  models  all  underpredict  areas 
at  the  7C°  isotherm  shown  by  the  dot  on  Fig.  82;  little  can  be  said  beyond  that 
isotherm,  since  no  data  exist.  The  Pritchard  curve  has  a  sharp  increase  in 
area  at  approximately  the  6C°  isotherm  due  to  Pritchard's  method  of  halting 
vertical  entrainment  when  a  predicted  plume  depth  of  3.05  m  (10  ft)  is 
reached. 

F.  Decay  of  Centerline  Temperature  and  Velocity  with  Depth 


Observation  of  the  Point  Beach  field  measurements  with  depth  reveals 
an  approximate  constancy  of  temperature  and  velocity  for  the  depths  measured 
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(to  a  depth  of  Z.  5-3.0  m).  It  is  expected  that,  beyond  the  2. 5-3.0  m  of  uniform 
profile,  there  exists  a  vertical  region  for  which  temperature  and  velocity  ex¬ 
cesses  decrease  to  zero  (see  Fig.  83). 
These  vertical  distributions  are  per¬ 
haps  reminiscent  of  fully  developed 
turbulent  flow  in  pipes,  which  consists 
of  a  large  uniform  temperature  and 
velocity  core,  combined  with  a  sharp 
approach  to  the  boundary  value  (wall 
of  pipe,  or  here,  the  lake  bottom). 

The  results  of  the  data-fitting 
procedure  for  Point  Beach  appear  in 
Figs.  84-87.  The  points  in  these 
four  figures  are  the  centerline  tem¬ 
perature  and  velocity  excesses  from 
the  fitting  procedure  at  150  m  from  the  outfall  along  the  centerline.  Model 
predictions  of  temperatures  and  velocities  with  depth  at  150  m  are  also 
plotted.  The  lake  depth  at  this  distance  is  about  3. 0-3. 5  m. 

The  chosen  distance  of  1  50  m  is  about  two  to  three  times  the  length 
of  the  region  of  flow  establishment;  only  the  two-dimensional  models  predict 
uniformity  of  vertical  profiles  at  that  distance.  Although  predictions  of  the 
Prych  and  Stolzenbach-Harleman  models  look  poor,  their  predictions  for  the 
centerline  values  are  too  low,  biasing  the  entire  vertical-distribution  pre¬ 
dictions.  More  accurate  centerline  predictions  would  perhaps  have  yielded 
better  agreement  with  the  data.  The  three-dimensional  models  evaluated  here 
predict  too  rapid  a  temperature  and  velocity  decay  with  depth.  The  flow- 
establishment  assumptions  for  vertical  temperature  decay  of  Stolzenbach  and 
Harleman  (constant-temperature  vertical  core,  with  small,  turbulent  boundary 
layer  below  it)  also  appear  to  be  appropriate  in  the  region  of  established  flow. 
The  Stolzenbach-Harleman  decay  with  depth  beyond  the  region  of  flow  estab¬ 
lishment  is  essentially  parabolic,  with  no  jet  effects  beyond  a  certain  depth. 

The  Prych  decay  is  Gaussian  in  the  vertical  direction.  The  Pritchard  and 
Motz- Benedict  models,  being  two-dimensional,  predict  no  change  in  behavior 
with  depth. 

The  Palisades  data  are  more  of  an  enigma  (see  Fig.  88).  The  center- 
line  distance  of  1  50  m  is  expected  to  border  the  end  of  the  region  of  flow 
establishment,  yet  the  vertical  temperature  and  velocity  decay  appears  quite 
sharp.  The  depth  of  constant  temperature,  if  any,  is  quite  small.  Complicating 
factors  here  are  the  irregular  shallow  bottom,  approximately  2  m  in  depth,  a 
diverging  surface  outfall  structure,  and  a  nonuniform  velocity  distribution. 

More  detailed  data  would  be  required  before  further  conclusions  could  be 
drawn. 


TEMPERATURE  WITH  DEPTH 
150  METERS  ALONG 
CENTERLINE 


VELOCITY  EXCESS  WITH  DEPTH 
150  METERS  ALONG  CENTERLINE 


Fig.  83.  Idealized  Decay  of  Temperature  Ex¬ 
cess  and  Velocity  Excess  with  Depth 


117 


<oD)  °Q  (39S/UU3)  3n 


2  oo 

n)  c 
co 

CO  r-H 

(U  9 
o  £ 
x  <U 
w  cd 


u 

o 

> 

X> 

a 

03 


D 

o 

E 

o 


<L)  ^ 

^  O 
r*  ^ 

g  - 

3  03 

«  .g 
<3  cx 
a.  ^ 

GQ 

<L)  '-+-< 

H  O 

<L)  C 
C  O 

•  f-H  «i-H 
i— H  4—* 

«-*  o 
<u  a 

a  rD 

(U  tL< 
U  03 


io 

00 


oo 

•i-H 

Lu 


<v 

E 


on 

03 

co 

CO 

<D 

U 

X 

W 

4-J 

o 

o 

*3 

> 

X> 

c 

03 


3 

u 

X 

w 

<u 

»-< 

D 

*-j 

03 

M 

0) 

a 

E 

<u 

H 

d> 

a 

a 


c 

a> 

O 


oo 

c 


3 

CO 

<D 

0^ 


=J 

O 


£ 

o 


E 

o 

uo 


03 

*-» 

Oh 

<D 

Q 


a 

o 


o 

a 

3 

Lu 


00 

00 

•H 

U-> 


(.D)  °e 


(39S/UJ3) 


from  the  Fitting  Procedure  and  Model  Calculations  for  from  the  Fitting  Procedure  and  Model  Calculations  for 

Point  Beach:  May  18,  1972  Point  Beach:  May  23,  1972 


118 


(oD)  °Q  (33S/W0)  3n 


oo 
c  , 


<u  3 

u  £ 

X  <u 

w  cd 


cj 

u 

> 

TP 

a 

03 


3 

o 

03 

u 

*3 

TP 

o 


<D 
CJ 
*  10 
W  rH 

g  - 

3  03 

3  -S 

Cl, 

E  Q 

<U  L-. 

H  ° 


TP 

C 

03 


<D  05 

9  m 
-5  V 
V  c 
y  E 
2  g 

Cl  Oh 
00  ,«> 
c  1/5 


u,  u 


0)  c 

C  r3 

a)  Uh 

U  03 


03 

(D 

CQ 


E  c 


o 

CL 


00 


oo 

•H 

Uh 


<v 

e 


03  00  , 

C 

00  .d 

to  n 

a) 
o 
X 

ui 


CJ 

o 

<D 

> 

TP 

C 

03 


<U 

cj 

X 


a 

o 

*d 

03 
*— 1 
3 
O 

03 

u 

*3 

TP 

o 


TP 

G 

03 


W  ^ 
0)  T“< 
§  3 
2  <5 

CL 
CL  ni 

E  Q 

« 

H  o 

<u  c 
c  o 

•  !—*  "H 
f—i  *-> 

5  c 

C  3 

(U  b-, 

U  «> 


<M 
C- 

«— t  OP 

TP  rH 
<D 

CJ  CO 
O  rH 

t-4 

CL  >N 

00  3 

C  — > 


Uh  CJ 

«  S 

3  CQ 
4-J  ^ 

E  s 

2  O 

id  CL 


CO 

CO 


00 

•  H 

Uh 


UP) 


(D3S/UU3)  3fl 


119 


Fig.  88.  Centerline  Temperature  Excess  and  Velocity  Excess  as 
a  Function  of  Depth  at  150  m  from  Outfall  Resulting 
from  the  Fitting  Procedure  and  Model  Calculations  for 
Palisades:  October  10,  1972 

G.  Variation  of  Temperature  and  Velocity  Width  with  Depth 


Each  mathematical  model  predicts  a  constant  velocity  and  temperature 
width  with  depth.  The  two-dimensional  models  do  so  by  the  nature  of  their  two- 
dimensionality,  the  three-dimensional  models  by  definition  of  their  lateral  and 
vertical  profiles.  The  Stolzenbach-Harleman  and  Prych  models  assume,  for 
any  cross  section  of  the  jet  normal  to  the  centerline,  the  same  decay  laterally 
independent  of  depth,  as  well  as  identical  decay  vertically,  irrespective  of 
lateral  distance.  Consequently,  the  widths  must  maintain  a  constant  value 
with  depth. 

The  data,  however,  show  a  more  lens-shaped  profile,  indicating  that  the 
temperature  and  velocity  widths  decrease  with  depth  (see  Figs.  89-92).  The 
local  centerlines  (with  depth)  do  not  coincide,  yet  are  sufficiently  close  to  indi¬ 
cate  that  such  vertical  profile  results  may  be  meaningful.  From  the  figures, 
the  temperature  width  appears  to  be  uniform  with  depth  for  about  1.5  m  and 
then  decreases  abruptly;  the  velocity  width  appears  constant  for  a  shorter 
depth  distance  (0.5-1  m)  before  decreasing  toward  zero. 
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Fig.  89.  Half-widths  of  Temperature  and  Velocity  Distributions  as  a  Function  Fig.  90.  Half-widths  of  Temperature  and  Velocity  Distributions  as  a  Function 
of  Depth  at  150  m  from  Outfall  Resulting  from  the  Fitting  Procedure  of  Depth  at  150  m  from  Outfall  Resulting  from  the  Fitting  Procedure 

and  Model  Calculations  for  Point  Beach:  May  18,  1972  and  Model  Calculations  for  Point  Beach:  May  23,  1972 
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Fig.  91.  Half-widths  of  Temperature  and  Velocity  Distributions  as  a  Function  Fig.  92.  Half-widths  of  Temperature  and  Velocity  Distributions  as  a  Function 
of  Depth  at  150  m  from  Outfall  Resulting  from  the  Fitting  Procedure  of  Depth  at  150  m  from  Outfall  Resulting  from  the  Fitting  Procedure 

and  Model  Calculations  for  Point  Beach:  July  13,  1972  and  Model  Calculations  for  Point  Beach:  September  9,  1972 


The  Palisades  data  (see  Fig.  93)  are  insufficient  and  irregular,  due 
again  perhaps  to  the  shallow  lake  depths  in  the  discharge  vicinity  as  well  a 
the  nonuniformity  of  the  diverging  discharge. 
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Fig.  93.  Half-widths  of  Temperature  and  Velocity  Distributions  as  a  Function 
of  Depth  at  150  m  from  Outfall  Resulting  from  the  Fitting  Procedure 
and  Model  Calculations  for  Palisades:  October  10,  1972 
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VIII.  SUMMARY  AND  CONCLUSIONS 

A.  Field-data  Acquisition 

The  experimental  technique  described  in  Sec.  II  for  making  meas¬ 
urements  in  the  jet  regime  represents  a  compromise  between  speed  and  com¬ 
pleteness.  More  data  in  a  shorter  period  of  time  are  needed  for  a  more 
accurate  comparison  with  model  predictions.  Given  the  manpower  and  equip¬ 
ment  limitations  under  which  the  work  was  done,  the  technique  reported  on 
in  this  report  gives  reasonable  results  in  terms  of  number  of  data  points 
versus  time  of  measurement.  The  data  obtained  in  this  manner  reflect  the 
gross  features  of  the  jet  regime,  but  have  two  limitations: 

1.  Since  the  temperature  and  velocity  at  a  point  in  the  jet  are  not 
constant,  a  single  measured  value  of  each  may  not  accurately  define  these 
parameters  at  a  point. 

2.  The  measurements  are  not  made  simultaneously.  Typically,  they 
are  collected  over  a  period  of  3-7  hr  during  which  time  factors  affecting  the 
jet  may  change. 

B.  Data- smoothing  Technique 

The  fitting  procedure  described  in  Sec.  V  was  used  to  extract  the 
gross  features  of  the  temperature  and  velocity  distributions  from  the  field 
measurements.  Some  procedure  of  this  type  is  necessary  because  the  data 
consist  of  measurements  at  isolated  points  in  the  temperature  and  velocity 
fields.  The  functional  forms  used  in  this  study  were  chosen  for  several  rea¬ 
sons.  One  reason  is  their  similarity  to  the  forms  resulting  from  the  simple 
yet  fairly  widely  used  jet  analyses  of  Pritchard17  and  Motz  and  Benedict.5 
Another  reason  is  that  preliminary  examination  of  the  Point  Beach  data  in¬ 
dicated  that  such  forms  might  represent  the  available  data  fairly  well.  Finally, 
it  was  thought  that  for  this  initial  attempt  at  data  fitting,  simple  forms  with 
limited  numbers  of  free  parameters  would  be  appropriate.  Clearly,  unless 
the  actual  functional  forms  are  known,  any  fitting  procedure  of  this  type  will 
have  limitations  and  biases. 

The  results  of  the  present  fitting  procedure  can  be  considered  to  be  a 
successful  first  attempt.  It  is  difficult  to  assess  the  goodness  of  fit  in  a  mean¬ 
ingful,  quantitative  way;  therefore  a  point- by- point  tabulation  has  been  included 
in  Appendix  D.  The  average  deviations  of  the  fitted  function  from  the  data  for 
any  particular  set  of  measurements  never  exceeded  1  . 1  C°  for  temperature  and 
8  cm/sec  for  velocity.  These  values  correspond  to  12  and  14%  of  the  initial 
excess  temperature  and  average  outfall  velocity,  respectively.  Certainly  at 
least  the  gross  behavior  of  the  near-field  temperature  and  velocity  distribu¬ 
tions  has  been  accounted  for  with  the  present  fitting  procedure. 


124 


C.  Analytical  Model;  Field-data  Comparisons 


The  results  of  the  model-data  comparisons  are  summarized  here  for 
each  model  evaluated. 

1.  Pritchard  Model 


The  Pritchard  model  compared  best  with  the  Point  Beach  Unit  1 
and  Palisades  jet  data  presented  in  this  report.  Model-data  discrepancies 
usually  indicated  conservative  predictions  on  the  part  of  the  model.  The  suc¬ 
cess  of  Pritchard's  model  for  these  cases  is  due  in  part  to  the  persistent 
efforts  of  the  author  in  calibrating  the  model  with  field  and  hydraulic  data. 
Also,  the  field  data  used  by  Pritchard  in  the  development  of  his  model  were 
generally  from  plants  with  similar  discharge  and  environmental  character¬ 
istics  to  Point  Beach  and  Palisades. 

This  basically  empirical  or  phenomenological  model  is  not  ex¬ 
pected  to  perform  well  when  buoyancy  is  significant  (initial  densimetric 
Froude  number  less  than  about  2)  or  when  the  dilution  capacity  of  the  re¬ 
ceiving  water  is  significantly  restricted  by  lateral  or  bottom  boundaries. 

The  model  is  also  weak  in  its  treatment  of  vertical  mixing;  the  model  has 
scarcely  been  tested  for  prototype  situations  in  which  vertical  entrainment 
is  significant. 

Although  the  model  is  too  simplistic  to  handle  the  major  compli¬ 
cations  that  may  be  important  in  certain  field  cases  (buoyancy,  cross  currents, 
etc.),  it  might  well  give  good  results  for  cases  of  the  highe r- velocity  dis¬ 
charges  (i.e.,  larger  densimetric  Froude  number  problems)  for  the  critical 
case  of  a  stagnant  lake.  Further  verification  work  for  this  model  is  desirable 
at  different  sites  describing  a  variety  of  geometric,  kinematic,  and  dynamic 

conditions . 

« 

2.  Motz- Benedict  Model 


The  Motz- Benedict  recommendations  for  the  choice  of  entrain¬ 
ment  coefficient  E  apparently  vary  over  too  large  a  range  for  the  data  ob¬ 
served  at  Point  Beach  and  Palisades.  The  recommended  value  of  0.05  for 
the  stagnant- lake  case  gave  a  reasonably  good  approximation  for  centerline 
temperature  and  velocity  decay.  Values  of  E  for  A  (ratio  of  ambient  velocity 
to  initial  jet  velocity)  between  0  and  0.2  were  determined  by  a  linear  interpo¬ 
lation  from  the  known  values  of  0.05  (A  =  0)  and  0.4  (A  =  0.2).  This  appeared 
to  us  to  be  as  good  an  approximation  for  E  as  any  in  the  presence  of  the  wide 
scatter  in  the  Motz- Benedict  data  from  which  E  is  to  be  determined.  The 
results  indicate  that  smaller  values  of  E  for  the  cases  of  nonzero  ambient 
crosscurrents  would  have  been  more  appropriate.  Temperature  and  velocity 
decays  would  not  have  been  as  rapid  and  plume  widths  would  have  been  larger. 
The  most  sensitive  parameters  are  E  and  r  =  b0,/bo.  Data  for  the  determination 
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of  E  have  much  scatter;  data  for  r  have  not  been  determined.  The  inade¬ 
quacies  of  this  simple,  two-dimensional  model  are  briefly  summarized  as 
follows: 


a.  Buoyant  forces,  as  well  as  vertical  entrainment,  are  ignored 
in  the  model  development.  For  densimetric  F roude  numbers  less  than  about  2, 
buoyant  spreading  is  quite  important;  for  densimetric  Froude  numbers  ex¬ 
ceeding  about  5,  a  substantial  amount  of  vertical  entrainment  may  be  expected. 
Consequently,  the  types  of  discharges  for  which  the  model  may  be  used  is 
limited.  The  entrainment  coefficient,  determined  from  data  fitting,  must 
actually  account  for  spreading  due  to  buoyancy  as  well  as  from  jet  entrainment. 

b.  The  choice  of  a  coefficient  of  entrainment  must  be  based  on 
inconsistent  data  as  expressed  above.  Data  presented  by  Motz  and  Benedict 
show  considerable  variation  in  the  value  of  entrainment  coefficient  from  situ¬ 
ation  to  situation.  Particularly  for  lake  data,  the  value  of  E  is  shown  to  be 
extremely  small,  -0.04,  in  relation  to  riverine  situations  where  E  is  about 
an  order  of  magnitude  larger.  A  particularly  interesting  contrast  is  afforded 
when  one  compares5  the  Motz- Benedict  laboratory  data  for  a  90°  discharge  to 
the  Rombe rg- Ayres  lake  data.  Although  most  of  the  initial  conditions  are 
similar,  the  entrainment  coefficients  are  different  by  about  a  factor  of  10. 

Part  of  the  problem  probably  lies  with  the  fact  that  Motz  and  Benedict  ana¬ 
lyzed  lake  data  from  locations  in  which  prominences  (breakwaters,  for  ex¬ 
ample)  were  present.  A  second  difficulty  is  that  the  limited  field  results 
indicated  that  the  entrainment  coefficient  remained  essentially  constant  at 

a  particular  location  with  changing  values  of  A;  this  behavior  is  consistent 
with  the  laboratory  findings  showing  E  to  be  relatively  independent  of  A  for 
each  site.  Here,  turbulent  intensity  of  the  ambient  current  (see  item  e  below) 
may  be  a  significant  factor  that  was  ignored. 

c.  We  believe  the  model  is  actually  restricted  to  small  ambient 
currents.  The  simulation  assumes  that  the  jet  velocity  approaches  zero  at 
large  distances  normal  to  the  jet  centerline.  The  assumption  of  Gaussian 
profiles  for  temperature  and  velocity  and  the  assumption  of  equal  rates  of 
entrainment  on  offshore  and  lee  sides  of  the  jet  are  not  valid  for  ambient 
currents  that  are  not  very  small. 

d.  The  model  does  not  simulate  unequal  rates  of  spread  for 
momentum  and  heat,  as  has  been  observed  to  be  significant  in  the  Point  Beach 
and  Palisades  data.  The  Prych  and  Pritchard  models  also  do  not  distinguish 
between  these  rates  of  spread. 

e.  The  model  does  not  treat  turbulent  intensity  in  the  crossflow, 
due  to  the  assumption  that  such  turbulence  and  its  effect  on  mixing  are  neg¬ 
ligible.  Although  some  account  of  this  effect  may  be  made  in  the  choice  of 
entrainment  coefficient,  any  such  treatment  would  be  only  qualitative  in  nature. 
The  other  models  do  not  treat  this  phenomenon  adequately  either;  the  Prych 
model  does  include  ambient  turbulence  in  terms  of  a  horizontal  and  vertical 
ambient  eddy-thermal  diffusivity. 


126 


Deficiencies  a  and  b  could  eventually  be  fatal  to  the  model.  Im¬ 
proved  predictions  might  be  obtained  by  determining  the  free  parameters 

E  =  entrainment  coefficient, 
r  =  b0/bo 

and 


s  =  ratio  of  lateral  spread  of  heat  to  that  of  momentum 

for  each  set  of  data  available,  with  the  hope  that  a  consistent  trend  (or  corre¬ 
lation)  might  develop  as  the  initial  densimetric  Froude  number,  initial  angle 
of  jet  discharge  with  the  current,  velocity  ratio  A,  and  w  (ambient  stream 
width/discharge  width)  vary. 

3.  Stolzenbach-Harleman  Model 


The  Stolzenbach-Harleman  model  generally  compared  poorly  with 
the  jet  data  of  Point  Beach  and  Palisades.  Centerline  temperature  and  velocity 
decay  were  predicted  to  be  too  rapid;  the  lateral  spread  of  the  jet  was  much 
too  great.  The  model  does  not  consider  jet  interaction  with  the  lake  bottom, 
which  does  occur  to  some  degree  near  the  outfall.  Such  interaction  should 
provide  some  increase  in  lateral  spread  due  to  restricted  dilution  at  the  jet- 
lake  bottom  interface.  The  predicted  lateral  spread,  however,  greatly  ex¬ 
ceeded  the  observed  spread,  even  with  bottom  effects  assisting  that  lateral 
growth.  Surprisingly,  the  model  tends  to  underpredict  lateral  spread  in  the 
region  of  flow  establishment.  These  poor  comparisons  of  model  to  data  may 
be  traced,  in  part,  to  the  model  assumption  on  late ral- spreading  velocity, 
which  was  based  more  on  physical  intuition  than  on  any  data. 

A  second  major  fault  of  the  model  lies  in  the  presumed  jet  struc¬ 
ture  based  upon  nonbuoyant  jet  theory.  First,  the  four- zone,  rectangular  jet 
structure  may  not  be  valid  for  buoyant  jets.  In  particular,  the  cross  section 
of  a  buoyant  jet  is  normally  taken  to  be  lens- shaped  rather  than  rectangular, 
as  supposed  by  the  model.  The  assumed  division  of  the  jet  into  four  distinct 
regions  necessitates  that  interregional  velocities  be  specified.  The  forms  of 
these  velocities  are  unknown  and  therefore  must  be  guessed.  Stolzenbach  and 
Harleman  also  require  that  no  turbulent  momentum  transfer  occurs  between 
regions  of  the  jet  or  between  the  jet  and  ambient  water.  This  is  tantamount 
to  dropping  the  Reynolds  stress  terms  in  the  equations  of  motion  or,  equiva¬ 
lently,  dropping  the  turbulent- diffusion  mechanism.  As  a  consequence,  turbu¬ 
lent  jet  diffusion  had  to  be  artificially  simulated  through  the  entrainment 
coefficient  and  similarity  forms  for  temperature  and  velocity.  In  any  case, 
some  calibration  of  the  model  to  actual  field  and  hydraulic  data  might  have 
provided  better  predictions. 

Aside  from  the  largely  theoretical  criticisms  of  the  model,  there 
are  practical  difficulties  in  actually  obtaining  a  numerical  solution  to  the  set 
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of  ordinary  differential  equations.  Due  to  the  complexity  of  the  set,  derivatives 
must  be  found  by  solving  a  linear  set  of  algebraic  equations  before  applying  a 
Runge-Kutta  scheme.  For  many  cases,  the  matrix,  which  must  be  reduced,  is 
nearly  singular  and  much  precision  is  lost  in  solving  for  the  derivatives. 

Among  the  problems  that  may  be  encountered  in  using  the  code  are: 

a.  Width  predictions  may  decrease  by  as  much  as  50%  with  each 
o  rde  r- of- magnitude  reduction  in  the  error  criteria  until  the  program  finally 
fails'. 


b.  The  program  may  not  run  for  cases  of 

(1)  Low  initial  densimetric  Froude  number. 

(2)  Low  aspect  ratio. 

(3)  Initial  angles  greater  than  90°. 

c.  Numerical  underflows  must  be  suppressed  for  successful 
completion  of  many  runs. 

d.  Differences  in  machine  precision  due  either  to  differences  in 
word  structure  or  to  differences  in  the  operating  system  may  cause  differences 
of  up  to  5%  . 

4.  Prych  Model 


The  Prych  predictions  have  the  same  problems  as  those  of 
Stolzenbach  and  Harleman: 

a.  Too  rapid  a  decay  in  centerline  temperature  and  velocity. 

b.  Too  great  a  lateral  spread. 

The  Prych  model  also  compares  poorly  with  the  Point  Beach  and  Palisades  data. 

We  suspect  that  a  major  difficulty  with  the  theoretical  development 
is  in  the  assumption  for  a  late  ral- spreading  velocity  based  upon  the  analogy  to 
the  celerity  of  a  density  front  of  a  uniform  depth  with  a  uniform  density  differ¬ 
ence.  This  model  for  lateral  spreading  is  apparently  incorrect  as  simulated. 
Also,  the  hydrostatic  pressure  force  is  simulated  to  act  in  the  longitudinal  di¬ 
rection  only.  Pressure  forces  in  reality  act  longitudinally  and  laterally,  with 
an  approximate  hydrostatic  distribution  vertically.  The  as  sumption  of  a  fictitious 
late  ral- s preading  velocity  was  made  to  remedy  that  omission.  The  adequacy  of 
the  Prych  simulation  of  ambient  turbulence  and  shear  stresses  has  not  been 
verified.  Calibration  of  the  model  and  its  empirical  coefficients  with  hydraulic 
or  prototype  data  might  have  improved  predictions. 

As  with  the  Stolzenbach- Ha rleman  model,  the  Prych  model  is  ap¬ 
plicable  for  small  or  zero  ambient  currents  due  to  the  assumption  of  similarity 
for  temperature  and  velocity  profiles,  as  well  as  equal  entrainment  on  offshore 
and  lee  sides  of  the  jet.  The  computer  code  developed  by  Prych  operates  well 
(model  equations  are  simpler  than  those  of  Stolzenbach  and  Harleman);  this 
makes  it  easier  for  future  alteration,  manipulation,  and  calib ration  of  the  model. 
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IX.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

A.  Field-data  Acquisition 

Other  techniques,  such  as  fixed  instrument  arrays,  aerial  photography 
and  aerial  infrared  imagery,  and  nonstationary  measurement  systems,  may 
prove  worthwhile  in  future  measurements  in  the  jet  regime.  These  techniques 
or  combinations  of  techniques  may  make  it  possible  to  overcome  the  limita¬ 
tions  inherent  in  the  fixed-boat  method  described  here,  but  will  probably  also 
result  in  increased  complexity  of  the  measurement  and  significant  monetary 
commitment. 

B.  Data- smoothing  Technique 

To  extract  as  much  information  as  possible  from  the  data,  the  most 
general  functional  forms  practicable  should  be  used.  This  would  require  more 
free  parameters  and  a  more  involved  and  lengthy  fitting  process.  Future  at¬ 
tempts  to  extend  this  method  of  data  analysis  might  include  some  of  the  following: 

1.  Instead  of  Gaussian  lateral  profiles,  a  function  that  allows  for  a 
flat  region  near  the  jet  centerline  might  be  chosen.  This  could  then  simulate 
the  core  region  included  in  the  Stolzenbach- Harleman  analysis. 

2.  The  form  of  the  centerline  temperature  excess  should  be  such  as 
to  require  that  it  extrapolate  to  the  measured  excess  at  the  outfall  (s  =  0). 

The  present  form  does  not  have  this  property. 

3.  In  the  present  procedure,  the  rate  of  dropoff  of  the  centerline 
temperature  and  velocity  excesses  is  fixed  as  being  inversely  proportional 

to  the  one-half  power  of  s,  the  distance  from  the  outfall.  Instead  of  this  being 
restricted  to  the  one-half  power,  it  could  be  left  as  a  free  parameter  to  be  de¬ 
termined  by  fitting  to  the  data.  The  more  recent  phenomenological  model  by 
Pritchard24  has  employed  this  form  for  centerline  decay  of  temperature. 

4.  Additional  parameters  could  be  added  to  the  expressions  for  the 
widths  to  allow  them  to  vary  quadratically  with  s  instead  of  linearly.  The 
alternatives  are  limitless,  and  only  through  repeated  attempts  at  a  fitting  pro¬ 
cedure  can  it  be  determined  whether  significant  improvements  are  possible. 

C.  Analytical  Model;  Field-data  Comparisons 


1.  An  attempt  should  be  made  to  calibrate  the  Motz- Benedict, 
Stolzenbach- Harleman,  and  Prych  models  to  field  and  hydraulic  data.*  The 
Pritchard  model  (No.  l)  should  be  further  tested  with  field  data  from  other 
sites  as  well  as  available  data  from  physical  hydraulic  models.  The  new  nu¬ 
merical  models  of  Brady  and  Geyer,  Paul  and  Lick,  and  Waldrop  and  Farmer 

*'Work  is  presently  underway  by  Dr.  M.  Shirazi  at  the  Pacific  Northwest  Environmental  Research  Laboratory 
at  Corvallis  to  improve  the  Stolzenbach-Harleman  and  Prych  models  by  calibration  with  data.  At  this 
writing,  a  successful  modification  and  calibration  of  the  Prych  model  appears  imminent. 
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look  promising,  and  attempts  should  be  made  at  verification  with  prototype 
field  data.  The  new  phenomenological  model  of  Pritchard  (No.  2),  based  upon 
52  sets  of  model  and  prototype  data,  also  looks  promising  and  should  be 
verified. 

2.  Considerably  more  data  are  required  from  more  ideal  or  classical 
types  of  discharge  structures  for  verification.  The  Palisades  data  had  too 
many  irregularities  (rough  shallow  bottom  and  diverging  discharge  channel) 
for  adequate  model  evaluation  for  those  integral-type  models  studied  in  this 
report.  Data  are  required  for  model  verification  (for  both  integral  and  nu¬ 
merical  models),  which  include  a  wide  range  of  densimetric  Froude  numbers, 
aspect  ratios,  bottom  slopes,  angle  of  discharge,  ambient  currents,  etc.,  to 
provide  a  fair  and  wide  variety  of  test  situations  for  the  models.  Only  when 
this  large  body  of  data  (physical  model  or  preferably  prototype  field  data)  be¬ 
comes  available  will  it  be  possible  to  fully  and  fairly  evaluate  and  improve,  or 
develop,  better  models. 

From  our  verification  efforts  supplemented  by  the  work  done  by 
Dr.  M.  Shirazi  at  the  Pacific  Northwest  Water  Laboratory  of  the  EPA,  we 
recommend  that: 

1.  The  Stolzenbach-Harleman  and  Prych  models  not  be  used  as  they 
exist  in  their  present  form  for  those  cases  when  significant  bottom 
interaction  is  expected.  However,  the  vast  majority  of  prototype 
situations  do  have  some  bottom  interaction. 

2.  The  Motz  -  Benedict  model  be  used  for  stagnant  ambient  water  case 
only  (with  an  entrainment  coefficient  on  the  order  of  0.04) 

3.  The  Pritchard  model  be  used  for  stagnant  receiving  water  only. 

Further  analytical  work  is  necessary  to  modify  the  Prych  and  Stolzenbach- 
Harleman  models  so  they  can  be  used  for  shallow  water  surface  discharges. 
Additional  work  in  model  calibration  is  necessary  before  a  viable  form  of  the 
Motz  -  Benedict  model  can  be  achieved  for  ambient  currents.  More  verifica¬ 
tion  work  on  the  Motz  -  Benedict  (no  current)  and  Pritchard  models  would  be 
useful. 
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APPENDIX  B 

Preliminary  Feasibility  Study 


On  November  3,  1971,  a  preliminary  feasibility  study  of  the  technique 
for  studying  the  temperature  and  velocity  profiles  was  made  near  the  outfall 
of  the  Point  Beach  Nuclear  Power  Plant  (Unit  1).  A  three-point  mooring  sys¬ 
tem  was  used  to  hold  the  boat  steady  while  simultaneous  temperature  and 
velocity  measurements  were  obtained  in  the  near-field  region.  Anchors  were 
located  on  either  side  of  the  plume  and  attached  to  the  stern  cleats  of  the 
boat,  and  a  bowline  was  attached  to  the  center  of  the  outfall.  Transects  across 
the  plume  centerline  were  then  made  at  about  8,  27,  and  73  m  from  the  outfall 
by  pulling  the  boat  from  one  side  anchor  to  the  other.  The  position  of  the  boat 
was  held  relatively  constant,  and  transits  were  used  to  obtain  the  location  of 
each  station.  (Station  locations  for  this  jet  study  are  shown  in  Fig.  94.)  A 
Bendix  Q-15  current  meter  with  a  YSI  thermistor  attached  was  used  to  measure 
current  velocity  and  water  temperature.  The  meter  was  lowered  over  the  side 
and  suspended  at  2-ft  intervals  to  a  depth  of  10  ft  or  to  the  bottom.  The  lake 
depth  was  4.1  m  at  the  outfall  and  decreased  in  depth  to  2  1  m  at  a  point  73  m 
from  the  outfall. 


Fig.  94.  Station  Locations  for  Jet-regime  Study:  November  3,  1971, 

1245-1605  Hours.  ANL  Neg.  No.  190-573. 

The  data  with  the  range  of  variability  are  plotted  in  Figs.  95-97  for 
three  different  depths.  Drawings  for  the  8-  and  1 0  -  ft  depths  were  not  made 
because  data  were  not  available  at  all  stations.  The  temperature  at  a  given 
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depth  and  location  exhibited  greater  variability  than  expected,  even  near  the 
plume  centerline.  As  a  result,  temperature  ranges  are  presented  in  the 
figures  with  the  lowest  temperature  in  parentheses.  Each  measurement  was 
made  over  a  period  of  about  1  min,  which  is  indicative  of  the  rapid  variability 
of  the  temperature.  Variations  in  the  velocity  were  not  as  easily  observed, 
because  of  the  current-meter  time  constant.  The  lengths  of  the  arrows  in  the 
figures  are  proportional  to  current  magnitude,  and  their  directions  indicate 
current  directions. 

A  more  complete  discussion  of  this  jet- regime  study  appears  in 
Ref.  19  of  Appendix  A.  The  results  of  this  preliminary  study  are  presented 
here  only  for  completeness  and  were  not  used  in  the  analysis  discussed.  In 
fact,  since  all  the  data  were  collected  within  75  m  of  the  outfall,  the  data  are 
more  representative  of  the  flow- e stablishment  region  than  of  the  established 
flow  regime  of  the  jet. 
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APPENDIX  C 

FORTRAN  Listing  for  Fitting  Procedure 


The  FORTRAN  listing  of  the  computer  code  JETFIT  used  for  the  jet 
data  fitting  procedure  appears  below. 


********* ******** ***************************************************** 
PROGRAM  JETFIT 


**************************************************************** ,******<3 


PROGRAM  FITS  A  FUNCTION  (CALLED  ‘JET-FUNCTION')  WITH  12 
PARAMETERS  TO  JET  STUDY  DATA.  THE  FUNCTION  YIELDS  TEMPERATURE 
AND  VELOCITY  CENTERLINESt  CENTERLINE  DECAYS  AND  WIDTHSU/2).  THE 
CJ  NCT  T  ON  ‘JET-FUNCTION*  USES  GAUSSIAN  PROFILES,  A  QUADRATIC  FORM 
FOR  THE  CENTERLINE  TRAJECTORIES,  A  CENTERLINE  DECAY  OF  TEMPERATURE 
AND  VELOCITY  THAT  FALLS  OFF  AS  THE  INVERSE  OF  THE  SQUARE  ROOT  OF 
THE  DISTANCE  FROM  THE  OUTFALL,  AND  WIDTHS  WHICH  INCREASE  LINEARLY 
WITH  DISTANCE. 

INPUT  - 

CARD  1  (20A4) 

TITLE 

CARD  2  (  21  5 , 5F  1 0.  5 ) 

NPTS  =  NO.  OF  DATA  POINTS  (STATIONS)  (<=60) 

NLEVEL  =  NO.  OF  LEVELS  (DEPTHS)  AT  WHICH  DATA  WAS  TAKEN  (<=6» 
ANGN  =  ANGLE  OF  NORTH  W.R.T.  +  X-AX I S  (DEG.) 

BO  =  FULL  WIDTH  OF  OUTFALL  (FT.) 

8ETA0  =  ANGLE  OF  OUTFALL  W.R.T.  +X-AXIS  (DEG.) 

TO  =  OUTFALL  TEMPERATURE  (DEG.-C) 

UO  =  AVERAGE  OUTFALL  SPEED  (CM. /SEC.) 

CARD  3  (6F10.5) 

AMTEMP(1  TO  NLEVEL)  *  AMBIENT  TEMPERATURE  AT  EACH  LEVEL 

(DEG.-C ) 


CARD  4  (6F10.5) 

AMCUR ( 1  TO  NLEVEL)  =  AMBIENT  CURRENT  (CM. /SEC.)  (ASSUMFD  TO 
BE  PARALLEL  TO  THE  +X-AXIS) 

CARD  5  (6F10.5) 

CON( 1  TO  NLEVEL)  =  FACTOR  WHICH  SETS  THE  CONVERGENCE 

CRITERION  FOR  THE  SEARCH  TYPE  FIT  BASED  ON  THE  DELTA'S 
FOR  EACH  OF  THE  12  PARAMETERS  TO  BE  ENTERED  BELOW. 
CONVERGENCE  IS  ASSUMED  TO  HAVE  OCCURRED  WHEN  CHANGES 
IN  THE  PARAMETERS  ARE  ALL  LESS  THAN 

CON*DELTA(PARAMETER )  • 


CARD  6  (6F10.5  ) 

ALIMT(1  TO  NLEVEL)  =  MAXIMUM  NUMBER  OF  ITERATIONS  IN  FITTING 
THE  TEMPERATURE  PART  OF  ‘JET-FUNCTION*  (STEPS  ARE 
PRINTED  OUT  IF  AL  IMT  <  0.0). 

CARD  7  (6F10.5) 

ALIMVd  TO  NLEVFL  )  =  SAME  AS  ABOVE  BUT  FOR  THE  VELOCITY  PART 
OF  THE  FUNCTION. 

CARD  8  (6F10  .5 ) 


o  o  o  o 


o  o  o  r>  o  o  o  to  ooonooooi-joonoooooooooooooooooooooooooooooooooooo 
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PCltl  TO  NLFVEL)  =  FIRST  GUFSS  AT  FIRST  PARAMETER  OF 
•JET-FUNCTION*  (DIMENSIONLESS). 

CARD  9  (6FI0.5) 

DELP( 1 » 1  TO  NLEVEL)  =  SMALL  CHANGE  IN  P(I)  TO  BE  USED  TO 
NUMERICALLY  CALCULATE  DERIVATIVES  AND  CONVERGENCE 
CRITERION. 

CARDS  10  THROUGH  31  (6F10.5) 

REST  OF  THE  12  DI MFNSI ONLFSS  PARAMETERS  AND  THEIR  DELTA'S. 

FOLLOWING  THESE  31  PRELIMINARY  CARDS  COMES  THE  DATA  DECK  (THIS 
DECK  IS  COMPATIBLE  WITH  THE  DATA  DECK  FOR  PROGRAM  JETDAT  EXCEPT 
THAT  THE  LAST  BLANK  CARD  MUST  BE  REMOVED). 

CARD  1  (3F10.5) 

XP  =  X-COORDINATE  OF  STATION  (FT.) 

YP  =  Y-CDORDINATF  OF  STATION  (FT.) 

DEPTH  =  DEPTH  OF  WATER  ( M. ) 

CARDS  2  (5F10.5)  (DATA  AT  DIFFERENT  LEVELS  -  TWO  LEVELS  PER  CARD) 
TPMP(K)  =  TEMPERATURE  AT  K-TH  LEVEL  (DEG.-C) 

VEL(K)  =  SPEED  (CM. /SEC.) 

DIR(K)  =  DIRECTION  OF  CURRENT  W.R.T  NORTH  (DEG.) 

TEMP(K+1)  = 

VEL(  K  4-1 )  = 

DIR(K4-1)  = 

REPEAT  FOR  FACH  STATION  UP  TO  »NPTS*  STATIONS  (<=60). 


DETAILS  OF  'JET-FUNCTION* 

A  COORDINATE  SYSTEM  IS  CHOSEN  SUCH  THAT  THE  +Y-AXIS  IS  DIRECTED 
IN  THE  OFF-SHOPE  DIRECTION.  THE  +Z-AXIS  IS  DIRECTED  VERTICALLY 
UPWARD,  AND  THF  +X-AXIS  IS  SUCH  AS  TO  BE  ORTHOGONAL  TO  THE  OTHER 
TWO  AND  SO  AS  TO  FORM  A  RIGHT  HANDED  COORDINATE  SYSTEM  IN  THE 
CONVENTIONAL  SENSE  (X,Y,Z). 

THE  PARAMETERS  OF  'JET-FUNCTION'  ALONG  WITH  TYPICAL  VALUES  FOR 
SOME  POINT  REACH  UNIT  NO.  1  DATA  ARE  GIVFN  BELOW. 

P(  1)  =  A  (  1.0) 

P ( 2  )  =  ALPHA  (4.0) 

P ( 3 )  =  C  (0.5) 

P( 4)  =  GAMMA  ( 0. 33) 

P ( 5 )  =  RT  (1.0) 

P  ( 6 )  =  KT  (RANGES  FROM  ABOUT  -1.0  TO  4-1.0) 

P(  7)  =  B  (1.0) 

P ( 8 )  =  BFT  A  (4.8) 

P( 9)  =  D  (0.5 ) 

P( 10)  =  DELTA  ( 0. 19) 

0(11  )  =  RV  (1.0) 

P  (  1 2 )  =  KV  (RANGES  FROM  ABOUT  -1.0  TO  4-1.0) 


TEMPERATURE  papt  of  • jft-funct ION' 

CFNTERLINE  TRAJECTORY - 

X=XSI*COS(RT*BETAO) -(0.01* KT/ BO) *XS 1**2 *S IN ( R  T  *BET  AO ) 
C  Y=XSI*SIN(RT*BETA0)4(0.01*KT/B0)*XSI**2*C0S(RT*BETA0) 
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WHERF  XSI  IS  INTRODUCED  ONLY  DUE  TO  THE  PARAMETRIC 
FORM  OF  THE  EQUATIONS  (PENDS  LEFT  FOR  KT>0,  BENDS 
RIGHT  FOR  K  T<0 ) 

CFNTERLINF  TEMPERATURE  EXCESS  RATIO - 

(TC-T A) /(TO-TA)=A  IF  S<ALPHA*BO 

=  A*S  QRT ( AL  PHA*B0  /  S )  IF  S>ALPHA*BO 
WHERE  S  IS  THE  DISTANCE  FROM  THE  OUTFALL  MEASURED  ALONG 
THE  CENTERLINE 

TEMPERATURE  WIDTH  (TO  1/2  THE  CENTERLINE  VALUE! - 

WT/BO=C+GAMMA*S/BO 


VFLOCITY  PART  OF  'JET  FUNCTION* 

THE  CENTERLINE  TRAJECTORY  EQUATIONS  HAVE  THE  SAME  FORM  AS 
THE  ABOVE  WITH  KT  REPLACED  BY  KV  AND  RT  REPLACED  BY  RV. 

CENTERLINE  VFLOCITY  RATIO  DECAY - 

UC/UO=B  IF  S<BETA*BO 

=  B*SQRT (BETA* BO /S )  IF  S  >BET A*BO 

WHERE  UC  IS  VELOCITY  IN  EXCESS  OF  AMBIENT  CURRENT 

VELOCITY  WIOTH  (TO  1/2  THE  CENTERLINE  VALUE) - 

WU/BO=DOELTA*S/BO 


WARN! NG: 

BE  VERY  CAREFUL  OF  ATTACHING  ANY  SIGNIFICANCE  TO  THE 
INDIVIDUAL  VALUES  OF  THE  PARAMETERS  (P(l)  -  P(12)),  THEY  ARE 
USUALLY  NOT  UNIQUELY  DETERMINED  BY  THE  DATA.  DO  NOT  USE 
THESE  PARAMETERS  TO  EXTRAPOLATE  TO  VALUES  OF  S  BEYOND 
(OP  BEFORE)  THE  DATA. 


DIMENSION  TITLE ( 20) , CON ( 6) , AL I MT( 6) , AL I  MV ( 6)  ,P ( 12 , 6)  ,DE PTH ( 60 ) , BT ( 
16)  ,BV(6) ,EPSB(6),AREA(6) 

COMMON/HAVE /DEL P( 12 » 6) , TEMP ( 60 ,6  » , VE L ( 60 ,6 ) , DI R(60 ,6 ) , XP (60 ) , Y P ( 60 
1  )  ♦  AMTEMP(6  ),  AM  CUP.  (6) ,  BET A 0, BO , ANGN , TO , UO , NPT S ,DR  fRD,K 
EXTERNAL  FTSIG 
EXTERNAL  FVSIG 
DR=0. 0174532925 
R  D  =  57. 2  95780 
IM AX=  27 

FORMAT  STATEMENTS  - 

FORMAT (20A4) 

FORMAT!  21  5  »5F1 0. 5 ) 

FORM  AT ( 6F10. 5) 

FORMAT ( 3  FI 0 • 5 ) 

FORMAT! • 1PP0GRAM  JETFIT*,//,*  JET  STUDY  -  INPUT  DAT  A • , / ♦ 5X , 20 A4, // 
1  ) 

FORMAT ( *  NUMBER  OF  STATIONS  =  *,  13,  /,  '  NUMBER  OF  LEVELS  =',  12, 

1/,  •  ANGLE  OF  NORTH  W. R. T.  +X-AXIS  =  •  ,  F8.2,  *  DEG.*,  /,  •  FULL  WI 
2DT  H  OF  OUTFALL  =*♦  F7.2,  *  FT.»,  /,  •  ANGLE  OF  OUTFALL  W.R.T.  +X-A 
3X1 S  =' ,  F7 . 2 , '  DEG.*,  /,  •  OUTFALL  TEMPERATURE  (TO)  =',  F7.2,  '  DE 
4G.-CS  /,  •  AVERAGE  OUTFALL  VELOCITY  (UO)  =•  ,  F7.2,  •  CM. /SEC.*,  / 
5/,  •  LEVEL*,  2X»  *  AM  B  •  TEMP.*,  2X,  ' A  MB •  CURRENT*,  2X,  'CONV.  FACT 

6  OP •  ,  2  X ,  ' LIMIT(TEMP)* ,  2X,  • L IM IT ( V  EL ) ' ,  /,  10X,  • (DEG.-C ) • ,  4X, 

7MCM./SFC.  )  *,  /) 
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96  F0PMAT(4X,  !2,  5X,  F7.2,  6X,  F8.2,  6Xf  F8.5,  6X,  F7.1,  5X,  F7.1) 

97  FORMAT ( /// » '  INITIAL  PARAMETERS  OF  JET-FUNCTION  (ALL  ARE  DIMENSION 
1LFSS  )',//*  *  LEVEL',8X,'P(1)',8X,'P(2)',8X,'P(3)',8X,'P(4)',8X,'P(5 
2  )•  ,8X,»  P(6)  '  ,/) 

98  FORMAT! 4X, I  2) 

99  F0RMAT(6X,6(2X, F10.6)  ) 

900  FORMAT! //,'  LEVEL • , 8X , • P ( 7 ) ' , 8X , • P ( 8 ) • , 8X , *P ( 9  )  • , 7X,  *P ! 101 ' » 7X, • P ( 
1111  *,7X, 'P( 12) • ,  /) 

901  FORMAT (• 10ATA  DECK ' // , 5X , 20A4, // , 1  STAT I  ON ' , 14X, • X ' , 1 1 X , » Y* , 3X, • DE 
1PTH*  »2X  » '  LEVEL* ,4X,' TEMP.' , 7X ,  • S PEED* , 2X ,  • DI RECT ION *,/♦ 1 8X, MET.) • 

2,  7X» '  ( FT. ) ' » 4X,  MM.) • , 8X , • ( DEG.-C 1 • ,2X,' (CM. /SEC. )  •  ,5X,' (DEG. )' ,/ ) 

902  FORM  AT ( 5X  t  13,  5X,  F10.2,  2X,  F10.2,  2X,  F6.2) 

903  F0RMAT(48X,  12,  5X,  * - ',  8X,  * - *,  7X,  ' - ') 

904  FORMAT!  48X,  12,  2X,  F  7.  2 , 5X ,  F  7.  2 , 4X ,  F  7.  2 ) 

905  FORMAT (• 1LEVEL* , 12, 3X, ' (TEMPERATURE  FIT ) *,//,5X, 20 A  4, // I 

906  FORMAT!//, •  FINAL  TEMPERATURE  RESULTS  F0R',I2,'-TH  LEVEL ',//, 8X, ' S 
1TGMA  =',  F8 • 3, '  AFTER*, 16,*  ITERATIONS.',//) 

907  F0RMAT(10X,» P(*  ,12  ,*  I  =*,F10.6) 

908  FORMAT! • 1C0MPARI SON  OF  DATA  AND  FIT  RESULTS  FOR  THE*,I2,'-TH  LEVEL 
1 .*,//, 5 X,20A4, //,2X, 'STATION', 1 IX, *X', 1 IX, • Y ' , 2X, ' TEMP. ( DA TA ) • ,2X, 
2 'TEMP. (CALC. )•  ,3X,*  BET  AT (CALC. ) • ,2X , • S PEED ( DAT A ) ' , 2X , 'SPEED! CALC. » 
3',3X, ' BET AV( CALC. ) • ,/, 16X, ' (FT. ) * ,7X,» (FT. ) • ,5X,' (DFG.-C )• ,6X,» (DE 
4G.-C ) ' ,°X , • (DEG. ) ',3X, • ( CM . / S EC . ) ' , 4X, ' ( CM . /SEC. ) ',9X, '(DEG. ) • ,/) 

909  FORMAT ( 6X,I 3,2X,F10. 2 , 2 X , F 1 0. 2 , 6 X , F 7. 2 , 7X , F7 . 2 , 8X , F7 .2 , 6X , F7. 2 , 7X , 
1F7.2,8X,F7.2) 

910  FORMAT ( // , •  TABLE  OF  JET-FUNCTION  (TEMPERATURE  PART )',//, 9X, • XSI ' , 
1 11X, 'S' ,10X, ' XC* ,10X, • YC' ,2X ,' ( TC-TA )/( TO-TA > • ,2X , ' H ( 1/ 2 ) / BO • , 5X , • 
2BFTA',/,7X,'(FT.)',7X, '(FT.) *,7X, *(FT.)',7X,'(FT.) • ,31 X, ' ( DEG. ) • ,/ 
3) 

911  FORMAT ( 2X,F10. 2,2 X,F 10.2 ,2X,F 10.2 ,2X  ,F 10.2  ,  9X ,F8. 4 , 3X , F8 . 3 , 2X , F7 . 
12) 

912  FORMAT!//, •  TABLE  OF  JET-FUNCTION  (VELOCITY  PART )',//, 9X, • XSI ', 1 IX 
1,  'S',  10X,  'XC  ',  10X,  'YC  •  ,4X,  'UC/UO*  ,2X,«  (UC  +  UA)  /UO'  ,2X  ,'  W(  1/2  )/B0'  ,5 
2X,' BETA* ,/,7X, * (FT. ) • ,7X, • (FT.) «,7X, '(FT. )' , 7X, '(FT. ) ' ,35X, '(DEG. ) 
3*  ,  /) 

913  FORMAT!  2X,F10.2,  2X,F 10.2, 2X,F 10. 2, 2X ,E 1 0. 2 , 2 X ,F 7. 4 ,4X ,F8 .4 ,3X ,F8 . 3 
1  ,2X,F7.2) 

914  FORMAT!/, •  TOO  MANY  DATA  POINTS.',/) 

915  FORMAT!/, •  TOO  MANY  LEVELS.',/) 

916  FORMAT (• 1LEVEL* , 12, 3X ,' (VFLOCITY  F IT )',//, 5X , 20A4, // ) 

917  FORMAT!//, •  FINAL  VELOCITY  RESULTS  FOR',  12,  '-TH  LEVEL',  //,  8X, 

l'SIGMA  =',  F8.3,  •  AFTER',  16,  •  ITERATIONS.',  //) 

918  FORMAT! • 1APPP0XIMATE  EXCESS  TEMPERATURE  ISOTHERM  AREAS  AND  JET  DEP 
1THS.',  //,  5X,  20A  4 ,  ///,  2X ,  'TEMP.  EXCESS',  6X,  ' AREA! ACRES  )  -  L 
2EVEL' , / ,6X, ' (C-DEG.) ',  IIX,  '1',  11X,  *2',  11X,  *3',  11X,  '4',  11X 

3,  '5',  11X,  *6',  /) 

919  FORMAT!  7X, F7.2, 6( 2X, F10. 4)  ) 

920  FORMAT ( / / / / ,  11X,  'S',  2X ,  'TEMP.  EXCESS',  4X,  'DEPTH',  6X,  'E  TA ( F 

IT.)  -  LFVEL',  /,  7X ,  '(FT.)',  6X,  '(C-DEG.)',  2X ,  '(LEVEL)',  9X,  • 
21',  9X »  *2',  9X,  *3',  9X ,  '4',  9X,  *5',  9X,  *6',  /) 

921  FORMAT (2X,F10.2 ,7X,F7 .2  ,  2X  ,  F7.3,  6(2X,  F8.2)) 

922  FORMAT!////,  11X,  'S',  3X,  'VEL.  EXCESS',  4X,  'DEPTH',  6X,  ' ET A ( FT 

1.)  -  LFVEL',  /,  7X ,  '(FT.)',  4X,  • ( CM. /SEC .  )  •  ,  2X,  '(LEVEL)',  9X, 

2*1',  9X ,  '2',  9X ,  '3',  9  X ,  *4',  9X ,  '5',  9X,  *6',  /) 

923  FORMAT!//,  '  GOODNESS  OF  FIT:',  /,  20X,  'TEMPERATURE',  F8.2,  •  *', 
1/,  20X ,  'VELOCITY',  F11.2,  •  ?•  ) 

C 

C  R  FAD  IN  DATA 

33  RFAD(5 ,90»END=30 ) (T I T L F ( I ) ,  1= 1 , 20 ) 

REAO( 5,91)  NPTS,NLEVFL , ANGN , 80 , BETAO , TO ,U0 
IF ( NPTS -60 ) 1 , 1,31 
1  IF(NLFVFL-6)7,7,32 
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7  RFAD(5,92)  ! AMTEMP! K ) , K  =  1 ,6 ) 

RF  AD  ( 5,92)  ( AMCUR ( K ) , K  =  1 , 6 ) 

RF  AD ( 5  »  92  )  (CON(K) ,K=1 ,6 ) 

READ  (  5,  92  )  ( AL IMT(K) ,K=1,6) 

RFAD(5,92)  ( AL  IMV ( K  ) , K  =  1 , 6 ) 

HO  2  J  =1 » 1 2 

R  EAP ( 5 »  °2 )  (P( J,K)  ,K=1,6) 

RF  AD ( 5 , 92 )  (DELP(J,K),K=1»6) 

2  CONTI NUF 

DO  6  1=1, NPTS 

3  R  FAD ( 5  *  93 )  X P(I) , YP ( I) , DEPT H( II 
IF  (  XPII  ).NE.  0.  0)  GOTO  4 
IF(YP(  I  ).NE.0.0)  GOTO  4 
IF(OEPTH(I).NE.O.O)  GOTO  4 
GOTO  3 

4  ANUM=NLEV  EL 

ANUM=ANUM/2.0+0.75 
M= ANUM 

N=2*N 

DO  5  K=1 » N ,2 

5  READ! 5, 92)  TPMP( I  , K  )  ,  VE L ( I ,K > ,D I R (  I  ,  K)  ,TE MP! I , K+l ) ,  VEL ( I , K+l ),DIR( 
1  I, K  +  l ) 

6  CONTINUE 

c 

C  PRINT  OUT  DATA 

C 

WR I TE ( 6 »  94 )  (TITLEU  )  ,1*1,20) 

WRITF(6,95  )  NPTS, NLEVEL, ANGN , BO* BETAOt TO, UO 
DO  8  K  =  1  ,  NLE  VE  L 

8  WR ITE ( 6  »  96 )  K , A M TE MP ( K ) , A MC UR ( K ) , CON ( K) , AL I MT ( K) , AL I  MV ( K I 
WRITEI6  ,97  ) 

DO  9  K  =  1  , NLEVEL 
WR ITE I  6, 98  )  K 

WR I T  F ( 6  ,99  )  (P(J,K),J=1,6) 

WR I TE ( 6  » 99 )  (DFLP! J,K) , J=1 ,6) 

9  CONTINUE 

WR  I  TF ( 6  ,900 ) 

DO  10  K=l, NLEVEL 
WR  ITE ( 6, 98 )  K 

WR  I  TE ( 6  ,99 )  (P(J,K),J=7,12) 

WRITE! 6, 99)  ( DE L P ( J , K )  , J  =  7 , 1 2 ) 

13  CONTINUE 

WR I T  E ( 6 , 901 )  (TITLE! I) ,1=1 ,20) 

DO  11  1=1, NPTS 

WRITE (6, 902  )  I, XP ( I  ) , YP ( I ), DEPTH! I ) 

DO  11  K=1 , NLEVEL 
IF!TEMP( I , K )  )  1 3 , 14, 14 

13  IE (VEL (  I  ,K ) )15,14,14 

15  WR  I  TF  (  6 , 903  )  K 

GOTO  11 

14  WR ITE! 6 ,904 )  K , TE MP ! I , K ) , V EL (  I , K ) , D IR ( I , K ) 

11  CONTINUE 

C 

C  START  LOOP  THROUGH  LEVELS  (COMPLETE  ONE  LEVEL  AT  A  TIME). 

C 

DO  16  K=l, NLEVEL 

C  TF(ALTMT(k  ).EQ.O.O.AND.ALIMV(K).EQ.O.O)  GOTO  16 

M  =  6 
C 

C  FIT  TO  TEMPERATURE  DATA. 

C 

S I  GT  =  0.  0 
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OELMAX=0.0 
no  17  J=1,M 
BT< J)=P( JtK) 

OELMAX  =  DELMAX  +  DELP(J,K  ) 

17  EPSBf J) =OELP( J,K)*CON(K  ) 

DELMAX=0. 16666667*DFLMAX 
L=AL  IMT  (K  ) 

C  IF(L.EQ.O)  GOTO  12 

WRITE< 6,905)  K ,  (  TI TIE ( I)  , I =1 , 20) 

CALL  GMIN(FTSIG,M, BT,S IGT,EPSB,L, ITER,DELMAX) 
CALL  L  A  BEL ( ITER) 

WP ITE( 6, 906 )  K,SIGT,L 
DO  18  J= 1 , M 
P( J,K)=BT( J) 

18  WP  ITE ( 6 , 907 )  J,BT(J) 

WP I T E ( 6  ,910  ) 

no  24  I =1 , 1  MAX 
ANUM=I-1 
XS I=ANUM*50.0 
ANUM=0.01*BT(6) /BO 
THETA=OR*BT( 5)*BFTA0 
STH=SIN(THETA) 

CTH=COS( THETA) 

A1=2.0*XSI*ABS(ANUM  ) 

I F ( A 1-0. 1)34,34,35 

34  S=XSI*( 1.0+0. 166666667*A1*A1-0.025*A1**4) 

GOTO  36 

35  A2=SQRTU.0+A1*A1  ) 

S=0.5*XSI*< A2+AL0GC A 1+A2) /A1 ) 

36  X=XST*CTH-ANUM*XSI*XSI*STH 

y=xsi*sth+anum*xsi*xsi*cth 

T=PT( 1 ) 

IF (S-BT (2)* BO) 25,25,26 

26  T=T*SQRT(BT( 2)*B0/S) 

25  A1=STH+2.0*ANUM*XSI*CTH 

A2=CTH-2.0*ANUM*XSI*STH 
BFTAT=RD*ATAN2( A1 ,A2) 

WBO=BT (3)+BT(4)*S/B0 

24  WRITF(6,911)  XS I  ,  S  ,  X  ,  Y  ,  T  ,  W  BO  ,  BET  AT 

C 

C  FIT  TO  VELOCITY  DATA. 

C 

12  SI  GV  =  0.  0 

HELM AX  =  0  .0 
00  19  J=1 ,M 
BV(J)=P(J+M,K) 

OELMAX=OELMAX+DELP( J+M,K ) 

19  EPSR( J)=DELP( J+M,K)*CON( K) 

PEL M A X= 0.1 666666 7+OEL MAX 
L=  AL I  MV ( K  ) 

C  T F  f L.EQ.O)  GOTO  52 

WRITE(6,916)  K, ( TITLE ( I ) ,1=1,20) 

CALL  GMIN(FVSIG,M,BV,SIGV,EPSB,L, ITER, DEL  MAX  I 
CALL  L A  BEL ( I TEP ) 

WRITEt6,917)  K,SIGV,L 
DO  20  J=7 , 1 2 
P(  J,K)  =BV(  J-M) 

23  WR!TE(6, 90  7  )  J,BV(J-M) 

WPITF(6,912 ) 

DO  27  I  =  1 , T  MA  X 
ANUM= 1-1 
XSI=ANtJM*50.0 
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ANUM=0.01*BV(6)/B0 
THET  A=  OR  *BV ( 5 ) *BETAO 
STH=  SI  N  (  THETA) 

CTH=COS(THETA) 

A1=2.0*XSI*ABS (ANUM) 

IE(A  1-0.1) 37t37,38 

37  S  =  XS  I*( 1.0+0. 16666667*A1*A1-0.025*A1**4) 

GOTO  39 

38  A2=S0RT< 1.0+A1*A1) 

S=0.5*XST*< A2+AL0G( A1+A2J/A1 ) 

39  X=XSI*CTH-ANUM*XSI*XSI*STH 

y=xsi*sth+anum*xsi*xsi*c th 

VC=BV(1  ) 

I E ( S-BV(2)*B0)28,28,29 
29  VC=VC*SQRT(BV( 2)*BO/S) 

28  A1=STH+2.0*ANUM*XSI*CTH 

A  2=C  TH-2.  0*A  NUM*  X  SI  *  STH 
BETAV=ATAN2(A1,A?) 

VD=VC+COS(BETAV)*AMCUR(K)/UO 

BETAV=RD*BFTAV 

W80=  BV ( 3 )  +  BV (4 ) *S/B0 

27  WR  I  TE  (  6 , 91  3  )  XS I  ,  S ,  X  ,  Y  ,  VC  ,  VD  ,WBO  ,  BET  AV 

C 

C  PRINT  OUT  TABLE  OF  COMPARISON  OF  CALCULATIONS  FROM  'JET-FUNCTION' 

C  TO  DATA  FOR  THIS  LEVEL. 

r 

52  WP  ITE (6 » 908  )  K , < T  ITL E (I) , I  *  1 , 20 1 
A  NUM=0. 0 
A1=0.0 
A2=0.0 

DO  21  I =1 1 NPTS 

IF  I  TEMP ( I ,  K  ) >  22,  23,23 

22  I F ( VCL ( I ,K) )21, 23,23 

23  CALL  TEMVEL(XP(I  )  ,YPII)  ,  BT  ,  BO  ,  BET  AO ,  RT  ,  BET  AT  ) 

T=RT*TO-RT*AMTEMP(K )+AMTEMP(K) 

BETAT=BETAT*RD 

A  1  =  A  1+ ( TEMP ( I ,K) -AMTEMP(K) 1**2 

CALL  T  EMVEL ( XP ( I ) , YP ( I ) , BV, BO,BETAO,RV,BETA VI 

OEVBET=( ANGN-DIRt I,K) )*DR-PETAV 

VD=VEL< I ,K)*COS(OEVBET) 

VC=RV*UO+AMCUR(K )*COS( BETAV) 

BETAV=RD*BFTAV 
A2=A  2+VD* VD 
ANUM= ANUM  +  1  .0 

WR  I  TE  (  6,909)  I  ,  X  P  (I)  ,  Y  P  ( I  )  ,T  EMP  ( I  ,  K )  ,  T  ,  BET  AT  ,  V  D,  VC ,  BETAV 
21  CONTINUE 

ANUM  =  A  NUM-6 .0 

IFIANUM.LT. 0.5)  ANUM=1.0 

A 1  =  100.0* ( 1.0-SI GT/SORTI A 1/ANUM) ) 

A2  =100.0* ( i ,0-S I GV/SORT ( A2/ ANUM ) ) 

WP  I  TE  (  6,  923)  A1,A2 
C 

C  END  OF  LOOP  THROUGH  LEVELS. 

C 

16  CONTINUE 

C 

C  CALCULATES  ISOTHERM  AREAS  AND  JET  DEPTHS. 

C  DEPTH  IS  DEFINED  AS  THE  EQUIVALENT  TWO-DIMENSIONAL  DEPTH. 

C 

WRITE!  6,918)  I  TI  TLE  (  I  )  »I  =1  ,20) 

IT  =TO-AMTEMP (NLEVEL ) 

41  TE  =  1 T 


145 


no  40  K=l,  NLEVEL 
40  AREA(K)=TARFA(TE,P) 

WRITE (6, 91 9)  tp,  (are  AC  K)  ,K=1, NLEVEL) 

I T= I T- 1 

IF ( IT.GE.l )  GOTO  41 
WR I TE ( 6  »  92  0 ) 

DO  42  1  =  1 » IMA X 
ANUM  = 1-1 
S=50.0*ANUM 

T=  P  (  1,  1)*{  TO-A  MTE  MP(  1)) 

IF(S.LE.P(2,1)*B0)  GOTO  43 
T=T*SQRT( P  (  2  » 1 )  *  BO/S ) 

43  T=0.5*T 

00  44  K=1 »  NLFVEL 

44  ARFA(K)=ETAT(S,T,P) 

A  1  =  0  .0 

IF(AREA(1).LT.(0.1*B0))  GOTO  46 
A  1=0. 5*ARE  A ( 1) 

00  45  K=2, NLFVEL 

45  A1=A1+AREA(K) 

A1=A1/AREA(1)+1.0 

46  WRITF(6,921)  S , T , A 1 ♦ < AR E A ( K ) , K=1 , NLE VF L ) 

42  CONTINUE 

WPTTFI 6, 922) 

00  47  1=1, IMAX 
ANUM= I - 1 
S=  50.0*ANUM 
VC  =  P  C 7 , 1 )*U0 

I F ( S.LF.Pl 8,1 )*B0)  GOTO  48 
VC=VC*5QRT(P( 8, l)*B0/S) 

48  V0=0.5*VC 

DO  49  K=l, NLFVEL 

49  AREAIK  )  =  ET A V ( S , VO , P ) 

A1=0.0 

!F(AREA( 1).LT. (0. 1*B0) )  GOTO  51 
A1=0.5*ARFA< 1) 

00  50  K=2 , NLEVEL 

50  A1=A14-ARFA(K) 

A1  =  A1/ARFA(1  14-1.0 

51  WR I TF ( 6 , 92 1 )  S , VO , A1 , < AR E A( K ) , K=1 , NL EVEL ) 

47  CONTINUE 

c 

GOTO  33 
C 

31  WR  ITF(  6,914) 

GOTO  30 

32  WRITE!  6,915) 

33  CALL  EXIT 
STOP 

END 

C  C 

SUBROUTINE  G  MI  N<  F  UNC  ,  M,X  ,  F  ,  EPS  ,  L I M  IT  ,  IT  EP ,  DELMAX  ) 

C  C 

£+ +4 **+****++*+++*++*+**+**++*+*+***+* *+*++**+*+*++++*++++**+++*+*+**++0' 

C 

C  ROUTINE  FINDS  THE  MINIMUM  BY  A  SEARCH  TECHNIQUE  OF  THE 

C  FUNCTION  (CALCULATED  IN  SUBROUTINE  FUNC)  WITH  RESPECT  TO  ITS 

C  PARAMETERS  X( I ) ,  1=1, M.  FIRST  THE  DIRECTION  OF  STEEPEST 

C  DECENT  IS  SELECTED  BY  EVALUATING  THF  GPADIFNT  OF  THE  FUNCTION 

C  WITH  RESPECT  TO  THE  M  PARAMETERS.  THEN  A  STEP  IS  TAKEN  ALONG 

C  THAT  DIRECTION.  THE  SIZF  OF  THE  STEP  IS  PRESET  BUT  MULTIPLES  OR 
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C  FRACTIONS  OF  THIS  STEP  SIZE  ARE  USED  IF  NECESSARY. 

C 

DIMENSION  XI 10) ,G( 10) ,EPS( 10) ,XZ(10) »0UMI10) 

EXTERNAL  FUNC 
C 

L  =  0 
LP=0 

IF(LIMIT.FO.O)  GOTO  35 
IP  =  1 

TFILIMIT.LT. 0)  IP  =  0 
LI  MI T= I ABS (LIMIT ) 

I F I  I P. FO. 0. OR. LIMIT.GE.100)  WRITE (6, 2) 

2  FORMAT ( 17X» *FUNC',5X, ’PARAMETERS',/) 

C 

I  I  TER  =1 
L=L  +  1 
LP=LP+1 

IFILP.GE. 100)  LP  =  0 
CALL  FUNC (M, X, F, G, ITER) 

IF(IP.EO.O.OR.LP.EQ.O)  WRITE(6,3)  F  ,  I  X  (  J  )  ,  J=  1,  M) 

3  FORMAT!  6X,  7F  15.  8) 

IFIITER.EQ.O)  GOTO  15 
I  TER  =2 

GN0RM=0.0 
DO  4  1=1, M 

4  GNORM=GNORM*G( I)*G! I ) 

GNORM=  SORT! GNORM ) 

C  IF(IP.FQ.O.OR.LP.FO.O)  WRITE(6»3)  GNORM ,( G t J ), J= 1 , M) 

TF(GNORM.LE.O.O)  GOTO  15 
C 

DO  5  1=1, M 

G I  I ) =G(  I)*DEL MAX /GNORM 

5  XZ(I)=X(I)-G(I) 

CALL  FUNC(M,XZ,FZ,DUM, ITER ) 

IFIFZ.LT. F)  GOTO  40 
C 

32  KK=0 

DO  33  I =1 , M 
G( I )=0.5*G( I  ) 

IFIABSIGI I  )  )  .GT.EPSII  )  )  KK=1 

33  XZ(I)=X(I)-G(I) 

CALL  FUNC (M,XZ»FZ,DUM,ITER) 

IFIFZ.LT. F)  GOTO  13 
TFIKK.EQ.l)  GOTO  32 
IT  ER  =  3 
GOTO  15 

13  DO  11  1=1, M 

II  X(I)=XZ(I) 

F  =  FZ 

GOTO  12 
C 

40  F=FZ 

DO  41  1=1, M 
XI  I  )  =  XZ(  I  ) 

41  xzm=xm-G(i) 

CALL  FUNC I M,XZ,FZ»DUM,ITER) 

IFIFZ.LT. F)  GOTO  40 
C 

12  IFIL.LT. LIMIT)  GOTO  1 

IFIITER.EQ.O)  GOTO  15 
TTER=4 
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15  L  T  MI T=L 

RFTURN 
C 

c 

35  ITER  =  2 

CALL  FUNC!M,X,F,DUM, ITER ) 

GOTO  12 
C 

FNO 

C ******************* ********************** *****************************0 


C  C 

SUBROUTINE  LABEL (ITFR) 

C  C 

c 

I=ITER+1 

GOTO! 5, 10, 15,20,25) , I 

5  WRITE!6,6) 

6  PORM  AT !  /  /  » '  ARGUMENT  OUT  OF  RANGE.') 

RETURN 

10  WP IT  E 1 6 , 1 1  ) 

11  PORMATC//,* 

RFTURN 

15  WR  ITE (6, 16 ) 

16  FORMAT!//, • 

RETURN 

20  WRITF(6,21) 

21  PORMAT!//,' 

RFTURN 

25  WRITE(6,26) 

26  FORMAT!//,' 

RETURN 


SUCCESSIVE  ITERATIONS  NOT  REDUCING  FUNCTION.') 

GRADIENT  TOO  SMALL  FOR  USEFUL  ITERATION.') 

ERROR  CRITERIA  SATISFIED.') 

MAXIMUM  NUMBER  OF  ITERATIONS  ATTAINED.') 


END 

C*  4c  *  *  *  **  *  *  *  *  **  *  *  ************  **  ************************  *****************C 

c  c 

SUBROUT INE  T PMVEL ! X, Y, B, BO,BETAO,R ,RBETA) 

C  C 


c 

C  ROUTINE  RETURNS  R  ILOCAL  EXCESS  TEMPERATURE  OR  VELOCITY  RATIO) 

C  AT  THE  POINT  IX,Y)  AS  DETERMINED  BY  'JET-FUNCTION*  WITH  PARAMETERS 

C  B ( I ) ,  OUTFALL  FULL  WIDTH  BO,  AND  OUTFALL  ANGLE  BE  TAO  W.R.T.  THE 

C  X-AXIS.  ROUTINE  ALSO  RETURNS  R  BET  A  !T  HE  ANGLE  IN  RADIANS  THE 

C  CENTERLINE  ASSOCIATED  WITH  THE  POINT  !X,Y)  MAKES  WITH  THE  X-AXIS). 

C 


DIMENSION  B I  6  » 


AK=0.01*B!6)/B0 
ABSA  K  =  AB  S ( AK ) 

THFTA=0. 0174532925*6! 5)*BETA0 
STH=SIN!THETA) 

CTH=COS( THETA) 


XO=X*STH-Y*CTH 

yo=x*cth+y*sth 

W=  27,  0*AK*AK*Y0*Y0+2.0* (  1. 0+ 2. 0* AK* XO ) **3 
IF  !Y0.LT.0.0 )  GOTO  3 
C 

C  FOLLOWING  FOR  CASE  WHERE  YO.GE.O.O 

C 

TFIW.LT.O.O)  GOTO  2 

A1=0,25*Y0*ABSAK 

A2=S0RT (0.002314R148148*W) 
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IF(A1.LT.0.03*A2)  GOTO  1 
XSI=(CUBRT(A1+A2)+CUBRT< A1-A2) ) /ABSAK 
GOTO  4 

1  XSI=YO*CUBRT (2 .O/W) 

GOTO  4 

2  A1=SQRT(-1.0-2.0*AK*X0) 

S=( 3.67423461 3*YO*ABSAK)/<Al*Ai*Al ) 

S=ARCOS( S)*0. 333333333 

XS  1= ( 0 • 81 64965809*A l*COS ( S ) I /ABSAK 

GOTO  4 
C 

r.  FOLLOWING  IS  FOR  THE  CASE  WHERE  YO.LT.O.O 

C 

3  IF(W.GT.O.O)  GOTO  12 
GOTO  2 

C 

C  NOW  HAVE  VALUE  FOR  XSI. 

C 

4  A1=2.0*XSI*ABSAK 
IF(  A  1-0.  1)  5,5,6 

5  S  =  XS  I*<  1.0  +  0. 166 666667* A 1*A 1-0. 02 5*A 1**4) 

GOTO  7 

6  A2=S0RT(  1.  0+Al*A  1) 

S=0.5*XS I*(A2+AL0G<A1+A2)/A1) 

C 

C  NOW  HA VF  VALUE  FOR  S  • 

C 

7  PTASQ=(XO+AK*XSI*XSI )**2*(Y0-XSI )**2 
A1=STH+2.0*AK*XSI*C TH 

A2=CTH-2 .0*AK*XS I *STH 
RBFTA=ATAN2( A1  ,A2 ) 

W=B( 3)*B0*B( 4)*S 
IF(W.LE.O.O)  GOTO  8 
A  1=0. 6931 471 806*ETASQ/(W*W) 

IF(A l.GE.25.0)  GOTO  8 
A1=FXP(-A1 ) 

GOTO  9 

8  A  1  =  0  .0 

9  IF<S-B(2)*B0)10,10,11 

10  R  =B(  1)  *A  1 
RETURN 

11  R=B(l)*Ai*SORT (B(2)*B0/S ) 

RETURN 

12  R  =  0 .0 

rbeta=theta 

RETURN 

END 

c *  * ** * * * * * * * * * * *  * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * c 

c  c 

FUNCTION  CUBRTIX) 

C  C 

r*************************************************$********************C 

C 

C  ROUTINE  CALCULATES  THE  CUBE-ROOT  OF  A  REAL  NUMBER. 

C 

Y=  X 

SI GN=1 • 0 
IF ( Y ) 1,3,2 

1  S I GN=- 1.0 
Y =-Y 

2  Y=0.  333333333*  ALOGI  Y) 

CUBRT=SIGN*FXP(Y) 


oooo  >->  is)  o  r>  o  o  r>  o  r»  o  t—  i\»  o  oc~»or»oooor>r» 
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RETURN 

3  CUBRT  =  0.  0 

RETURN 
FND 

*  *  *  *  4c  *  *  *  *  *  *  *  *  *  *  4c  *  *  *  *  *  *  **  *  *  *  ***  *  *  *  *  *  *  *  *  *  *  *  *  **  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  C 

c 

SUBROUT  I NF  FTS I G ( M, BT , S I GT , G, ITER ) 

* ************ ************** ****** ************** ft********** ************ 

ROUTINE  CALCULATES  THE  STANOARD  DEVIATION  OF  THE  TEMPERATURE  DATA 
FROM  ’JET-FUNCTION'  WITH  PARAMETERS  BT.  ON  CERTAIN  CALLS,  THE 
DERIVITIVFS  OF  THIS  DEVIATION  WITH  RESPECT  TO  EACH  PARAMETER  IS 
ALSO  CALCULATED. 

DIMENSION  BT(6) ,G<6) ,BB<6) 

COMMON /HAVE /DELP< 12,6) , TEMP (60,6) , VE L (60  ,6 ) , 01 R C 60 ,6 ) , XP <60 ) , YP ( 60 
1 ), AMTPMPC& ), AM  CUP  ( 6 ) » BET AO, BO, ANGN, TO, UO,NPTS,DR »RD  »K 

CALL  TSIG(BT,SIGT) 

IFt  ITER.EQ.2  )  RETURN 
DO  I  J  =  1 , M 
DO  2  1=1, M 
BB(I)=BTt I ) 

BR(  J)=BT<  JHDELPC  J,K) 

CALL  TSIGI BB,SIGP > 

BB ( J ) =BT ( J )-DEL  P ( J  »  K  ) 

CALL  TSIG(BB,SIGM) 

G( J)=0.5*( SIGP-STGM) /DELPI J,K) 

RFTURN 

END 

***************************************************************** ***+*£ 

C 

SUBROUTINE  FVSIG(M,BV,SIGV,G, ITER) 

C 

*********************************************************** ******** ***^ 

ROUTINE  DOES  THE  SAME  THING  AS  FTSIG  EXCEPT  FOR  THE  VELOCITY 
PART  OF  'JET-FUNCTION*  AND  THE  VELOCITY  DATA. 

DIMFNSTON  BV<6) ,G<6) ,BB< 6) 

COMMON/HAVE/DFLP( 12,6) ,TEMP( 60,6), VEL (  60,  6  ) ,  D  IR  (  60 , 6 )  ,  XP  <  60 )  ,  YP (  60 
1 ) ,AMTFMP( 6) ,AMCUR (6) , BET  AO , BO , ANGN , TO , UO , NPTS , DR , RD, K 

CALL  VS  I G( BV ,S IGV  ) 

IF ( I TER.EQ. 2)  RFTURN 
DO  1  J=1,M 
DO  2  1=1 ,M 
BB  (  I  )=BV(  I  ) 

BB(J)=BV(J  )+DELP(J*-M,K  ) 

CALL  VSIG( BB,S IGP) 

BB( J)=BV( J)-DELP< J+M,K) 

CALL  VSIGCBB,SIGM) 

G( J) =0.5* (S IGP-SIGM)/DELP(J+M,K) 

return 

FND 

********************************************************************* *c 

c 

SUBROUTINE  TS I G( BT»  S I GT ) 

C 

******************************************** ** *********************** *(^ 


o  o 


r>  r>  o  o  o  o  o  r>  i\»  t—  r>  o  <->  r>  o  o  o  r>  i\>  »—  o  oooo 
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ROUTINE  CALCULATES  THE  STANDARD  DEVIATION  OF  THE  TEMPERATURE 
DATA  FROM  'JET-FUNCTION*  OF  PARAMETERS  BT. 

DIMENSION  BT (6  I 

COMMON /HAVE /DEL P( 1 2 , 6 ) ,TE MP ( 60 ,6  I , VE L (60 , 6 ) , DI R < 60 ,6 ) , XP ( 60 1 ,  YPt  60 
l ) ,  AMT  EMP ( 6  ) , AMCURI6 ), BET  AO , BO , ANGN, TO, UO, NPTS , DR ,RD,K 

SIGT=0.0 

ANUM=0.0 

DO  2  1=1, NPTS 

I F ( TEMP ( I , K ) >2,1,1 

CALL  TFMVFL(XP( I >,YP( I ) , BT , BO, BETAO, RT , RBETA ) 

SIGT=SIGT*(TEMP(  I  ,  K) -RT*TO  +  PT* AMT EMP <K ) -AMT FMP <K  )  1**2 

ANUM  =  A  NUM  + 1 • 0 

CONTTNUF 

DEGFR=ANUM-6.0 

IF(DEGFR .LT.0.51  DEGFR= 1 . 0 

SI GT=SQRTCSIGT/DEGFR) 

RF  TURN 
END 

*********************************************************************0 

c 

SUBROUTINE  VSI G( B V, SIGV) 

C 

******************************************* **************************0 


ROUTINE  OOES  THE  SAME  THING  AS  TSIG  ONLY  FOR  VELOCITY, 

DIMENSION  B V( 6 ) 

COMMON/ HAVE/DELP< 12, 6  I, TEMPI  60,6) ,VELt  60,6), DIR ( 60,6) , XPI60)  ,YP(60 
1)  ,AMTEMP(6)  ,  AMCUR  (6  )  ,  BET  AO  ,  BO  ,  ANGN,  TO ,  UO,  NPTS  ,  DR  ,  R  D,  K 

S I GV  =  0  .0 
ANUM=0, 0 
DO  2  1=1, NPTS 
IF (VEL ( I , K )  ) 2 , 1 , 1 

CALL  TEMVFLIXPU ) ,YP(I) , BV , BO, BET  AO, RV, RBETA ) 

DEV BET  = ( A NGN-0 IR ( I,K) ) *DR-R  BETA 

SIGV=SIGV+(VEL(I,K)*COS( DEV BET )-R V*UO~AMCUR ( K ) *COS ( R BETA ) )**2 

ANUM=ANUM*1. 0 

CONT  TNUE 

DEGFR=ANUM-6,0 

IF (DEGFP.LT ,0,5)  DEGFR=1 ,0 

S IGV=SQRT (SIGV/DEGFR ) 

RETURN 

CNQ 

****  ************************************************** ***************Q 

C 

FUNCTION  TAR  EA ( TE , P ) 

C 

************************************************************** *******0 

FUNCTION  CALCULATES  THE  APPROXIMATE  ISOTHERM  AREAS  OF  'JET- 
FUNCTION'  NEGLECTING  THE  CURVATURE  OF  THE  CENTERLINE. 

DIMENSION  P ( 12 , 6  ) 

COMMON/HAVE/DE LP( 12 ,6) , TEMP < 60 ,6 ) , V EL ( 60 , 6 ) , D I P < 60 , 6 ) , XP ( 60 ) , YP (  60 
1  ), AMTEMPC6 ), AMCUR <  6) , BET AO, BO, ANGN , TO, UO , NPTS , DR  »RD  »  K 

TARE A=0, 0 

R=P( 1,K )*( TO-AMTEMP(K) ) /TE 
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IF  f  R  .LT  .1.01  GOTO  10 

SMAX=P<  2 » K )  *BO*R*P 

NSTP=0.1*SMAX 

I F  (  NST  P.  GT  .  100  )  NSTP=100 

A I =NSTP 

DELS=SMAX/AI 

S=0.0 

T ARE  A=  0.  5*ETAT(  S,TE  ,P) 

NSTP=NSTP-1 
00  11  1=1 » NSTP 
A  1=  I 

S=AI*DELS 

11  TAREA=TAREA+ETAT(S,TE,P) 

TARFA=DELS*TAREA*4.592c-5 

10  RETURN 

END 

£**********************************************************************0 

c  c 

FUNCTION  ETATI  S,  Tf  ,  p) 

C  C 

c  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ♦  *  *  *  c 

c 

C  FUNCTION  CALCULATES  THE  LATERAL  DISTANCE  OUT  TO  THE  TE  TEMPERATURE 

C  EXCESS  ISOTHERM. 

C 

DIMENSION  P< 12,6) 

C0MMON/HAVF/DELP<12,6) ,TEMP(60,6 ),VEL( 60,6),DIR<  60,6) ,XP( 601 ,YPC  60 
1) ,AMTEMP( 6) ,AMCUR (6) , BETAO , BO , ANGN ,T0 ,U0 , NPTS , DR , RD, K 
C 

FT  AT=0 .0 

TC=P(1,K)*(  TO- A  MTE  MP(  K)  ) 

IF(S.LF.P< 2,K)*B0)  GOTO  12 
TC=TC*SORT (P(2,K)*B0/S) 

12  IF(TE.GT.TC)  GOTO  13 

ETAT=1.44269504*AL0G(TC/TE) 

ETAT  =  SQRTIFTAT)MP(3,K)*B0+P<4,K  )  *S  » 

IFIFTAT.LT. 0.0)  F  TAT=0.  0 

13  RETURN 
FND 

£*******  ****************************< c**********************************Q 

c  c 

FUNCTION  ETAV(S,VE,P) 

C  C 

C ********************* ********  ********************************** *******c 

c 

C  FUNCTION  CALCULATES  THE  LATERAL  DISTANCE  OUT  TO  THE  VE  VELOCITY 

C  FXCESS  POINT. 

C 

DIMENSION  P( 12,6) 

COMMON/ HAVE/ DEL P(  12, 6 ) , TFMP ( 60, 6 ) , VFL ( 60, 6 ) ,DI R ( 60 ,6 ) , XP < 60) , YP( 60 
1)  , AMTEMPI6) , AMCUR (6) , BET  AO , BO , ANGN, TO , UO , NPT S , OR , R D, K 

C 

FT AV  =  0 .0 
VC  =  P (7  ,K)*UO 

IF ( S.LF.Pl 8,K)*B0)  GOTO  10 
VC=VC*SQRT ( P ( 8, K )*B0/S) 

10  IF(VE.GT.VC)  GOTO  11 

ETA V=  1.  44269504*AL0G( VC/VE) 

CT  AV  =  SOPT (FT  AV  )  *(  P  (  9  ,  K )*BO+P( 10,K)*S) 

IF (FTAV.LT. 0.0)  FT  AV=0 .0 

11  RETURN 
FND 
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APPENDIX  D 

Fitting  Parameters  and  Results 

This  Appendix  contains  the  final  values  of  the  sets  of  12  parameters 
(six  parameters  for  the  fit  to  the  temperature  measurements  and  six  for  the 
fit  to  the  velocity  data)  obtained  from  the  fitting  procedure  described  in  Sec.  V. 
The  results  for  the  five  jet  studies  analyzed  in  this  manner  are  included. 

Values  are  tabulated  (Tables  D.1-D.5)  for  each  depth  at  which  sufficient  data 
were  available  to  make  the  fitting  procedure  meaningful.  Also  given  are  outfall 
and  ambient  temperatures  and  velocities  used.  In  addition,  the  differences 
between  the  measured  values  at  each  station  and  the  values  calculated  from 
the  final  fitted  functions  are  tabulated.  The  quantity  AT  is  defined  as  the  cal¬ 
culated  temperature  excess  minus  the  measured  temperature  excess  at  each 
data  location;  the  quantity  Au  is  defined  as  the  difference  between  the  calculated 
jet  centerline  component  of  the  local  jet  velocity  for  each  data  point  and  the 
component  of  the  measured  velocity  parallel  to  the  fitted  velocity  centerline. 
Also  included  are  the  root-mean- square  deviations  dp  and  ctu  defined  in  Sec.  V 
and  AT  and  Au,  which  are  the  average  of  the  absolute  values  of  the  tabulated 
deviations.  The  quantities  ATajq  and  Aua^  are  the  average  deviations  for  all 
depths . 

TABLE  D.l.  Results  of  the  Fitting  Procedure  for  May  18,  1972,  Data 

T  =  17. 7°C 
o 

U  =54.7  cm/sec 
o 

=  0.0  cm/sec 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0 

(°C) 

9.2 

8.5 

8.3 

7.7 

7.6 

7.5 

(C°) 

8.5 

9.2 

9.4 

10.0 

10.6 

10.2 

A 

0.84 

0.71 

0.78 

0.68 

0.80 

a 

3.8 

4.8 

3.9 

3.0 

3.3 

c 

0.8 

1.3 

0.4 

2.0 

0.2 

Y 

0.39 

0.35 

0.41 

0.13 

0.19 

Rf 

0.77 

0.75 

0.81 

0.72 

0.61 

Kj. 

1.41 

1.29 

0.85 

2.85 

3.27 

B 

1.47 

1.32 

1.28 

1.24 

0.83 

6 

4.6 

3.7 

3.4 

3.0 

1.4 

D 

0.4 

0.8 

0.5 

0.9 

0.9 

6 

0.20 

0.12 

0.12 

0.06 

0.01 

Ru 

1.04 

0.96 

1.02 

1.03 

0.85 

K 

0.55 

1.24 

0.63 

0.42 

1.98 

A 


u 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

CJrji 

AT 


153 


TABLE  D.l  (Contd.) 

AT  =  Tcalc .  -  Tdata  ^ 

0.5  m  1.0  m  1.5  m  2.0  m  2.5  m  3.0  m 


0.0 

0.5 

— 

-0.7 

— 

-0.6 

-0.6 

-0.4 

-0.5 

-0.2 

— 

1.0 

— 

-0.9 

0.2 

-0.7 

-0.4 

-0.2 

-0.7 

-0.4 

-1.8 

-1.0 

-0.7 

-0.8 

0.0 

0.4 

0.2 

0.7 

0.5 

-0.5 

0.5 

0.5 

1.0 

— 

— 

0.4 

0.5 

-0.3 

-0.8 

0.2 

0.9 

0.0 

0.0 

1.1 

0.1 

0.6 

0.6 

— 

1.1 

— 

-1.4 

0.1 

0.0 

— 

— 

0.7 

0.9 

1.7 

— 

— 

-1.5 

-2.0 

-1.9 

-1.4 

— 

0.2 

0.4 

0.6 

0.3 

— 

-0.2 

0.0 

0.4 

0.5 

— 

1.6 

0.8 

1.0 

1.1 

0.0 

-1.0 

-1.2 

-2.0 

-0.7 

— 

1.2 

1.0 

1.3 

1.1 

0.5 

0.8 

0.6 

0.8 

0.8 

0.2 

=  0.7  C° 


TABLE  D.l  (Contd.) 


Station  # 


AU  =  U 


calc. 


'  Udata  fcm/sec) 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0  m 

1 

5 

3 

-- 

0 

-- 

2 

0 

0 

-2 

-3 

-5 

3 

-- 

-2 

-- 

0 

0 

4 

0 

2 

3 

1 

0 

5 

1 

3 

-1 

-2 

-1 

6 

-5 

-1 

3 

0 

0 

7 

-- 

-- 

-- 

-- 

-- 

8 

1 

1 

0 

-- 

— 

9 

-2 

-4 

-7 

-2 

1 

10 

1 

-4 

4 

10 

0 

11 

6 

3 

-- 

1 

-- 

12 

o 

0 

0 

-- 

-- 

13 

2 

3 

4 

-- 

-- 

14 

-6 

-4 

-5 

-8 

-- 

15 

9 

9 

9 

3 

-- 

16 

-10 

-8 

-11 

-5 

-- 

17 

-3 

0 

-3 

3 

-1 

18 

-8 

1 

0 

0 

-- 

°u 

6 

5 

7 

5 

4 

AU 

4 

i 

! 

3 

AUaii  =  3 

4 

3 

1 
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TABLE  D.2.  Results  of  the  Fitting  Procedure  for  May  23,  1972,  Data 

T  =  21.6  °C 
o 

U  =55.5  cm/sec 
o 

UA  =  0.0  cm/sec 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0  m 

Ta  (°C) 

14.3 

12.6 

12.1 

11.6 

11.3 

10.8 

eo 

7.3 

9.0 

9.5 

10.0 

10.3 

10.8 

A 

0.84 

0.91 

0.87 

0.85 

0.85 

0.96 

a 

5.6 

4.4 

3.9 

3.7 

3.8 

2.7 

C 

1.0 

0.7 

0.7 

0.9 

0.9 

0.5 

X 

0.27 

0.37 

0.32 

0.21 

0.14 

0.20 

R? 

1.28 

1.26 

1.33 

1.22 

1.28 

1.58 

Kj. 

-2.56 

-2.58 

-3.56 

-1.88 

-2.64 

-7.06 

B 

1.04 

0.98 

0.92 

0.93 

0.80 

0.93 

6 

4.8 

4.7 

5.1 

3.6 

2.6 

3.7 

D 

0.8 

0.9 

0.8 

0.8 

0.9 

0.6 

6 

0.17 

0.09 

0.04 

0.03 

0.02 

0.01 

Ru 

1.29 

1.24 

1.25 

1.28 

1.38 

1.49 

*u 

-2.11 

-1.76 

-1.90 

-2.03 

-2.84 

-3.93 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

a™ 


TABLE  D.2  (Contd.) 


AT  “  Tcalc .  '  Tdata  ^ 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0  m 

-0.6 

-2.2 

-2.6 

— 

— 

— 

0.8 

0.8 

1.3 

1.9 

1.6 

— 

-0.3 

0.0 

-0.5 

-0.8 

-0.5 

-0.6 

-0.9 

-0.7 

-1.1 

-1.7 

-1.5 

0.3 

-1.4 

-1.8 

-2.3 

-2.9 

-2.9 

-0.8 

-2.2 

-3.1 

-3.0 

-1.9 

-1.2 

0.3 

1.4 

2.8 

3.4 

3.1 

3.0 

— 

0.4 

0.2 

0.6 

2.3 

— 

— 

1.5 

1.8 

1.0 

— 

— 

— 

-0.3 

0.2 

0.5 

0.6 

0.6 

— 

0.0 

0.6 

1.4 

0.0 

1.2 

— 

-0.1 

0.2 

-0.8 

0.1 

-1.2 

— 

-0.4 

-1.3 

-1.7 

-2.1 

-2.6 

— 

-0.1 

1.5 

0.0 

1.1 

— 

— 

1.0 

-0.1 

0.4 

-0.2 

— 

— 

0.2 

0.2 

0.1 

-0.3 

0.7 

— 

0.1 

-0.2 

0.8 

0.0 

0.8 

— 

0.0 

-0.6 

-1.0 

-2.1 

— 

— 

0.4 

0.5 

0.7 

-0.1 

— 

— 

0.6 

-0.4 

-0.2 

-0.6 

-0.6 

-0.8 

-0.3 

-1.2 

-1.4 

-0.9 

-0.1 

-0.3 

1.4 

1.8 

1.8 

1.4 

1.6 

1.8 

-1.2 

-0.7 

-0.6 

-0.6 

0.1 

-0.6 

0.0 

0.0 

-1.4 

0.0 

0.0 

0.0 

-0.6 

-0.3 

-0.9 

-1.1 

-1.6 

-0.8 

0.2 

-1.1 

-1.5 

0.0 

-1.1 

-1.2 

-0.3 

-0.1 

0.4 

1.0 

1.0 

0.3 

-0.6 

-0.5 

-0.3 

-0.2 

-0.8 

-0.7 

0.2 

0.1 

0.4 

0.4 

0.3 

0.6 

0.9 

1.3 

1.6 

1.5 

1.7 

1.0 

0.6 

0.9 

1.1 

1.0 

1.1 

0.6 

aT^^j  =  0.9  C° 


TABLE  D.2  (Contd.) 


AU  =  Ucalc.  ■  Udata  fcm/sec) 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0  m 

1 

-2 

0 

0 

— 

— 

— 

2 

33 

28 

21 

24 

16 

— 

3 

-4 

0 

2 

-2 

5 

-1 

4 

-12 

-11 

-14 

-12 

3 

2 

5 

-28 

-21 

-25 

-31 

-30 

-6 

6 

-13 

-21 

-19 

-7 

-2 

7 

7 

13 

13 

11 

13 

18 

— 

8  ! 

39 

33 

23 

21 

— 

— 

9 

-2 

-2 

12 

— 

— 

— 

10 

-11 

-16 

2 

6 

-13 

— 

li 

-11 

4 

10 

-4 

-5 

— 

12 

-1 

0 

1 

-1 

-4 

— 

13 

1 

-2 

2 

3 

1 

— 

14 

-6 

-5 

-6 

-6 

— 

— 

15 

-7 

0 

-12 

-1 

— 

— 

16 

2 

2 

-1 

2 

8 

— 

17 

-1 

-8 

-4 

1 

2 

— 

18 

4 

6 

5 

2 

— 

— 

19 

1 

1 

0 

2 

— 

— 

20 

6 

-2 

-5 

-10 

-10 

-6 

21 

2 

0 

7 

4 

7 

0 

22 

0 

-5 

-5 

-1 

3 

0 

23 

-7 

0 

0 

0 

3 

3 

24 

0 

0 

0 

0 

0 

0 

25 

-8 

-8 

-8 

0 

0 

0 

26 

5 

2 

-3 

0 

4 

0 

27 

2 

5 

5 

6 

0 

0 

28 

0 

0 

0 

0 

0 

0 

29 

6 

-2 

0 

0 

0 

0 

°U 

14 

13 

11 

11 

11 

4 

AU 

8 

AUall 

7 

=  6  cm/ sec 

7 

6 

6 

2 
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TABLE  D.3.  Results  of  the  Fitting  Procedure  for  July  13,  1972,  Data 

T  =  20.3  °C 
o 

UQ  =  54.7  cm/sec 
=  5.7  cm/sec 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0  m 

Ta  (°C) 

13.0 

12.2 

11.4 

11.0 

10.6 

10.3 

6o  (C°) 

7.3 

8.1 

8.9 

9.3 

9.7 

10.0 

A 

0.92 

0.89 

0.70 

0.89 

0.81 

a 

6.2 

4.4 

4.6 

4.3 

4.9 

C 

0.3 

0.5 

4.1 

0.7 

0.8 

Rf 

0.97 

0.97 

0.88 

0.94 

0.98 

Kp 

-0.74 

-0.73 

0.07 

0.05 

-0.26 

B 

1.26 

1.22 

1.05 

1.26 

0.95 

8 

7.4 

5.3 

4.8 

7.5 

4.3 

D 

0.1 

0.4 

0.3 

0.7 

0.6 

6 

0.24 

0.18 

0.16 

0.03 

-0.01 

Ru 

0.94 

0.93 

0.98 

0.96 

0.96 

Ku 

-0.09 

-0.11 

-0.18 

-0.24 

-0.07 

iti 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

CJrp 
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TABLE  D. 3  (Contd.) 


AT  '  Tcalc.  -  Tdata  ^ 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

0.1 

0.3 

1.3 

0.8 

0.0 

0.3 

0.0 

0.6 

0.0 

1.4 

-0.6 

-0.6 

-0.9 

-0.3 

-0.6 

0.4 

0.6 

-1.2 

-0.1 

— 

0.0 

-0.1 

0.9 

— 

— 

-0.4 

-0.8 

-0.6 

-2.1 

-1.4 

0.8 

1.0 

0.8 

1.0 

-0.4 

-0.5 

-1.1 

-0.8 

-1.0 

-0.9 

-0.9 

-1.1 

-1.3 

-0.5 

-1.6 

0.8 

0.4 

0.4 

-1.0 

-0.7 

-0.3 

-0.6 

-0.1 

— 

— 

-0.1 

0.2 

0.2 

1.1 

-0.7 

0.1 

0.4 

0.5 

0.8 

1.1 

0.1 

0.6 

1.1 

1.1 

1.5 

0.4 

0.7 

-0.2 

0.0 

-0.7 

-0.2 

0.3 

-0.8 

-0.6 

-0.4 

0.1 

1.1 

0.2 

0.8 

0.2 

0.6 

0.8 

1.0 

1.2 

1.3 

0.4 

0.6 

0.7 

0.8 

0.8 

ATau  =  0.6  C° 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

CTTI 


TABLE  D.3  (Contd.) 


AU  =  Ucalc.  '  Udata  (cm/sec) 


1.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

9 

5 

-2 

5 

1 

-4 

-4 

5 

6 

0 

3 

2 

-7 

-7 

2 

-4 

-1 

-5 

-5 

— 

-6 

-4 

4 

— 

— 

2 

6 

4 

2 

1 

-3 

-5 

0 

-1 

5 

2 

-3 

2 

9 

0 

-9 

-5 

-12 

-13 

-5 

2 

-2 

-3 

0 

-1 

0 

3 

7 

— 

— 

7 

5 

1 

2 

3 

7 

10 

16 

2 

2 

2 

6 

10 

6 

3 

-4 

-2 

-1 

3 

4 

-1 

4 

3 

8 

7 

10 

6 

5 

4 

3 

7 

6 

8 

8 

4 

4 

4 

5 

5 

3 

AUall  =  4  cm/sec 
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TABLE  D.4.  Results  of  the  Fitting  Procedure  for  September  9,  1972,  Data 

T  =  24. 5°C 
o 

U  =54.7  cm/sec 
o 

UA  =  2.2  cm/sec 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0  m 

Ta  (°Q 

16.3 

15.2 

14.4 

14.1 

13.9 

13.8 

e  (C°) 

8.2 

9.3 

10.1 

10.4 

10.6 

10.7 

A 

0.79 

0.68 

0.66 

0.88 

0.31 

a 

4.5 

3.7 

2.1 

2.0 

0.4 

C 

0.1 

1.1 

2.1 

2.0 

0.8 

X 

0.64 

0.58 

0.45 

0.26 

0.61 

Rp 

0.98 

1.07 

1.22 

1.34 

1.21 

Kr 

-0.44 

-1.25 

-2.18 

-3.08 

-0.45 

B 

1.03 

0.85 

1.00 

0.99 

6 

4.4 

3.0 

4.0 

3.9 

D 

0.8 

0.6 

1.2 

1.2 

6 

0.22 

0.20 

-0.02 

-0.06 

Ru 

0.98 

0.98 

0.97 

1.01 

Ku 

0.58 

1.24 

0.82 

0.60 
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TABLE  D.4  (Contd.) 


Station  # 


"  ■  Tcalc .  -  Tdata  ^ 


0.5  m 

1.0  m 

1.5  m 

2.0  m 

2.5  m 

3.0  m 

1 

-0.9 

0.0 

— 

— 

— 

2 

0.6 

0.9 

1.0 

— 

— 

3 

-0.7 

-1.4 

-0.2 

0.6 

— 

4 

-1.2 

0.0 

-1.5 

-1.4 

— 

5 

1.6 

0.7 

2.2 

1.3 

— 

6 

-0.4 

-0.3 

-0.6 

-0.1 

— 

7 

0.2 

-0.2 

-0.8 

-0.9 

0.0 

8 

0.1 

-0.2 

-0.7 

-0.2 

0.2 

9 

0.1 

0.6 

-0.6 

0.3 

-0.5 

10 

0.6 

0.0 

0.2 

0.1 

-0.6 

11 

0.5 

0.1 

-0.4 

0.6 

0.2 

12 

-0.2 

0.3 

0.2 

-0.6 

0.1 

13 

-0.2 

-0.1 

0.2 

-0.4 

0.2 

14 

-0.3 

0.3 

0.6 

0.8 

0.5 

15 

-0.6 

0.0 

1.0 

1.4 

0.4 

16 

-0.4 

0.1 

0.0 

0.8 

0.3 

17 

0.0 

-0.9 

-1.6 

0.3 

-0.4 

0  rp 

0.8 

0.7 

1.2 

1.0 

0.5 

AT 

0.5 

0.4 

0.7 

0.7 

0.3 

"all 

=  0.5  C° 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

0TT 
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TABLE  D. 4  (Contd.) 

AU  =  Ucalc.  ■  Udata 


0.5  m  1.0  m  1.5  m 

6  3 

3  2  3 

-5  -6  -3 

-4  2  2 

12-1 
-4  -1  6 

3  5  0 

3  3  0 

6  2  0 

0  -3  -4 

0  -10  -6 

0  13 

7  4  0 

0  2  0 

-3  -3  1 

-4  -3  1 

10  0 

5  5  4 

3  3  2 


AUau  =  3  cm/sec 


2.0  m  2.5  m  3.0  m 


-4 

2 

-2 

7 

0 

1 

0 

-7 

0 

4 

0 

0 

1 

0 

1 

4 

2 


164 


TABLE  D.5.  Results  of  the  Fitting  Procedure  for  October  10,  1972,  Data 

T  =  22.2  °C 
o 

UQ  =  42.3  cm/ sec 
UA  =  0.0  cm/sec 


7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

CJrp 
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TABLE  D.  5  (Contd.) 

iT  ■  Tcalc .  -  Tdata  ^ 


0.5  in 

1.0  m 

1.5  m 

0.2 

-0.7 

0.5 

0.3 

0.9 

0.8 

-0.1 

-0.6 

-2.1 

1.0 

2.2 

1.3 

0.3 

0.2 

— 

-1.0 

0.0 

— 

0.8 

2.1 

— 

-0.3 

-2.0 

-1.2 

-1.2 

— 

— 

-0.7 

-1.1 

1.0 

0.7 

-0.2 

-0.5 

0.0 

-0.6 

-0.8 

-0.1 

0.9 

-0.2 

0.2 

-1.0 

-0.5 

-0.2 

-0.9 

0.7 

-0.2 

-0.1 

0.7 

-0.4 

0.8 

0.1 

1.0 

0.7 

-0.6 

-0.4 

-0.6 

-0.8 

0.7 

1.3 

1.2 

0.5 

0.9 

0.8 

AT  ,, 
all 

=  0.7  C° 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

aT  i 


TABLE  D.5  (Contd.) 


AU  =  U 


0.5  m 

1.0  m 

1.5 

0 

1 

2 

-2 

3 

2 

1 

0 

2 

9 

2 

3 

1 

5 

— 

5 

4 

— 

7 

6 

— 

-10 

7 

-5 

-2 

/ 

-7 

0 

-1 

-7 

0 

1 

2 

2 

2 

-1 

0 

-1 

-9 

-2 

-2 

2 

2 

-1 

-3 

1 

0 

2 

-2 

0 

1 

-1 

4 

-4 

3 

0 

6 

3 

2 

4 

2 

2 

calc.  '  Udata  (an/sec) 
2.0  m  2.5  m 


3.0  m 


AUan  =  3  an/sec 
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